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ABSTRACT 


This  report  presents  four  geotechnj  1  engineering 
programs  for  use  on  personal  computing  sterns.  An  Apple 
II-Plus  operating  with  DOS  3.3  Applesof  language  was  used. 
The  programs  include  the  solution  of  the  signpost  problem, 
the  cantilevered  sheet  pile  problem,  the  slope  stability 
problem,  and  the  flexible  pavement  design  program. 

Each  chapter  is  independent  and  does  not  rely  upon 
theories  or  data  presented  in  other  chapters.  A  chapter 
outlines  the  theory  used  and  also  presents  a  users  guide, 
a  program  list,  and  verification  of  the  program  by  hand 
calculation. 

This  report  assimilates  the  product  a  practicing  engi¬ 
neer  would  expect  to  receive  when  procuring  software  services 


CHAPTER  I 


INTRODUCTION 

Micro-computers  are  rapidly  becoming  the  work-horse 
of  the  small  business  and  engineering  world.  All  types 
of  businesses  are  finding  the  micro-computer  an  invaluable 
tool.  Uses  range  from  cost  accounting  to  word  processing 
to  an  engineering  calculating  machine.  The  key  advantages 
of  personal  computers  over  previous  methods  are  its  easy 
access,  speed,  reliability,  and  accuracy. 

Engineering  firms  whose  availability  to  main  frame 
computing  facilities  have  been  limited  by  economics  or  demo 
graphics  can  now  acquire  personal  computers  and  software 
for  a  fraction  of  the  capi'  al  outlay.  Software  can  be 
designed  to  fit  the  precise  needs  of  the  firm  whereas  firms 
used  to  cater  their  needs  around  an  established  software 
base.  This  is  especially  advantageous  to  specialized  firms 
whose  services  demand  repetitive  engineering  problems. 
Rather  than  expending  labor  on  iterative  problem  solving, 
a  personal  computer  matched  with  properly  designed  software 
can  now  provide  engineering  solutions  at  a  fraction  of 


the  cost. 


The  introduction  of  the  personal  computer  into  the 
engineering  firm  provides  a  domino  effect.  The  time  which 
was  once  expended  on  iterative  problem  solving  can  now 
be  devoted  to  more  productive  activity  such  as  consideration 
of  more  complex  problems  and  bidding  strategy.  Accordingly, 
the  small  engineering  firm  can  now  bid  more  jobs,  accomplish 
more  work  and  subsequently,  increase  the  firm's  net  worth. 

It  is  becoming  commonplace  for  short  and  long-range  business 
plans  to  include  the  purchase  and  use  of  personal  computers. 

The  personal  computer  is  a  godsend  to  the  geotechnical 
engineer.  Due  to  the  nature  of  soil,  the  engineer  is  con¬ 
stantly  dealing  with  lower  and  upper  bounds  of  possible 
problem  solutions.  Unlike  concrete  and  steel,  soil  has 
variable  engineering  properties  and  cannot  be  relied  upon 
to  perform  in  a  consistent  manner;  consequently,  the  engineer 
must  consider  several  possible  combinations  of  soil  behavior 
in  order  to  provide  a  safe  design.  Once  the  soil  charac¬ 
teristics  have  been  normalized,  the  engineer  uses  mechanics 
of  particulate  matter  to  best  approximate  the  response 
and  behavior  of  soil  acted  upon  by  external  forces  and 
natural  phenomena  such  as  flow  of  water  through  the  medium 
under  a  hydrostatic  head. 

At  this  point,  the  personal  computer  comes  into  play. 

The  computer  will  calculate  quantities  and  values  according 
to  a  predetermined  sequence.  If  the  same  imput  variables 


are  used,  the  computer  will  calculate  identical  values 
and  quantities.  The  engineer  then  varies  the  input  variables 
according  to  his  evaluation  of  the  possible  conditions 
that  may  exist  for  a  particular  problem.  The  output  will 
then  represent  a  range  of  expected  behavior  of  which  the 
engineer  will  use  for  his  design.  Problems  such  as  cal¬ 
culating  the  required  penetration  of  a  cantilevered  wall 
can  take  up  to  four  hours  to  calculate.  This  time  represents 
several  calculation  iterations  with  one  particular  set 
of  soil  data.  If  no  math  errors  were  made,  one  possible 
bounding  answer  would  be  established.  The  same  procedures 
would  be  repeated  to  establish  another  bound.  Two  trial 
boundary  conditions  have  been  established,  but  what  if 
intermediate  values  are  not  linearly  related  to  the  boundary 
values?  The  prudent  engineer  would  make  intermediate  value 
calculations.  The  time  involved  can  be  enormous.  The 
personal  computer  can  calculate  iterative  problems  in  a 
fraction  of  the  time  required  to  hand  calculate  the  problems 
and  without  the  math  errors  associated  with  hand  calculations. 

The  scope  of  this  special  report  is  to  program  several 
iterative  problems  of  interest  to  geotechnical  engineers. 

The  programs  includes  the  calculation  of  the  required  depth 
of  a  vertical  post  subjected  to  lateral  loads  (SIGNPOST  1) , 
the  required  penetration  of  a  cantilevered  wall  (CANTWALL  1) , 
and  the  design  of  flexible  pavement  (AASHTO  1) .  Additionally, 


a  slope  stability  program  was  translated  and  modified  by 
the  author  (BISHOP  1) . 

No  report  of  this  nature  would  be  complete  without 
a  warning  about  the  ignorant  use  of  computing  software. 

Two  major  problems,  separate  or  combined,  can  render  the 
software  and  subsequent  solutions  totally  useless.  The 
user  must  understand  the  problem  that  the  software  is  pre¬ 
sumed  to  solve.  In  general,  there  are  several  methods 
or  algorhythms  that  can  be  used  to  solve  engineering  problems, 
but  each  method  is  best  suited  for  a  particular  variation; 
furthermore,  the  solutions  can  significantly  vary  from 
one  method  to  another.  This  is  particularly  evident  with 
dynamic  pile  driving  formulas.  The  user  must  understand 
the  use  and  limitations  of  the  program  software.  Although 
the  user  may  understand  the  problem  and  its  methods  of 
solution,  an  erroneous  entry  or  a  program  option  inadvertently 
exercised  can  invalidate  the  computed  solution.  The  solution 
should  be  scrutinized  against  past  experience  and  sound 
engineering  judgment.  As  a  final  check,  a  hand  calculation 
of  the  final  solution  should  be  made.  If  the  user  allows 
the  computer  to  perform  the  iterations  and  hand  checks 
the  final  iteration,  the  user  can  be  assured  as  to  the 
validity  of  the  computed  solution. 

It  is  the  author's  intent  to  convey  a  concise  description 
of  the  software's  use  and  limitations.  It  would  behoove  the 


CHAPTER  II 

SIGNPOST  PROBLEM 

2 . 1  Problem  Definition 

The  computer  program  "SIGNPOST  1"  calculates  the  minimum 
safe  embedment  of  a  cantilevered  pole  subjected  to  lateral 
overturning  loads.  Figure  2.1  defines  the  input  variables 
and  general  geometry  of  the  problem.  Plastic  theory  is 
utilized,  thus  considerable  deflections  are  anticipated. 

All  loads  are  laterally  applied  to  the  pole;  overturning 
moments  are  resisted  by  passive  earth  pressure. 


Diameter 

Figure  2.1.  General  Problem  Definition 


2 . 2  Background  Theor 


Development  of  the  laterally  loaded  cantilevered  pole 
problem  was  initiated  by  J.  F.  Seiler  in  cooperation  with 
the  American  Wood-Preservers'  Association  in  the  early 
1930 's  [1].  Seiler  correctly  diagrammed  the  earth  pressures 
about  the  embedded  pole  (Figure  2.2). 


Figure  2.2.  Stress  Distribution  on  Foundation  [3]. 

Seiler's  objective  was  to  classify  embedment  depths  for 
particular  timber  pole  classes,  as  industrial  demand  was 
increasing  and  the  association  felt  an  economic  need  existed 
for  the  correlation  between  pole  class  and  required  embed¬ 
ment.  Timber  poles  are  classified  according  to  the  strength 
in  bending;  Class  1  being  the  strongest,  Class  6  the  weakest 
in  bending. 


Seiler  began  his  research  on  the  premise  that  a  particu¬ 
lar  pole  class  was  best  suited  for  a  particular  soil  type. 

The  most  economic  use  of  timber  coupled  with  decreased 
labor  of  excavation  lead  to  the  conclusion  that  if  a  proper 
pole  was  used,  it  would  develop  its  full  bending  strength 
just  prior  to  soil  failure.  This  conclusion  is  rational 
and  warrants  further  investigation.  Seiler,  like  many 
engineers,  was  unclear  of  the  definition  of  soil  failure. 
Although  he  properly  perceived  the  earth  pressure  diagrams, 
Seiler's  analysis  focused  on  pole  rotation  when  laterally 
loaded  and  ignored  the  soil  pressures  mobilized  by  the 
pole  rotation. 

A  majority  of  the  equations  developed  by  this  premise 
are  contingent  upon  the  angle  of  rotation  the  pole  would 
undergo  when  laterally  loaded.  Although  indirectly,  Seiler 
was  alluding  to  plastic  theory,  but  allowed  his  stress  analysi 
to  go  beyond  the  stress  which  would  cause  plastic  failure. 

His  analysis  never  made  its  debut  in  the  engineering  litera¬ 
ture.  Seiler  adequately  described  the  earth  pressures 
and  respective  depths  about  the  embedded  pole. 

In  the  early  1940's,  Professor  P.  C.  Rutledge  was 
requested  to  devise  a  system  by  which  embedment  depths 
for  signposts  could  be  estimated.  In  association  with 
the  Outdoor  Advertising  Association  of  America,  Professor 
Rutledge  devised  the  nomograph  as  presented  in  Figure  2.3  [2] 
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U  2 


where , 

D  =  embedment  depth  (ft) 

P  =  lateral  load  (#) 

2 

S-^  =  average  passive  soil  pressure  (#/ft  ) 

H  =  height  of  P  above  grade  (ft) 

B  =  diameter  (ft) 

is  dependent  upon  depth  as  it  is  the  passive  soil 
resistance  at  .34D;  thus,  the  above  equation  must  be  iterated 
To  calculate  a  safe  embedment  depth,  the  soil  pressure 
depths  presented  by  Seiler  and  later  used  by  Professor 
Rutledge  were  corroborated  by  Professors  W.  L.  Shilts, 


L.  D.  Graves,  G.  F.  Driscoll  of  Notre  Dame  University  and 
by  Dr.  J.  0.  Osterberg  of  Northwest  University  [2] . 
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Due  to  the  limited  ability  to  test  soils,  and  the 
lack  of  standardized  soil  classification,  Rutledge  devised 
a  testing  device  which  could  be  used  to  determine  the 
in-situ  average  soil  pressure  (S^)  [2] .  It  consisted  of 

a  l-Js"  hand  auger  which  after  being  rotated  into  the  soil 
would  be  pulled  up.  The  force  required  to  pull  out  the 
auger  was  correlated  to  S-^  (Figure  2.3).  A  scale  for 
cohesionless  soils  and  a  scale  for  silts  and  clays  are  pro¬ 
vided.  The  nomograph  is  limited  to  embedment  depths  of 
10',  post  diameters  of  6"  to  24"  and  a  load  height  of  24'. 

The  above  equation  must  be  used  for  any  parameters  beyond 
these  boundaries. 

In  1957,  D.  Patterson,  being  dissatisfied  with  Rutledges 
soil  test  method,  modified  the  nomograph  to  include  five 
general  soil  type  categories;  i.e.,  very  soft,  poor  soil, 
average  soil,  good  soil,  and  very  hard  soil  [2].  To  augment 
this  general  classification,  the  following  table  was  also 
provided ; 


Table  2.1.  Generalized  Soil  Classifications 


Clay,  in  lumps,  dry 

Poor  soil 

Clay,  damp,  plastic 

Poor  soil 

Clay  and  gravel,  dry 

Average  soil 

Clay,  gravel  and  sand,  dry 

Average  soil 

Earth,  loose,  perfectly  dry 

Average  soil 

Earth,  packed,  perfectly  dry 

Average  soil 

Earth,  loose,  slightly  moist 

Average  soil 

Earth,  packed,  more  moist 

Very  hard  soil 

Earth,  soft  flowing  mud 

Very  soft  soil 

Earth,  soft  mud,  packed 

Poor  soil 

Gravel,  one  inch  and  under,  dry 
Gravel,  two  and  one-half  inches 

Good  soil 

and  under,  dry 

Average  soil 

Sand,  clean  and  dry 

Average  soil 

Sand,  river,  dry 

Average  soil 

Patterson  contended  that  in  the  absence  of  better 
soil  data  the  above  table  would  yield  satisfactory  embed¬ 
ment  depths. 

In  an  effort  to  refine  the  soil  data  input,  D.  L. 

Ivey  and  L.  Hawkins  [3]  applied  Rankine's  formula  for  passive 
soil  resistance: 


Pp 


yzN 

<P 


+ 


2c 


N  =  tan2  (45°  +  * / 2 ) 
<P 


With  this  formula,  can  be  calculated  using  soil  strength 


data,  C  and  <p ;  furthermore,  introduction  of  a  ground  water 
table  with  subsequent  bouyant  forces  can  be  accounted  for. 


Ivey  and  Hawkins  extended  the  design  process  to  include 
checking  the  lower  stresses  (S9  &  SL)  against  the  allowable 
stresses  calculated  by  Rankine  formula.  This  is  especially 
critical  when  the  ground  water  table  is  at  or  near  0.68D 
as  the  lower  allowable  stresses  will  be  reduced.  Ivey 
and  Hawkins  recommended  applying  a  safety  factor  to  the 
design  by  dividing  the  ultimate  passive  resistance  by  a 
factor  of  safety  prior  to  checking  the  working  stresses 
(SI,  S2,  and  SL) . 

2 . 3  Programming  Rationale 

The  ultimate  purpose  of  "SIGNPOST  1"  is  to  provide 
the  user  a  design  depth  capable  of  resisting  a  specified 
lateral  load.  The  "Rutledge"  method  is  used  as  modified 
by  Ivey  and  Hawkins  as  described  in  Section  2.2.  The  program 
is  user  oriented  in  that  all  input  is  prompted  by  clear, 
concise  questions  displayed  on  the  monitor  prior  to  the 
appropriate  input  command.  Output  format  is  intended  to 
provide  easy  identification  of  the  calculated  solutions 
and  supplemented  by  intermediate  values  calculated  during 
the  program  routine.  The  flow  chart,  Figure  2.4,  is  provided 
as  a  skeleton  of  basic  routines  and  conditionals  as  well 
as  the  general  sequence  of  events  performed  between  the 
input  of  the  problem  parameters  and  the  printed  output. 


Program  Use  and  Limitations 


2.4.1  General 


"SIGNPOST  1"  was  programmed  on  an  Apple  II-Plus  with 
64K  random  access  memory.  The  disk  operating  system  was 
version  3.3  (DOS  3.3).  Prior  to  using  this  program,  the 
user  should  be  familiar  with  the  system  control  features 


of  the  Apple  II-Plus.  Namely,  the  user  needs  to  know  how  j 

to  LOAD,  RUN  and  use  the  return  key;  all  other  commands 
and  options  are  integrated  into  the  program. 

As  with  any  computer  run,  the  user  should  have  reviewed  I 

the  required  input  parameters  and  have  them  available  (prefer¬ 
ably  in  order  of  response)  prior  to  running  the  program. 

If  more  than  one  set  of  input  parameters  are  to  be  used,  I 

it  is  recommended  that  a  table  be  set  up  to  prevent  erroneous 
input  and  reduce  input  time. 

2.4.2  Input  I 

The  following  is  a  list  of  input  parameters  for 
SIGNPOST  1: 


a) 

b) 

c) 

d) 

e) 

f) 

g) 

h) 

i) 

j) 

k) 


number  of  runs  | 

angle  of  internal  friction  (degrees) 

cohesion  (PSF) 

wet  soil  weight  (PCF) 

saturated  soil  weight  (PCF)  ’ 

wind  pressure  (PSF)  (see  Fig.  2.5) 

height  to  sign  centroid  (ft)  | 

sign  area  (SF) 
post  load  (#) 

ground  water  table  depth  (ft) 
number  of  diameters 

I 


i 


RECOMMENDED  MINIMUM  DESIGN  WIND  PRESSURES 

(  Pounds  per  Sipiare  Foot) 

lased  upon  data  obtained  from  U.S.  WE  ATI  I  Ell  BUREAU  and 
AMERICAN  STANDABD  BUILDING  BEQUIBEMENTS 


Figure  2.5.  Recommended  Minimum  Design  Wind  Pressures [2] 


k]_)  diameter  #1  (ft) 
k2>  diameter  #2  (ft) 
kjj  diameter  #i  (ft) 

l)  safety  factor 

m)  depth  calculation  tolerance 

n)  passive  coefficient  of  earth  pressure* 

*Note:  This  quantity  can  be  input  or  the  program  can  calcu¬ 
late  the  passive  coefficient  of  earth  pressure  at  the  option 
of  the  user. 

2.4.3  Apple  II  Start  Up 

A  brief  discussion  of  the  steps  required  to  run  a 
program  on  the  Apple  II-Plus  will  be  provided  here. 

a)  Plug  the  computer  and  the  monitor  into  any  llOv, 

60  Hz,  power  outlet. 

b)  Turn  the  computer  on  by  flipping  the  switch  located 
on  left  back  side  to  the  on  position.  Hit  the  "reset" 

key  to  stop  the  disc  drive  from  turning. 

c)  Turn  the  monitor  on  by  pulling  out  the  brightness 
control . 

d)  Insert  the  system  master  diskette  into  disk  drive 

#1. 

e)  Type  PR# 6. 

f)  After  the  disk  drive  stops  (light  off)  and  the 
blinking  cursor  appears,  remove  the  system  master  diskette. 

g)  Insert  the  slave  diskette  with  the  program  to  be 
run  into  the  disk  drive. 

If  the  user  is  unsure  about  this  procedure,  refer 


to  the  operation  manuals  provided  with  tre  computer. 


If  a  printer  is  attached  to  the  computer  type  RUN 
SIGNPOST  1.  Include  a  space  between  RUN  and  SIGNPOST  and 
a  space  between  SIGNPOST  and  1.  Type  the  return  key.  If 


a  printer  is  not  available,  the  user  must  type  LOAD  SIGNPOST 
1.  Once  the  flashing  cursor  appears,  type  DEL  1080,  1080. 
This  will  delete  a  DOS  command  which  switches  the  output 
to  the  printer.  Type  RUN. 

2.4.4  Program  Use 

An  introduction  will  begin;  the  last  line  will  be 
the  first  user  question.  See  Section  2.7  for  listing  of 
all  introductory  statements,  user  prompting  questions, 
user  option  questions,  and  output. 

Each  run  is  associated  with  a  particular  set  of  input 
data.  Within  each  run  the  user  may  specify  up  to  a  maximum 
of  ten  diameters  for  which  a  required  embedment  depth  will 
be  calculated.  Refer  to  Section  2.7  for  an  example.  This 
allows  the  user  to  observe  the  change  in  required  depth 
associated  with  a  change  in  post  diameter  without  having 
to  repetitiously  input  the  bulk  of  the  input  data. 

SIGNPOST  1  allows  the  user  to  input  an  additional 
load  on  the  post  (see  Figure  2.1).  Rutledge's  depth  equation 
was  modified  to  include  an  additional  load  and  moment  arm. 


D  = 


2.37(P+P1)  + 


\T( 2~. 


37 (P+P1 ) 


+  4  x  2.64  (PH+P-^-^B 


2S1B 


A  practical  use  of  this  feature  involves  considering  wind 
pressure  against  wide  portions  of  the  cantilevered  post. 

The  user  must  input  force  and  height  (Figure  2.1).  If  no 
force  is  to  be  included,  the  user  should  input  zero  for 
both  force  and  height. 

SIGNPOST  1  provides  the  user  the  option  to  input  the 
coefficient  of  passive  earth  pressure  or  the  program  can 
calculate  the  coefficient.  Rankine's  formula  is  used. 

K  =  N  =  tan2  (45  -  4> / 2 ) 

P  <P 

The  user  may  know  from  experience  the  value  of  the 
coefficient  and  may  want  to  input  it  rather  than  have 
it  calculated.  The  computer  will  ask  the  user  which  option 
the  user  desires.  If  the  user  responds  to  input  the  coef¬ 
ficient,  the  computer  will  respond  with  the  question,  "What 
is  KP?"  The  user  should  at  that  time  input  the  desired 
coefficient  of  passive  earth  pressure. 

All  data  must  be  input  in  units  of  degrees,  pounds, 
and  feet.  In  order  to  change  the  units,  the  unit  weight 
of  water  must  be  changed  in  line  5030  from  62.4  PCF  to 
whatever  units  are  to  be  input.  The  output  will  be  in 
terms  of  the  new  unit,  but  the  print  statement  will  still 
print  out  degrees,  PCF,  etc.  behind  the  variables  and  solu- 


Refer  to  Figure  2.2  for  interpretation  of  the  output 
data.  S  is  the  maximum  pressure  of  the  parabolic  stress 
distribution  on  the  upper  2/3  of  the  embedded  post.  SI 
is  the  average  pressure  of  the  parabolic  stress  distribution. 
S2  is  calculated  from  SI  and  represents  an  average  stress 
on  the  lower  1/3  of  the  embedded  post.  S2  is  compared 
to  S2  ALLOW  to  insure  that  it  is  less.  S2  ALLOW  is  calculated 
from  Rankine's  formula  of  earth  pressure.  SL  is  the  maximum 
pressure  mobilized  at  the  bottom  of  the  embedded  post. 

SL  is  a  function  of  S2,  thus  a  function  of  SI.  SL  ALLOW 
is  calculated  by  Rankine 1 s  formula  of  earth  pressure  but 
unlike  S2  ALLOW,  the  safety  factor  is  taken  to  be  one.  Ivey 
and  Hawkins  [3]  contends  that  due  to  local  plastic  failure 
at  the  butt  of  the  post,  the  stress  would  distribute  upward; 
thus,  the  ultimate  pressure  versus  the  safe  pressure  should 
be  used  for  comparison. 

Following  the  output  of  the  calculated  data,  a  list 
of  the  input  data  will  follow.  This  list  will  serve  as 
a  verification  of  proper  data  input  as  well  as  a  record 
of  data  used  to  produce  the  calculated  output. 
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“l 

SPEED® 

150 

10 

PRINT 

II 

4.  4- 4+ 4- -*■*■■*»  4- 4- 

11 

PRINT 

II 

♦SIGNPOST* 

12 

PRINT 

II 

*•*•*■*■  ■*■  *■*■*■** 

13 

PRINT 

14 

PRINT 

15 

PRINT 

20  PRINT  "DANA  K.  EDDY,  578-80-8373 " 

22  PRINT  “GO.  INSTITUTE  OF  TECHNOLOGY " 

24  PRINT  "SCHOOL  OF  CIUIL  ENGINEERING" 

26  PRINT  "DEPARTMENT  OF  GEOTECHNICAL  ENGINEERING" 

30  PRINT 

31  PRINT 

32  PRINT 

35  PRINT  "SYSTEM  HARDWARE:  APPLE  II  PLUS  <84K>" 

37  PRINT  "SYSTEM  HARDWARE:  DOS  3.3,  APPLESOFT  BASIC  LANGUAGE" 

39  PRINT  "PROGRAM  DATE:  APRIL,  1983" 

42  PRINT 

43  PRINT 

44  PRINT 

45  PRINT 

50  PRINT  "SIGNPOST  ESTIMATES  THE  MINIMUM  EMBEDMENT  DEPTH  OF  A  SINGLE  CANT 

I LEUERED  POST  FOUNDATION.  THE  CLASSICAL  APPLICATION  IS  A  SIGN  OR  MAR 
DUE  SUBJECTED  TO  WIND  LOADS." 

51  PRINT 
32  PRINT 
53  PRINT 

i 00  3PEED=  255 

320  PRINT  "HOW  MANY  PROBLEM  SETS  DO  YOU  WANT  TO  RUN?  THE  USER  MAY  INPUT 
SEUERAL  POST  DIAMETERS  PER  PROBLEM  SET." 

330  INPUT  N 

331  PRINT 

332  PRINT 

340  DIM  A<N,28> 

350  FOR  I  =  1  TO  N 

3S0  PRINT  "WHAT  IS  THE  ANGLE  OF  INTERNAL  FRICTION  OF  THE  SOIL  FOR  RUN  *"I 
"?  C DEGREES)" 

370  INPUT  AC  1,1 > 

330  PRINT  "WHAT  IS  THE  COHESION  OF  THE  SOIL  FOR  RUN  #"I"?  C  PSF 
390  INPUT  A-:  1,2!:' 

■+00  PRINT  "WHAT  IS  THE  WET  WEIGHT  OF  THE  SOIL  FOR  RUN  #"I"?  CPCFV 
410  INPUT  AC  1,3 > 

+20  PRINT  "WHAT  IS  THE  SATURATED  WEIGHT  OF  THE  SOIL  FOR  RUN  #" I "?  CPCFV1 
430  INPUT  AC  1,4) 

440  PRINT  "WHAT  IS  THE  WIND  PRESSURE  AGAINST  THE  SIGN  FOR  RUN  #"I"?  C  PSF  > 

II 

+50  INPUT  AC  1,5) 
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460  PRINT  "WHAT  IS  THE  HEIGHT  ABOVE  GRADE  TO  THE  CENTROID  OF  THE  SIGH  FOR 
RUN  #"I"?  C  FEET  > " 

470  INPUT  ACI#6> 

430  PRINT  "WHAT  IS  THE  AREA  OF  THE  SIGN  FOR  RUN  #"I"?  CSF>" 

496  INPUT  AC  1,7 > 

500  PRINT  "WHAT  IS  THE  LOAD  ON  THE  POST  FOR  RUN  #"I"?  <  POUNDS  >" 

510  INPUT  AC  1,8 > 

520  PRINT  "WHAT  IS  THE  HEIGHT  ABOUE  GRADE  OF  THE  POST  LOAD  FOR  RUN  #"  I 7 
C  FEET  >" 

530  INPUT  AC  I ,3> 

540  PRINT  "WHAT  IS  THE  DEPTH  BELOW  GRADE  OF  THE  GROUND  WATER  TABLE  FOR  RU 
N  #"I"?  C  FEET  .V 
550  INPUT  AC  1,10 > 

5S0  PRINT  "HOW  MANY  POST  HOLE  DIAMETERS  DO  YOU  WANT  TO  INPUT  FOR  RUN  4"!" 

*7* 11 

570  INPUT  AC  1*11) 

5S0  FOR  J  —  1  TO  AC  1*11!' 

530  PRINT  "WHAT  IS  DIAMETER  #"J"  C FEET )" 

600  J  =  J  +  12 
610  INPUT  AC  I , J  > 

620  J  =  J  -  12 
630  NEXT  J 

640  PRINT  "WHAT  IS  THE  SAFETY  FACTOR  FOR  RUN  #"1"?" 

650  INPUT  AC  I  j.  12  > 

651  PRINT  "INPUT  THE  TOLERANCE  FOR  DEPTH  CALCULATION.  C RECOMMEND  .5  -*1. 
0  FT  )" 

652  INPUT  TL 

660  PRINT  "DO  YOU  WANT  TO  INPUT  THE  COEFFICIENT  OF  PASSIVE  EARTH  PRESSURE 

CKP)  # C  YES  > #  OR  HAVE  SIGNPOST  CALCULATE  KP  FOR  YOU  CNO>  ?" 

S70  INPUT  6$:X  =  ASC  CB$>:  IF  X  <  84  GOTO  700 
680  PRINT  "WHAT  IS  KP?" 

690  INPUT  AC  1*23) 

695  GOTO  710 

700  AC  I  #23  >  =  C  TAN  C  C  45  +  AC  I  »  1  )  /  2  >  *  .01745  )  >  2 

7 1 0  AC  I  #  1 1  >  =  AC  I  #  1 1  >  +  12 
720  FOR  K  =  13  TO  AC  I # 1 1 > 

730  AC  I #21 >  =  20 
740  Z  =  .34  *  AC  I  #21!) 

750  GOSUB  5000 

760  AC  I #22  >  =  PS 

770  AC  I #23 )  =  .6667  *  AC I#22> 

730  GOSUB  6080 
730  IF  D  <  AC  I #21 >  GOTO  820 
300  AC  1,21 >  =  AC  I #21 >  +  TL 
310  GOTO  740 

320  IF  D  >  AC  I #21^  -  TL  GOTO  850 
830  AC  I #21 >  =  AC  I #21 >  -  TL 
331  Z  =  .  34  *■  AC  I  #21') 

832  GOSUB  5000: AC  I #22 >  =  PS 
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333  fi<  I #23 >  =  .6667  *  fi<I,22>:  60SUB  6630 
834  GOTO  820 

■-'30  fi'.  I  >24  )  =  fit.  I  >23  )  /  <  <  .  28  *  fi<  I  >21  )  /  <  fid  >6  )  +  .  34  fi<  I  >21  )  )  >  +• 
860  Z  =  .68  *  fid, 21) 

870  G0SU8  5000 
888. fi< 1,25)  =  PS 

680  IF  ft<  I  >25  )  >  fid,  24)  GOTO  968 
808  fid >21  >  =  fid, 21)  +  TL 
918  Z  =  .34  *  fid, 21  > 

928  GOSUB  5080 

930  fi<  I  >  22  =  PS 

948  fid, 23 >  =  .6667  *  fid  ,22) 

950  GOTO  850 

960  fid  ,26)  =  2  *  fid ,24 > 

970  Z  =  fid, 21  > 

980  GOSUB  5080 

998  fid, 27 >  =  PS  *  fid,  12) 

1800  IF  fid ,27)  >  fid ,26 )  GOTO  1080 
1010  fid ,21)  =  fid, 21)  +  TL 
1020  Z  =  .34  *  fid, 21) 

1030  GOSUB  5000 
1040  fid, 22)  =  PS 


1650  fid ,23 >  *  .6667  *  fid, 22) 

1000  fid, 24)  -  fid, 23)  /  <<.28  *  fid ,21)  /  <fi<I»6)  +  .34  *  fid, 21) 
1070  GOTO  960 

1080  L$  =  CHR$  < 4):  PRINT  L*i"PR#r*:  FOR  L  =  1  TO  N 

1081  K  =  K  -  12 

1 088  PRINT 

1089  PRINT 

1890  PRINT  "OUTPUT  FOR  OlfiMETER  #"K"  ,RUN  #"1"." 

1091  PRINT 

1092  PRINT 

1093  PRINT 
1100  K  =  K  +  12 

1110  PRINT  "OlfiMETER  ="fi<I,K>"  FEET  ", "DEPTH  ="fi<I,2D"  FEET" 

1111  PRINT 

1120  PRINT  "S  ="fi< 1,22)"  PSF  ","S1  *"fl<I,23>"  PSF" 

1121  PRINT 

1130  PRINT  "S2  ="fi< I ,24 >"  PSF","S2  fiLLOH  ="fi< I ,25 )"  PSF" 

1131  PRINT 

1 1 40  PRINT  "SL  ="fi<  1,26  )"  PSF" , "SL  fiLLOH  ="fi<I,27)"  PSF  < IJLTIMfiTE 

1141  PRINT 

1142  PRINT 
1150  NEXT  K 

1158  PRINT  "**********************" 

1159  PRINT  " ********************** « 

1 160  PRINT  "INPUT  FOR  RUN  #"I"." 

1161  PRINT  " ********************** •• 

1162  PRINT 
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1170  PRINT  "ANGLE  OF  INTERNAL  FRICTION  ="ACI#1>"  DEGREES" 

1171  PRINT 

1 1 80  PRINT  " CONES I ON  ="A< I #2  > "  PSF " 

1181  PRINT 

1190  PRINT  "KP  ="ACI,28> 

1191  PRINT 

1200*  PRINT  "WET  WEIGHT  OF  SOIL  ="ACI#3>"  PCF" 

1201  PRINT 

1210  PRINT  "SATURATED  WEIGHT  OF  SOIL  ="A(I,4)"  PCF" 

1211  PRINT 

1220  PRINT  "WIND  PRESSURE  =  "ACI#5>"  PSF" 

1221  PRINT 

1230  PRINT  "HEIGHT  OF  SIGN  CENTROID  =  “A<I,6>"  FEET" 

1231  PRINT 

1240  PRINT  "AREA  OF  SIGN  ="ACI#?>"  SF" 

1241  PRINT 

1250  PRINT  "LOAD  ON  POST  =  "AC  I  #8  >"  #“ 

1251  PRINT 

1280  PRINT  "HEIGHT  OF  POST  LOAD  ="A<I,9>"  FEET" 

1261  PRINT 

1270  PRINT  "DEPTH  OF  GWT  =  "A< I # 10 >"  FEET" 

1271  PRINT 

1 280  PR I NT  " SAFETY  FACTOR  ="A< I # 12 ) 

1281  PRINT 

1284  PRINT  "TOLERANCE  =+/-"TL"  FEET" 

1285  PRINT  L$; "PR#0" 

1286  PRINT 

1287  PRINT 
1290  NEXT  I 
1300  END 

5000  IF  2  >  AC  I,10>  GOTO  5030 

f.010  PS  =  CCAC  1,3)  *  Z  *  A<  I  #28  >  >  +  C2  *  A<  I  #2  >  *  AC  I  #23)  A  .  5  >  >  AC  I  .12 
5020  RETURN 

5030  PS  =  <  <  AC  I  ,3)  *  AC  1, 10)  *  A<  I  #28  >  >  +  <2  *  A<  I  #2  >  *  AC  I  #23  >  .  5>  f  < 

A<  I #4  >  -  62. 4  *  <Z  -  AC  I  #10))  *  AC  I#28>))  -y  AC  1/12  > 

5040  RETURN 

£.000  D  =  1.13  *  CCAC  1,5)  *  AC  I  #7))  +  AC  I  #8  >  >  /  ACI#K)  /  AC  I  #23  >  +  CCi.iS 
C  <  AC  1 ,  5  >  *  AC  1 , 7  >  )  +  AC  I  # 8  >  >  ✓  ACI#K)  AC  I  #23))  2  +  <  AC  1,5  >  *  AC  I 

7)  *  AC  I  #6 )  +  AC  I  #3)  *  AC  I  #9))  *  2.63  ✓  AC  I#K>  ✓  AC  I #23) >  .5 

601O  RETURN 


2.6  Variable  List  (SIGNPOST  1) 


Input 


N 

= 

#  of  runs 

A  (1 , 1 ) 

= 

Phi  angle 

A  (1 , 2 ) 

= 

Cohesion 

A  (1 , 3 ) 

= 

Wet  soil  weight 

A  ( 1 , 4  ) 

= 

Saturated  soil  weight 

A  (1 , 5) 

= 

Wind  pressure 

A  (1 , 6 ) 

= 

Height  of  sign  centroid 

A  (1 , 7) 

= 

Sign  area 

A  ( 1 , 8 ) 

= 

Post  load 

A  ( 1 , 9 ) 

= 

Load  height 

A ( I , 10 ) 

= 

GWT  depth 

A  ( 1 , 11 ) 

= 

#  of  diameters 

A ( I , 12 ) 

= 

Safety  factor 

A (I , J) 

= 

Diameter 

A (I / 28) 

= 

Passive  earth  pressure  coefficient 

TL 

= 

Tolerance 

Flow  Control 

X  =  Question  input 

Counters 

I  =  Run  # 

J  =  Diameter  # 


K  =  Output 


Miscellaneous 


Z  =  Depth 

PS  =  Earth  Pressure 

A  ( 1 , 22 )  =  S 

A  ( 1 , 2  3 )  =  SI 

A  ( 1 , 2  4 )  =  S2 

A  (1 , 25)  =  S2  allowable 

A  ( 1 , 26 )  =  SL 

A(I,27)  =  SL  allowable  (ultimate) 
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2 . 7  Program  Verification 


*3 1 GNFuST  i 


DANA  K.  EDOV,  578-89-8378 
GA.  INSTITUTE  OF  TECHNOLOGY 
SCHOOL  OF  CIUIL  ENGINEERING 
DEPARTMENT  OF  GEOTECHNICAL  ENGINEERING 
OR.  RICHARD  0.  BARKSDALE,  ADUISOF; 


5 VST EM  HARDWARE:  APPLE  II  PLUS  <  6410 
S  VST  EM  HARDWARE:  DOS  3.3,  APPLESOFT  BASIC  LANGUAGE- 
PROGRAM  DATE:  APRIL,  1383 


SIGNPOST  ESTIMATES  THE  MINIMUM  EMBEDMENT  DEPTH  OF  A  SINGLE  CANT I 
NDATION.  THE  CLASSICAL  APPLICATION  IS  A  SIGN  OR  MARQUE  SUBJECTE! 


HON  MANY  PROBLEM  SETS  DO  VOU  WANT  TO  RUM?  THE  USER  MAV  INPUT  SE 
TERS  PER  PROBLEM  SET. 


THAT  IS  THE  ANGLE  OF  INTERNAL  FRICTION  OF  THE  SOIL  FOR  RUN  *1 •  ’ 

WnAT  IS  THE  COHESION  OF  THE  SOIL  FOR  RUN  #1?  ■:  PSF  '■ 

•"  SOU 

mi-iAT  IS  THE  WET  WEIGHT  OF  THE  SOIL  FOR  RUN  #  1  7  PCF  • 

7  1  ].  0 


mHAT  IS 
7 1 22 .  4 
mm  AT  IS 
•  4U 

<.-hT  IS 
7 1 07 . 5 
hmhT  IS 
7  4 SO 
mm A 7  IS 

7  i  000 

mm AT  IS 
7  50 

■  -AT  IS 


THE  SATURATED  WEIGHT  OF  THE  SOIL  FOR  RUN  81 7  1  PCF  • 

THE  WIND  PRESSURE  AoA  I  NsT  THE  S I  ON  FuR  RUN  ffl  ■  1  F'sF  ■ 

THE  HEIGHT  ABODE  GRADE  TO  THE  CENTROID  OF  THE  >IGN  F i.'m  - 
I  hit  At.tA  ijF  i  Hfc.  IuN  FuR  RUN  **  1  1  1 

THE  LOAD  ON  THE  POST  FuR  RUN  si  ■  ’  PuUNu'r-  ■' 

THE  HEIGHT  AEiOUE  GRADE  uF  THE  Fl'-T  LUi-iU  FuR  RUN  ttl  ■  •  t- ti 
THE  DEPTH  BELOW  GRADE  OF  THE  GROUND  HATER  TABLE  FOR  RIJN 


LEUERtD  PO 
D  TO  in  ml 


UERAL  PO  :■ 


iEl-R 


» 1 


HOW  MANV  POST  HOLE  DIAMETERS  DO  VOIJ  WANT  TO  INPUT  FOR  RUN  #1? 

WHAT  IS  DIAMETER  #1  <  FEET '• 

7  2 . 0 

WHAT  IS  DIAMETER  #2  FEET  > 

WHAT  IS  DIAMETER  #3  < FEET > 

•?  3 . 0 

Wi-tAT  IS  THE  SAFETV  FACTOR  FOR  RUN  #1? 

INPUT  THE  TOLERANCE  FOR  DEPTH  CALCULATION.  <  RECOMMEND  .5  -  1.0  FT 
?.  5 

DO  VOU  WANT  TO  INPUT  THE  COEFFICIENT  OF  PASSIVE  EARTH  PRESSURE  <  KF 
HUE  SIGNPOST  CALCULATE  KP  FOR  VOU  «'  NO  ? 

'"'NO 


OUTPUT  FOR  DIAMETER  #1  ,RUN  #1. 


DIAMETER  =2  FEET 
S  =2654.32639  PSF 
32  =2932.34826  PSF 
SL  =5364.63651  PSF 


DEPTH  =45. 5  FEET 

51  =1769.37275  PSF 

52  ALLOW  =3339.17206  PSF 

SL  ALLOW  =10295.936  PSF  < ULTIMATE 


JTFTJT  FOR  DIAMETER  #2 


* RUN  #1. 


DIAMETER  =2.5  FEET 
S  =2430.34245  PSF 
i'2  =2776.  661 16  Ps.F 


SL  =c 


DEPTH  =41.5  FEET 

51  =1653.64431  PSF 

52  ALLOW  =3713.35304  PSF 
SL  ALLOW  =9631 .  r  6239  FsF  1  UL  i  IMA  1 1- 


•1 

3 


FSF 


s 


S.6W  A £EA 
^  450  fiz 


107.6 


--  40  Psr 


P,  =  I  ooo 


FPQ'D  bEPfH 


&  =  liopcp 

-  izz.4  pep 
^  --  ZS° 
a  ^  Zoo  psf 


//9a ;  *■  f45"  a  %)  -  taW  z‘^5" a 


-  2.77 


S.7.  =Z.O 


r<0/e  2  a  15" 


(  p)  ^  Z  C  /z 


/TJAt  Z  -  I  S’  ' 


(W)^  -y^o)  (£-6v'0)kf>  +2c  /k. 


Assui/ue.  d  «  41. S’ 


CfiLducATE:  5  €>  S  ~  .  34  D  =  .34(4  1-5 
S  =  *  (wo  ( i4. n')  2..H+  5  (zoo)  (TTis)  / 

=  54ft 


Cfi  LtULATE.  bi 

>L  -  73  s  -  73  iz<i%i.sr) 

-  /6>sc.0 

CAlCU  Cflre  t> 

i. ia(p+p.)  J(  MsCprP.Vi2-.  ft7HH7u.')'W. 


■  P  --  456 ffc1  -  40 1  /Pt2-  -  Z&.ooo.o 

1,1$  (/ooot  16,000) 

U-s~)  i&ss  i~ 


j  ( /./&  (/00C>^/^000)\Z  (flSOOO- /Ol-s)  -i-^000  ■  go) 
v2-£)  /&&S-  j  . 4)  /(.SVT 


5".  4/<9  ^  J25363  ^  fZiet.lL 


4/.3S  ‘  ~  41.5  7 


CffeCK  <  SU 


,  2%D 
H-  f  ,3 45 


.  2 &(4/.3S") 
/01S+  .14(11-35') 


■  SIS 


5  z  "  s  1  /.  S'!  5" 

-  ^ 

-  Si  SO -4 

SZ.  -  S-f?  &  ’Z-  -  >  (oS>V  -  2 &./ , 

^z/um  =[(i\o(is)  +(\z iA-biA)(iz  )1  -15)) 

fsTiJ  /  2. 

-  3108.4  *- 

S  2  ^  £2  /\  LL-OtJ 


sl  -  z(sz)  =  a  (z-7%0.4) 
-  SSL  0.<?)  * }f^ 


^^-ni_\.trvO  -  -ftjo  (ultimate.  G  2  =  b  - 

211)  Lund  '  faoOs)*(n2.4  ^  ?.4)^/.3S“-/5-)) 
a.  (zoo)  iz.11 

-  4b  /S .b  * /g>  ~ 

S>  L-  ^  Si  L  A  llq  tO 

R.e.Q\jigEr>  ~DzprH-  -  41. 35  ' 
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CHAPTER  III 


CANTILEVERED  WALL 

3 . 1  Problem  Definition 

CANTWALL  1  calculates  the  required  embedment  depth 
of  a  cantilevered  wall.  Although  cantilevered  wall  heights 
are  limited  by  structural  constraints  due  to  high  bending 
moments  in  the  wall,  CANTWALL  1  can  calculate  the  theoretical 
penetration  depth  required  to  support  any  height  of  wall. 

The  limitations  of  wall  height  are  discussed  in  more  detail 
in  Section  3.4.  As  depicted  in  Figure  3.1,  the  vertical 
wall  penetrates  through  two  soils.  The  soil  characteristics 
are  specified  in  terms  of  angle  of  internal  friction,  co¬ 
hesion,  saturated,  and  wet  unit  weight.  A  friction  angle 
for  the  soil-wall  interface  must  also  be  specified.  The 
ground  water  table  can  be  specified  to  exist  anywhere  from 
the  top  of  Soil  #1  to  any  depth  below  grade. 

CANTWALL  1  satisfies  the  static  summation  of  horizontal 
forces  and  moments. 


y  (Gl)  MAI) 
Yg  (G2 )  ka(Kl) 

c  (G2 )  6 (B1 ) 


Soil  #1 


Wall  Height 


YH20  =  62.4  PCF  (G5 ) 


Penetration  Depth 


yw(g3> 

Ys(G4) 

c  (C2 ) 
kp(K2) 

(k3) 

a 


<j>  (A2 ) 
6  (B2 ) 


Soil  #2 


Figure  3.1.  General  Problem  Diagram  Cantilevered  Wall. 


3 . 2  Background  Theory 
3.2.1  General  Definition 

A  cantilevered  sheetpile  wall  depends  upon  its  embed¬ 
ment  depth  to  develop  resistance  against  the  overturning 
effect  of  a  soil  backfill.  Cantilevered  walls  develop 
their  strength  through  passive  pressure  in  the  lower  soil 
thus  counteracting  the  active  earth  pressure  in  the  back¬ 
fill.  These  walls  do  not  depend  upon  an  anchor  in  the 
backfill  for  support. 

As  depicted  in  Figure  3.2(A),  cantilevered  walls  rotate 


about  a  point  in  the  lower  soil  [1] . 
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F igure  3.2  («)  Assumed  elastic  line  of  the  sheetpiling;  (h)  probable  and  as  obtained  in  finite-element 
solution  qualitative  soil-pressure  distribution;  (r)  simplified  pressure  diagram  Tor  computational  pur¬ 
poses  (granular  soil  and  no  water  as  show  n).  [1] 


Through  model  testing  and  field  experience,  the  earth 
pressures  mobilized  as  the  wall  rotates  are  shown  in 
Figure  3.2  (b)  [1],  For  ease  of  calculation,  these  pressures 

have  been  simplified,  Figure  3.2  (c) . 

The  classic  solution  of  the  cantilevered  wall  involves 
assuming  a  trial  embedment  depth  and  varying  the  passive 
pressures  in  the  lower  soil  until  the  summation  of  horizontal 
forces  approximate  zero.  Moments  are  summed  about  the 
point  of  zero  shear  (Point  o,  Figure  3.2  (c) )  in  the  lower 
soil.  If  the  net  moment  indicates  the  wall  will  overturn. 


a  deeper  depth  is  assumed  until  a  safe  condition  is  calculated 


Two  methods  for  applying  a  safety  factor  have  been  used. 
The  passive  pressures  in  the  lower  soil  can  be  reduced 
by  a  factor,  or  the  calculated  depth  increased  by  20%  to 
40%  [2]. 

Although  there  are  other  methods,  the  Rankine  theory 
of  earth  pressure  is  used  to  calculate  the  coefficients 
of  earth  pressure. 


k  = 


.  2 

sm  (<*  +  $) 


2 

sin  a  sin(a-6) 


i  -  v 


sin  ( <j>  +  6  )  sin  (<j>- 8) 
sin  (a- 6  )  sin  (ct  +  B) 


k  = 


sin  (  ~  ) 


sin  a  sin(ct  +  6) 


1  _  /sin  ( <j>  +  6  )  sin  ( $  +  8) 
sin  ( ct  +  <5 )  sin (a  +  8 ) 


where , 

a  =  wall  inclination  from  horizontal 
8  =  backfill  inclination  from  horizontal 
6  =  wall-soil  friction  angle 
<f>  =  angle  of  internal  friction 


Similarly,  the  Rankine  equations  for  plastic  soil  behavior 
are  used  to  calculate  the  active  and  passive  states  [3]. 


where 
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a  =  active  earth  pressure 
Op  =  passive  earth  pressure 
y  =  unit  soil  weight 
Z  =  soil  depth 
C  =  cohesion 

Effective  stresses  are  considered  using  bouyant  soil  weights 
in  the  above  equations.  Water  pressure  is  superimposed 
on  the  earth  pressures  when  the  soil  is  saturated. 

In  cohesive  soils,  tension  cracks  will  develop  when 
the  soil  is  allowed  to  expand.  This  is  the  case  in  active 
pressure  2ones.  The  backfill,  or  soil  #1  as  referred  to 
in  this  text,  is  an  active  zone.  The  depth  of  these  tension 
cracks  are  calculated  as  [3] : 


where , 

Zq  =  tension  crack  depth 
Y  =  soil  unit  weight 
c  =  cohesion 

k  =  passive  coefficient 
P 

Any  water  that  may  accumulate  in  the  tension  cracks  is 
considered  in  the  computation  of  the  active  force  in  the 
backfill.  Figure  3.3(A)  illustrates  the  active  pressures 
mobilized  behind  the  wall  for  cohesive  soils.  The  cohesive 
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is  distributed  along  the  entire  height  of  the  wall  (A2) . 

Figure  3.3  (c)  illustrates  the  final  pressure  distribution 
of  the  equivalent  active  force  concept.  In  addition  to 
disregarding  the  negative  cohesive  force  and  its  contribution 
to  moment,  the  equivalent  active  force  method  increases 
the  lever  arm  distance  thus  increasing  the  overturning  moment 
and  ultimately  the  required  penetration. 

3.2.3  Pressure  Calculation,  Soil  #2 

To  complete  the  soil  pressure  in  the  lower  soil  (Soil 
#2),  active  and  passive  pressures  are  calculated  for  each 
side  of  the  embedded  wall.  Figure  3.4  illustrates  the 
various  components  of  the  earth  pressures.  As  previously 
mentioned,  a  factor  of  safety  can  be  applied  by  reducing 
the  passive  pressures.  The  active  pressure  of  one  side 
is  subtracted  from  the  passive  pressure  of  the  other  side. 

The  resulting  combined  pressure  diagrams  are  shown  in  Figure 
3.5.  In  general,  a  granular,  noncohesive  soil  will  have 
a  pressure  diagram  similar  to  Figure  3.5  (b) ;  conversely, 
a  cohesive  soil  will  have  a  combined  pressure  diagram  similar 
to  Figure  3.5  (a).  P2  represents  the  pressure  at  the  soil 
interface  and  will  be  discussed  in  more  detail  in  Section  3.4. 

The  pressure  against  the  embedded  depth  of  wall  is 
varied  by  altering  Line  Ll  which  is  drawn  from  the  lower  right 
hand  side  pressure  diagram  up  to  an  arbitrary  point  on 
the  left  hand  side  pressure  diagram.  The  actual  pressure 


-'V  V 


Active 


Passive 
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distribution  is  approximated  as  in  Figure  3.2.  Line  LI 
is  varied  until  the  summation  of  horizontal  forces  approaches 
zero.  Although  the  active  forces  in  soil  #1  is  emitted 
from  the  diagram,  the  forces  are  considered  in  the  summation 
of  forces  and  moments.  Moments  are  summed  about  the  point 
where  the  line  crosses  the  embedded  wall.  The  net  moment 
indicates  whether  the  assumed  depth  needs  to  be  increased 
or  decreased.  The  entire  sequence  is  repeated  for  a  new 
depth  until  static  equilibrium  is  achieved. 

3 . 3  Programming  Rationale 
3.3.1  Program  Flow 

The  program  flow  is  executed  similar  to  a  manual  calcu¬ 
lation.  As  depicted  in  Figure  3.6,  the  problem  necessitated 
four  main  branches.  Each  branch  is  contingent  upon  the 
location  of  the  ground  water  table.  If  the  water  table 
is  in  soil  #1,  the  program  will  iterate  within  one  branch 
exclusively.  This  is  typically  true  if  the  water  table 
is  specified  to  be  deeper  than  the  originally  assumed  wall 
penetration  and  the  calculated  wall  penetration.  The  program 
must  iterate  between  branches  when  the  water  table  is 
within  the  embedded  wall  depth. 

As  previously  defined,  the  pressures  in  soil  #2  are 
varied  until  static  equilibrium  is  achieved.  The  pressures 
in  soil  #1  do  not  change  during  iterations;  therefore, 
the  program  minimizes  the  number  of  times  the  pressures 
are  calculated.  To  vary  soil  #2  pressures,  the  depth  to 


the  maximum  pressure  on  the  left  hand  side  of  the  combined 
pressure  diagram  is  increased  or  decreased  depending  upon 
the  previously  calculated  net  horizontal  force.  If  the 
net  sum  of  horizontal  forces  is  positive  (to  the  left) , 
the  depth  of  maximum  left  side  pressure  is  increased  thus 
increasing  the  negative  pressure  (left  side  pressure  to 
the  right).  As  described  in  Section  3.4,  the  depth  to 
the  maximum  pressure  is  defined  as  Z4. 

3.3.2  Iteration  by  Slope-Intercept 

To  reduce  the  number  of  iterations  required  to  balance 
horizontal  forces,  a  slope  intercept  method  was  employed 
which  would  increase  or  decrease  Z4  toward  a  projected 
new  Z4  which  corresponds  to  a  net  summation  of  horizontal 
forces  of  zero.  Figure  3.7  illustrates  this  slope-intercept 
concept . 


Figure  3.7.  Slope-Intercept  Search  Method. 


After  using  the  above  routine,  it  was  found  that  only 
three  iterations  were  required  to  balance  the  horizontal 
forces  of  the  system.  This  indicates  a  linear  relationship 
between  the  summation  of  horizontal  forces  and  Z4.  The 
same  method  is  used  for  establishing  the  best  new  embedment 
depth  (Zl)  compared  against  the  summation  of  moments.  The 
program  iterates  up  to  15  times  prior  to  balancing  the  moments 
thus,  Zl  is  not  linearly  related  to  the  summation  of  moments. 
3.3.3  Sign  Convention 

To  simplify  calculations,  the  sign  convention  was 
established  such  that  all  forces  to  the  right  are  negative 
and  forces  to  the  left  are  positive.  Counterclockwise  moments 
are  positive;  therefore,  all  lever  arms  are  positive  except 
the  lever  arm  associated  with  the  lower  right  hand  side 
force  in  soil  #2.  This  lever  arm  is  negative  because  the 
summation  of  moments  is  about  point  0  (Figure  3.2c)  and  as 
previously  mentioned,  forces  to  the  left  are  positive. 

Subroutines  are  used  as  much  as  possible  to  reduce 
repetitious  program  lines.  The  calculation  of  pressures 
in  soil  #2,  the  summation  of  forces  (5000),  and  the  summation 
of  moments  (6000)  are  the  main  subroutines. 
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3 . 4  Program  Use  and  Limitations 
3.4.1  Input 

CANTWALL  1  is  a  user-oriented  program.  All  input 
and  options  are  prompted  by  statements  and  questions  which 
instruct  the  user  that  a  particular  input  is  necessary. 

As  with  all  programs,  the  user  should  be  familiar  with 
the  input  variables  prior  to  beginning  the  run.  This  will 
prevent  inputting  incorrect  or  mistaken  variables.  A  sample 
of  the  prompting  questions  are  presented  in  Section  3.7. 

The  user  must  input  the  following  variables  (refer  to 
Figure  3.1): 

a)  number  of  runs 

b)  wet  weight,  soil  #1 

c)  saturated  weight,  soil  #1 

d)  cohesion,  soil  #1 

e)  Phi  angle,  soil  #1 

f)  wall-soil  friction  angle,  soil  #1  (see  Table  3.1) 

g)  wet  weight,  soil  #2 

h)  saturated  weight,  soil  #2 

i)  cohesion,  soil  #2 

j)  Phi  angle,  soil  #2 

k)  wall-soil  friction  angle,  soil  #2  (see  Table  3.1) 

l)  wall  height 

m)  ground  water  table  depth 

n)  tolerance 

o)  safety  factor 


Interface  Materials 


Mass  concrete  on  the  following  foundation  materials: 

Clean  sound  rock . . 

Clean  gravel,  gravel-sand  mixtures,  coarse  sand... 
Clean  fine  to  medium  sand,  silty  medium  to  coarse 

sand,  silty  or  clayey  gravel........... . . 

Clean  fine  sand,  silty  or  clayey  fine  to  medium 

sand . . . . 

Fine  sandy  silt,  nonplastic  silt . . . 

Very  stiff  and  hard  residual  or  preconsolidated 

clay . . . . . . 

Medium  stiff  and  stiff  clay  and  silty  clay . 

(Masonry  on  foundation  materials  has  same  friction 
factors. ) 

Steel  sheet  piles  against  the  following  soils: 

Clean  gravel,  gravel-sand  mixtures,  well-graded 

rock  fill  with  spalls . 

Clean  sand,  silty  sand-gravel  mixture,  single  size 

hard  rock  fill . . . 

Silty  sand,  gravel  or  sand  mixed  with  silt  or  clay 

Fine  sandy  silt,  nonplastic  silt . . . 

Formed  concrete  or  concrete  sheet  piling  against  the 
following  soils: 

Clean  gravel,  gravel-sand  mixture,  well-graded 

rock  fill  with  spalls . . . . . 

Clean  sand,  silty  sand-gravel  mixture,  single  size 

hard  rock  fill . 

Silty  sand,  gravel  or  sand  mixed  with  silt  or  clay 

Fine  sandy  silt,  nonplastic  silt . 

Various  structural  materials: 

Masonry  on  masonry,  igneous  and  metamorphic  rocks: 

Dressed  soft  rock  on  dressed  soft  rock. . 

Dressed  hard  rock  on  dressed  soft  rock . . 

Dressed  hard  rock  on  dressed  hard  rock . 

Masonry  on  wood  (cross  grain) . . . 

Steel  on  steel  at  sheet  pile  Interlocks . . 


Friction 
factor, 
tan  S 

Friction 
angle, S 
degrees 

0.70 

35 

0.55  to  0.60 

29  to  31 

0.45  to  0.55 

24  to  29 

0.35  to  0.45 

19  to  24 

0.30  to  0.35 

17  to  19 

0.40  to  0.50 

22  to  26 

0.30  to  0.35 

17  to  19 

0.40 

22 

0.30 

17 

0.25 

14 

0.20 

11 

0.40  to  0.50 

22  to  26 

0.30  to  0.40 

17  to  22 

0.30 

17 

0.25 

14 

0.70 

35 

0.65 

33 

0.55 

29 

0.50 

26 

0.30 

17 

p)  option,  tension  cracks  in  soil  #1 

q)  option,  equivalent  active  force  in  soil  #1 

r)  option,  input  coefficient  of  earth  pressure  or 
calculate 

s)  assumed  penetration  depth 

3.4.2  Units 

All  input  is  in  the  units  of  feet,  pounds,  and  degrees 
The  units  may  be  altered  by  changing  line  744  (G5  =  62.4) 
which  establishes  the  unit  weight  of  water  as  62.4  PCF. 
Angles  must  be  input  in  degrees.  Although  the  units  of 
the  input  variables  may  be  changed,  the  print  statement 
will  label  all  output  in  the  original  units. 

3.4.3  Repeat  Runs 

Each  run  is  totally  independent  of  the  previous  run; 
all  input  will  be  required  again.  This  function  saves 
time  by  eliminating  the  introductory  statements  and  the 
system  commands  the  user  must  execute  to  run  the  program. 

Input  variables  b)-l)  are  self  explanatory  and  are 
depicted  in  Figure  3.1.  The  depth  of  the  ground  water 
table  is  taken  from  the  surface  of  soil  #1. 

3.4.4  Tolerance 

The  tolerance  for  depth  calculation  is  recommended 
as  0.1  to  0.01.  The  tolerance  is  used  during  the  summa¬ 
tion  of  forces  and  moments  routines.  The  tolerance  is 
the  maximum  difference  between  a  calculated  depth  and  the 
new  assumed  depth.  For  force  summations,  Z4  is  compared 


and  for  moment  summations,  Z1  is  used. 


3.4.5  Safety  Factor 


CANTWALL  uses  the  safety  factor  to  decrease  the  passive 
pressures  in  soil  #2  prior  to  creating  a  combined  pressure 
diagram.  This  method  was  previously  described  in  Section 
3.2,  when  specifying  a  soft  cohesive  soil  for  soil  #2, 
it  is  recommended  to  use  a  safety  factor  of  1  as  the  resulting 
combined  pressure  diagrams  will  indicate  wall  instability. 

This  is  discussed  in  more  detail  later  in  this  section. 

3.4.6  Options 

The  user  may  exercise  three  options;  i.e.,  specify 
tension  cracks  in  soil  #1,  specify  an  equivalent  active 
force  in  soil  #1,  and  input  the  coefficients  of  earth  pressure 
versus  using  the  values  calculated  by  the  program. 

The  theoretical  depth  of  tension  cracks  as  described 
in  Section  3.2  are  used  when  this  option  is  exercised. 

With  this  option,  the  force  due  to  water  pressure  in  the 
crack  is  included.  The  force  is  calculated  from  hydrostatic 
pressure  for  the  calculated  tension  crack  depth.  The  problem 
checks  the  crack  depth  against  the  ground  water  table  depth 
and  the  wall  height.  The  crack  depth  can  not  exceed  the 
wall  height  or  the  depth  of  the  water  table. 

The  equivalent  active  force  option  is  calculated  as 
described  in  Section  3.2  (Figure  3.3).  Because  the  tensile 
effect  of  cohesion  is  ignored  and  the  lever  arm  of  the 
positive  force  is  increased,  the  overall  effect  of  exercising 
this  option  is  to  increase  the  required  wall  penetration. 


The  program  allows  the  user  to  input  the  coefficients 
of  earth  pressure.  If  the  user  exercises  this  option, 
the  program  will  not  calculate  the  coefficients  as  defined 
in  Section  3.2.  The  user  will  be  required  to  input  the 
active  coefficient  for  soil  #1  and  the  active  and  passive 
coefficient  for  soil  #2.  All  input  will  be  prompted  by 
questions  if  the  option  is  used. 

3.4.7  Output 

Printed  output  consists  of  a  list  of  input  variables 
and  the  required  penetration  depth.  A  supplementary  data 
list  is  available  giving  the  value  of  the  variables  used 
in  the  programs.  After  the  required  penetration  depth 
is  printed,  the  user  is  prompted  by  a  question  asking  if 
the  list  of  variables  is  desired.  Figures  3.8  to  3.13 
serve  as  a  guide  for  interpreting  the  supplementary  variables 
listing. 

To  locate  the  proper  diagram,  the  user  must  know  the 
ground  water  table  depth  (Z3)  ,  the  wall  height  (Z2) ,  and 
from  the  supplementary  list,  know  the  values  of  Zl,  Z4, 
and  P2.  If  Z4  is  less  than  Z2,  Figures  3.8  or  3.9  apply; 
if  P2  is  positive.  Figure  3.9  applies.  Table  3.2  is  provided 
to  easily  identify  the  proper  figure. 


GWT  in  Soil  #1  (-P2) 


Figure  3.12.  Forces  and  Lever  Arms,  GWT  in  Soil  #2  (-P2) 
( Z 3  >  Z4  +  Z2) . 


Figure  3.13.  Forces  and  Lever  Arms,  GWT  in  Soil  #2  (+P2) 
( Z 3  >  Z4  +  Z2 )  . 


Table  3.2.  Figure  Identification 


GWT  Depth 


Z3 

< 

Z2 

Z2 

< 

Z3 

<  Z4  +  Z2 

Z4 

+ 

Z2 

<  Z3  <  Zl 

Z3 

> 

Zl 

(soil  #2) 

Z3 

> 

Zl 

(soil  #1) 

A  separate  figure  is  not  included  for  Z3  >  Zl  because 
the  correlation  of  data  to  the  figures  is  identical  to 
those  already  shown.  When  Z3  >  Zl  use  the  forces  and  lever 
arms  of  Figures  3.8  and  3.9  for  soil  #2;  similarly,  use 
Figures  3.10  and  3.11  for  soil  #1.  Forces  and  lever  arms 
consistently  correspond  to  each  other;  i.e.,  LA  is  always 
the  lever  arm  for  FI,  LE  is  the  lever  arm  for  F5,  and  LJ 
is  the  lever  arm  for  F0 . 

CANTWALL  will  many  times  calculate  a  lever  arm  that 
is  not  used  but  will  be  printed  in  the  supplementary  list. 
As  long  as  the  corresponding  force  is  zero,  the  user  can 
summize  that  the  particular  force  and  lever  arm  was  not 
a  part  of  the  calculation. 

In  the  first  problem  in  Section  3.7,  the  supplementary 


data  list  indicates  FG  =  0#  and  LF  =  24.23  ft.  Referring 
to  Figure  3.8,  the  user  finds  that  F6  is  the  water  pressure 


due  to  tension  cracks  and  that  since  FG  =  0#,  tension  cracks 
were  not  specified  and  that  although  LF  was  calculated, 
it  did  not  affect  the  summation  of  forces  or  moments. 


3.4.8  Error  Warning 

CANTWALL  uses  warnings  to  prevent  the  user  from  specifying 
problem  parameters  that  can  yield  incorrect  solutions. 

The  warning,  "Assumed  depth  must  be  increased"  followed 
by  the  prompting  statement,  "Input  assumed  depth"  will 
be  printed  on  the  monitor  when  Z4  exceeds  Zl.  Referring 
to  Figure  3.8,  it  is  evident  that  Z4  cannot  be  greater 
than  Zl  as  the  combined  pressure  diagram  will  not  reflect 
the  actual  pressure  distribution  of  Figure  3.2(b).  The 
user  must  input  a  new  depth,  preferably  two  to  three  times 
greater  than  originally  assumed.  Input  an  estimated  depth 
according  to  Table  3.3. 

Table  3.3.  Recommended  Assumed  Wall  Depths. 

N (Sand)  Clay  (Cu,  PSF)  Depth  (SF  =  1) 


0  - 

4 

Soft 

(250- 

500) 

2.  OH 

5  - 

10 

Firm 

(500  - 

1000) 

1.5H 

11  - 

30 

Stiff 

(1000  - 

1500) 

1. 25H 

31  - 

50 

Very  Stiff 

(1500  - 

2000) 

l.H 

50  + 


Hard 


(>2000) 


.25  H 


Although  the  number  of  moment  iterations  may  be  reduced 
in  half  by  specifying  a  depth  close  to  the  required  depth, 
it  will  many  times  reduce  actual  user  time  by  inputting 
a  seemingly  large  assumed  depth.  Doing  so,  the  warning 
will  not  appear  and  the  program  will  iterate  the  required 
depth  in  a  relatively  short  time. 

When  the  condition: 

4c  -  q  <  0 

where , 

q  =  surcharge  on  soil  #2 

c  =  cohesion  of  soil  #2 

is  satisfied  the  passive  pressure  below  the  dredge  line 
is  always  less  than  the  active  pressure  (Figure  3.14). 
Therefore,  an  equilibrium  condition  cannot  be  achieved 
no  matter  how  deep  the  sheet  piling  is  driven  below  the 
dredge  line.  This  condition  may  occur  in  soft  and  very 
soft  clays  typically  having  undrained  shear  strengths  less 
than  about  500  psf. 

The  program  identifies  this  condition  when  PI  is  cal¬ 
culated  as  a  positive  number  (see  Figure  3.8).  Upon  identi¬ 
fication  of  this  condition,  the  statement  "‘Warning*  Computed 
Wall  Instability"  is  printed  on  the  monitor.  The  user 
is  given  the  option  to  end  the  program  or  to  re-enter  the 
input  variables.  The  first  correction  would  be  to  reduce 
the  safety  factor  to  one.  If  the  program  successfully  runs, 


Figure  3.14.  Quick,  Undrained  Loading  in  Cohesive  Soil. 

increase  the  depth  by  20%  to  40%  [2] .  The  second  correc¬ 
tion  should  be  to  decrease  the  wall  height  thus  decreasing 
the  value  of  "q."  If  all  else  fails,  the  user  must  change 
the  soil  conditions  in  soil  #2  or  abort  the  design  completely 
The  final  warning  advises  the  user  that  soil  #1  is 
in  tension  and  does  not  mobilize  an  active  force  on  the 
wall.  This  will  occur  when  soil  #1  is  specified  to  have 
large  cohesive  characteristics.  As  described  in  Section 

3.2,  this  phenomena  occurs  when  the  cohesive  effects  of 
« 

the  soil  are  great  enough  to  allow  the  soil  to  stand 


unsupported  to  a  height  equal  to  or  greater  than  the  speci¬ 
fied  wall  height. 

The  only  way  to  avoid  this  condition  is  to  increase 
the  wall  height  thus  increasing  the  active  force  on  the 
wall  due  to  soil  weight,  or  decrease  the  cohesive  charac¬ 
teristic  of  the  soil. 

3.4.9  Limitations 

The  limitations  of  CANTWALL  relate  primarily  to  the 
physical  description  of  the  conditions.  The  most  severe 
limitation  is  not  being  capable  to  specify  a  surcharge 
load  on  the  backfill  (soil  #1) .  Many  practical  applications 
would  warrant  a  surcharge.  If  a  surcharge  exists,  the 
user  can  decrease  the  cohesion  in  soil  #1  and  soil  #2 
by  an  appropriate  amount  such  that  the  net  effect  on  the 
wall  is  the  same.  Since  the  cohesion  creates  a  rectangular 
pressure  diagram  as  does  a  surcharge,  this  could  be  done. 

A  negative  cohesion  may  be  input  to  simulate  a  surcharge 
in  a  cohesionless  backfill.  The  user  must  convert  the 
surcharge  to  a  proper  value  of  cohesion  by  the  following 
analogy: 


aC  =  -2c'  /ka 
a 


cohesion 


=  q  ka 


surcharge 


where , 

c 

a  =  active  earth  pressure  due  to  cohesion 

cl 

=  active  earth  pressure  due  to  surcharge 
c'  =  equivalent  cohesion 
q  =  surcharge 

k  =  active  earth  pressure  coefficient 

Ci 

In  a  cohesive  soil,  the  equivalent  cohesion  is  added 
to  the  cohesion  of  soil  #1  and  soil  #2.  The  user  must 
maintain  a  consistent  sign  convention.  See  Figure  3.15. 


M  H 

M 

c  q 

,  c  q. 

°a  •  a 

a  a 

(a)  Active  pressures  due  to  (b) 

cohesion  and  surcharge 

Equivalent  cohesion 

Figure  3.15.  Surcharge  Represented 
in  a  Cohesive  Soil. 

by  Equivalent  Cohesion 

Another  limitation  is  the  predetermination  of  soil 
layers.  Although  the  user  may  easily  specify  a  homogeneous 
condition,  the  user  is  limited  to  two  soil  types  separated 
at  a  predetermined  depth.  Engineering  judgment  is  the 
only  guide  in  this  case. 

The  calculated  method  used  in  this  program  assumes 
a  rigid  wall  and  therefore  no  effects  of  moment  redistribution 
attributed  to  wall  flexibility  are  considered. 

The  final  limitation  is  contingent  upon  the  fact  that 
any  cantilevered  structure  is  limited  in  length  or  depth 
by  the  large  bending  moments  developed  under  load.  In 
general,  heights  exceeding  15  to  20  ft  are  infeasible. 
Construction  materials  with  section  moduli  high  enough 
to  withstand  the  bending  moments  become  uneconomical  compared 
to  the  cost  of  material  used  in  alternate  construction 
techniques  such  as  anchored  walls  [1] .  This  is  an  engineering 
decision.  The  supplementary  data  list  contains  the  necessary 
data  to  calculate  the  bending  moments  in  the  wall.  Figure 
3.16  provides  an  insight  into  the  section  moduli  required 
to  resist  loads  imposed  under  three  cases. 


Case  A  Case  B  Case  C 

*=■30°  or  35°  ♦  *  30°  or  35°  c  *  500,  1000  or  1500  lb/ft5 

Dry  soil  throughout  Undisturbed  clay  throughout 


q  is  surcharge  in  lb/ft 


Figure  3.16.  Cantilevered  Walls  [5] 
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3 . 5  Program  List 


5  SPEED*  150 

10  PRINT  11  *##*******■“ 

PRINT  "  *CANTNALL*" 

PR  I  NT  "  '* 

PRINT  :  PRINT  :  PRINT 
PRINT  "  DANA  K.  EDDY" 

PRINT  "  60.  INSTITUTE  OF  TECHNOLOGY" 

PRINT  "  SCHOOL  OF  CIUIL  ENGINEERING" 

PRINT  "  DEPORTMENT  OF  GEOTECHNICOL  ENGINEERING" 

PRINT  :  PRINT  :  PRINT 

PRINT  "  SYSTEM  HARDWARE:  APPLE  II  PLUS  <6410 
PRINT  "  SYSTEM  SGFTNORE:  DOS  3.3,  APPLESOFT  BASIC  LANGUAGE" 


11 

12 

13 

2w 

22 

24 

26 

30 

35 

36 


’•I 


37  PRINT  "  PROGRAM  DATE:  MAY,  1933" 


40 

4-5 


PRINT  :  PRINT  s  PRINT  :  PRINT 

PRINT  "  CANTWALL  ESTIMATES  THE  EMBEDMENT  DEPTH  OF  A  CANTILEUERED  NALL. 


50 

PRINT  : 

PRINT  :  1 

55 

SPEED* 

255 

470 

PRINT 

"HOW  MANY 

430 

INPUT 

Q 

435 

DIM  A< 33 DIM 

436 

DIM  C< 

88  >1  DIM 

490 

FOR  R 

=  1  TO  0 

431 

PRINT 

:  PRINT 

495 

AD  =  0 

500 

PRINT 

"WHAT  IS 

510 

INPUT 

61 

511 

PRINT 

:  PRINT 

520 

PRINT 

"WHAT  IS 

530 

INPUT 

62 

531 

PRINT 

:  PRINT 

534 

PRINT 

"WHAT  IS 

533 

INPUT 

Cl 

533 

PRINT 

:  PRINT 

542 

PRINT 

"WHAT  IS 

546 

INPUT 

A1 

547 

PRINT 

:  PRINT 

550 

PRINT 

EES  >" 

"WHAT  IS 

554 

INPUT 

ei 

555 

PRINT 

:  PRINT 

553 

PRINT 

"WHAT  IS 

562 

INPUT 

G3 

PRINT  :  PRINT  :  PRINT  :  PRINT 


< uEbREEs 


DEGR 


.%  ,\ 


’ 


.%  ■  -i  A. 


64 


563 

PRINT 

:  PRINT 

566 

PRINT 

"WHAT  IS  THE 

570 

INPUT 

64 

571 

PRINT 

:  PRINT 

574 

PRINT 

"WHAT  IS  THE 

578 

INPUT 

C2 

579 

PRINT 

:  PRINT 

532 

PRINT 

"WHAT  IS  THE 

586 

INPUT 

A2 

587 

PRINT 

:  PRINT 

590 

PRINT 

II 

"WHAT  IS  THE 

595 

INPUT 

B2 

596 

PRINT 

:  PRINT 

600 

PRINT 

"WHAT  IS  THE 

605 

INPUT 

22 

606 

PRINT 

:  PRINT 

610 

PRINT 

"WHAT  IS  THE 

615 

INPUT 

Z3 

616 

PRINT 

:  PRINT 

£;20 

PRINT  "WHAT  IS  THE 
•  -  .01)?" 

625 

INPUT 

TL 

626 

PRINT 

:  PRINT 

630 

PRINT 

"WHAT  IS  THE 

635 

INPUT 

SF 

636 

PRINT 

:  PRINT 

640 

PRINT 
D  IN  : 

"DO  VOU  WANT 
SOIL  #1?  (YES 

645 

INPUT 

AT-:  X  =  ASC 

646 

IF  X 

(  84  GOTO  671 

647 

PRINT 

:  PRINT 

650  27  =  2  *  Cl  G1 
655  IF  27  <  23  GOTO  665 
660  27  =  22 

665  IF  27  <’  22  GOTO  671 
6 10  2."  —  22 

671  PRINT  :  PRINT  :  PRINT  "DO  VOU  WANT  TO  USE  AN  EQUIUALENT  AOTIUE  FORCE 
IN  SOIL  #1?  < VES  OR  NO  >" 

672  INPUT  CT:W  =  ASC  <CT) 

673  PRINT  :  PRINT 

675  PRINT  "DO  VOU  WANT  TO  INPUT  THE  COEFFICIENTS  OF  EARTH  PRESSURE  (.YES'* 
OR  HAUE  THEM  CALCULATED  FOR  VOU  <NO>?" 

680  INPUT  BT :  X  —  ASC  (.  BT  >:  IF  X  C  84  GOTO  725 

681  IF  X  C  84  GOTO  725 

682  PRINT  :  PRINT  :  PRINT 

685  PRINT  "WHAT  IS  K  ACTIUE  FOR  SOIL  #1?" 

690  INPUT  K1 

691  PRINT 


635  PRIM.  "WHAT  IS  K  PASSIUE  FOR  SOIL  #2?" 

760  INPUT  K2 
701  *  PRINT 

705  PRINT  "WHAT  IS  K  ACTIUE  FOR  SOIL  #2?" 

71Q  INPUT  K3 
715  GOTO  744 
725  RA  =  .01745 

730  K1  =  <  SIN  <1.571  +  A1  *  Rfl))  a  2  /  SIN  <1.571  -  B1  *  RA>  <1  +  <  SII 
<  <01  +  Bl)  *  RA  >  *  SIN  <A1  *  RA)  /  SIN  <1.571  -  B1  *  RA ) )  .5)  -  2 

735  K2  =  <  SIN  <1.571  -  A2  *  RA>>  ~  2  /  SIN  <1.571  +  B2  *  RA)  ✓  <1  -  <  Sli 

<<A2  +  B2)  *  RA >  *  SIN  <A2  *  RA)  ✓  SIN  <1.571  -  B2  *  RA>>  .5)  2 

740  K3  =  <  SIN  <1.571  +  A2  *  RA))  -  2  ✓  SIN  <1.571  -  B2  *  RA)  ✓  <1  +  <  Sli 

< < A2  +  B2 )  *  RA)  *  SIN  <A2  *  RA)  ✓  SIN  <1.571  -  B2  *  RA))  '•  .5)  •  2 

744  G5  =  62.4 

745  PRINT  :  PRINT 

746  IF  Z3  >  Z2  GOTO  747:Z3  =  2  *  Cl  *  K1  ^  .5  *•  Z2  /  <<G1  *  Z3  +  <  G2  -  G5 
)  *  <Z2  -  Z3>)  *  Kl):  IF  ZS  <  Z2  GOTO  750:  GOTO  748 

747  ZS  =  2  *  Cl  /  <61  *  Kl  A  .5):  IF  ZS  <  Z2  GOTO  760 

74S  PRINT  "SOIL  #1  IS  IN  TENSION,  REEUALUATE  THE  COHESION  OF  SOIL  #1  OR  T 
HE  HEIGHT  OF  THE  WALL.  TYPE  <1>  TO  RESTART  THE  PROGRAM. " 

74S  INPUT  AB:  GOTO  491 

760  PRINT  "INPUT  ASSUMED  DEPTH  OF  WALL  PENETRATION  < REFER  TO  USERS  HANUA 
L>." 

762  INPUT  Z1 

763  PRINT  :  PRINT 

764  N  =  1 

770  Z4  =  .72  *  Zl:  GOTO  774 

771  PRINT  "ASSUMED  DEPTH  MUST  BE  INCREASED":  PRINT  :  PRINT 

772  AZ  =  0:  GOTO  760 

774  FI  =  0: F2  =  0:F3  =  0:F4  =  0:F5  =  0:F6  =  0:F7  =  0:F8  =  0:F9  =  9:F0  =  0: 

LA  =  0:LB  =  0:LC  =  0:LD  =  0:LE  =  0:LF  =  0:LG  =  0:LH  =  0:LI  =  0:LJ  =  0 

775  M  =  1 

730  IF  Z3  >  Z2  GOTO  930 

781  IF  W  <  84  GOTO  735 

782  FI  *  <<G1  *  Z3  +  <G2  -  85)  *  <Z2  -  Z3 ) )  *  Kl  -  2  Cl  *  Kl  .5  .•  *  22 

-  2  *  Cl  *  Kl  .5  *  22  /  <G1  *  Z3  +  <G2  -  65)  *  <Z2  -  23  > )  1 ) i 


783  GOTO  805 


'35 

FI 

S 

G1 

it 

Z3  2  *  Kl  /  2 

30 

F2 

= 

61 

it 

Z3  *  Kl  *  <Z2  - 

Z3  > 

'95 

F3 

= 

<  G 1 

*♦ 

•  Z3  +  <  G2  -  G5 ) 

*  <<l  2  ■ 

•00 

F4 

- 

•w. 

*  Cl  *  Kl  -  .5 

it  Tv 

.05 

F5 

G5 

<Z2  -  Z3)  2  ✓ 

-1 

<- 

>10 

F6 

= 

G5 

it 

Z7  -  2  2 

■:  1 5 

FA 

= 

FI 

+ 

F2  +  F3  +  F4  + 

F5  +  F6 

*  Kl  *  v  Z2  - 


THE  NET  ACTIUE  FORCE  AGAINST  THE  WALL  IN  SOI 


816  IF  FA  >  0  GOTO  845 

817  PRINT  "**  WARNING  ** 

L  #1  IS  NEGATIUE. 11 

818,  PRINT  :  PRINT 

619  PRINT  "TO  RESTART  THE  PROGRAM  TVPE  Cl >;  TO  EXERCISE  SOIL  #1  ACTIUE  OP 
TIONS  TVPE  <  2 >. " 

820  INPUT  AA:  IF  AA  =  1  GOTO  491:  GOTO  636 

821  GOTO  636 
845  Z  =  21 

850  GGSUB  6200 
855  Z  =  Z4 
860  GOSUB  6100 

865  PI  =  PL 

866  GOSUB  6800 

867  IF  AO  —  1  GOTO  491 
870  Z  =  0 

875  GOSUB  6108 
880  P2  =  PL 

885  GOSUB  4000 

886  IF  AZ  *  1  GOTO  771 

887  IF  M  =  1  GOTO  894 

690  IF  TL  >  =  ABS  <100  *  < 24  -  C<M>>  /  24  >  GOTO  896 

894  M  =  M  +  1 
893  GOTO  855 

896  IF  W  <  84  GOTO  900 

897  LA  =  25  +  Z2  /  3 

898  GOTO  917 

900  LA  -  Z5  +  Z2  -  2  *  23  /  3 

910  LB  »  25  +  <22  -  23)  /  2 

915  LC  =  25  +  <22  -  23 >  ✓  3 

916  LD  =  25  +  22  /  2 

317  LE  =  25  +  <22  -  23 >  ✓  3 

918  LF  =  25  +  22  -  2  *  27  ✓  3 

320  GOSUB  6000 

921  IF  N  =  1  GOTO  923 

922  IF  TL  >  =  ABS  <  1O0  *  <21  -  A<H>>  Z1  >  GOTO  3O00 

923  N  =  N  +  1 
925  GOTO  770 

930  IF  23  >  22+21  GOTO  2331 
331  IF  W  <  84  GOTO  935 

9j2  FI  =  <22  *  61  +  K1  —  2  *  Cl  *  K1  .  5  >  <22  -  2  +  Cl  /  <G1  +  K  l  .5 

>  •••'  2 

J33  GOTO  345 

935  FI  =  Gi  +  22  2  *  K1  /  2 

940  F2  =  —  2  *  Cl  +  K 1  .5  *  22 

943  F3  —  G5  +  27  2  2 

350  FA  =  FI  +  F2  +  F3 

951  IF  FA  >  0  GOTO  965 

952  PRINT  ”*+  WARNING  **  THE  NET  ACTIUE  FORCE  AGAINST  THE  WALL  TN  9.V 
L  41  IS  NEGATIUE." 


67 


353 

954 


355 

356. 

365 


980 

335 


PRINT  :  PRINT 

PTIONs"tVPER^P»T  THE  PR,;"3RPH  TVPE  <l>;  T0  EXERCISE  SOIL  #1  RCTIUE  0 

INPUT  RB:  IF  SB  =  1  GOTO  491:  GOTO  63R 
GOTO  636 

IF  24  <  =  Z3  -  22  GOTO  2010 

970  2  *  21 
375  GOSUB  6400 
2  =  24 
GOSUB  6300 
990  PI  =  F'L 
331  GOSUB  6800 
992  IF  NO  *  1  GOTO  431 

395  ^5  =  Zt  -  <  PR  *  <  24  -  21  >  ✓  <Pi  -  PR )  , 

1000  2  =  23  -  22 
1 005  GOSUB  6300 
1010  P3  =  PL 
1015  2=0 
1020  GOSUB  6700 
1025  P2  =  PL 

1030  IF  P2  >  0  GOTO  1120 
1035  F4  =  P2  *  <23  -  22) 

1040  F5  =  <P3  -  P2>  *  <23 
1045  F6  =  <P1  -  P3>  *  <24 
105O  F7  =  P3  *  <24  -  23  + 

1055  F8  =  PI  *  <25  -  24)  ✓ 

1060  F9  =  PR  *  <21  -  25  ) 

1065  GOSUB  5000 

1066  IF  fl2  =  1  GOTO  771 

IF  H  =  1  GOTO  1074 
IF  TL  >  =  flBS  < 100 

1074  M  =  M  +  1 
1073  GOTO  380 


22)  /  2 
23  +  22)  / 


22  ) 


1067 
1 070 


<24  -  C<  M  )  > 


24)  GOTO  108S) 


1080 

Lfl 

=  25 

+  22  3 

iUSi 

IF 

M  > 

34  GOTO  1086 

1065 

LB 

=  25 

+  22  ✓  2 

1 036 

LC 

=  25 

+  22  -  2  *  27  ✓  3 

1087 

LO 

-  25 

-  <23  -  22)  /  2 

1 083 

Lc 

=  25 

~  2  *  <23  -  22)  3 

1089 

LF 

=  25 

~  <23  -  22  +  2  *  24) 

i  030 

LG 

=  25 

~  <  ii  j  —  c2  +•  24 )  /  2 

1031 

LH 

=  2  *. 

<  25  -  24  )  s  3 

i  03  "i 

LI 

ZZ  — 

2  *  <21  -  25)  3 

1110 
1111 
1  L 12 
1113  H  =  N 
l  1. 1 


1  1 20 
1 122 


GOSUB  6000 
IF  H  =  1  GOTO  1113 
IF  TL  >  =  mBS  <100  *  <21  -  m<H 

+  1 
hut o  770 

=J_ZZ  ~  22)  *  P2  /  <  F'2  -  P3) 

IF  i-6  >  <23  -  22)  GOTO  1200 


21)  GOT 0  3O0O 


1125 

F4  = 

P2 

* 

Z6  2  2 

1130 

F5  = 

P3 

* 

<Z3  - 

Z2  -  Zi 

S  !> 

2  2 

1135 

F6  = 

<P1 

.  - 

•  P3>  i 

*•  CZ4  - 

23  +  Z 

2 

>  2  2 

114© 

F7  = 

P3 

* 

<Z4  - 

Z3  +  Z; 

2  > 

1145 

F8  = 

PI 

* 

<Z5  - 

Z4>  2 

2 

1150 

F9  = 

PP. 

<Z1  - 

Z5 >  2 

C. 

1155 

60SUB  5000 

1156 

IF  AZ  = 

=  1 

.  GOTO 

771 

1157 

IF  M  = 

1 

GOTO 

1164 

1 160 

IF  TL  > 

=  AB3  < 100 

*  1 

CZ4  - 

•  i 

C<  H  >  >  2 

Z4  ) 

GOTO 

1170 

1164 

M  =  M  + 

1 

1165 

GOTO  980 

1170 

LA  = 

Z5 

+ 

Z2  /  3 

1171 

IF  W  > 

84  GOTO 

1173 

1172 

LB  = 

Z5 

+ 

Z2  2  ; 

1173 

LC  = 

Z5 

+ 

Z2  -  ; 

2  *  Z? 

/  : 

si 

1174 

LD  = 

Z5 

- 

Z6  2  3 

1175 

LE  = 

Z5 

- 

<Z6  + 

2  *  23 

;  - 

2  * 

Z2)  2  3 

1160 

LF  = 

Z5 

- 

2  *  <  Z4  -  Z3 

+ 

Z2> 

/ 

sJ> 

1185 

LG  = 

Z5 

- 

<Z4  - 

Z3  +  22; 

2  2 

1180 

LH  = 

<2 

* 

Z5  +  ; 

Z4)  2  3 

1185 

LI  = 

- 

2 

*  CZl 

-  Z5) 

3 

1 197 

GOTO  2000 

i  200 

F4  = 

P3 

* 

<Z3  - 

Z2> 

1205 

F5  = 

<P2  -  P3> 

*  <Z3  - 

Z 

2)  2 

<*•, 

<1 

1210 

F6  = 

P3 

* 

<  ZG  — 

Z3  +  Z2  :> 

2  2 

1215 

F7  = 

PI 

* 

<Z4  - 

Z6>  2 

2 

1 220 

F3  = 

PI 

* 

<Z5  - 

Z4 ')  2 

2 

1225 

F9  = 

PR 

* 

<Z1  - 

Z5>  2 

2 

1238 

GOSUB  ' 

5000 

1232 

IF 

AZ  : 

=  ; 

l  GOTO 

771 

1235 

IF 

TL 

> 

=  Ae 

8  <  100 

if 

•■.  Z4  ■ 

* 

C<  H  >  >  2 

Z4  > 

GOTO 

1250 

1240 

H  = 

M  + 

1 

1245 

GOTO  9; 

30 

1250 

LA  = 

Z5 

+ 

•CiL  S 

1255 

IF 

H  > 

o, 

4  GOTO 

1265 

1260 

L6  = 

Z5 

+ 

Z2  2 

U. 

1 265 

LC  = 

Z5 

+ 

Z2  - 

2  *  Z7 

/ 

1270 

LD  = 

Z5 

- 

•1 Z3  - 

Z2>  2 

o 

i  275 

LE  = 

‘  Z5 

- 

<Z3  - 

•  Z2  >  2 

"7 

1280 

LF  = 

Z5 

+ 

•12  * 

ii2  ^ 

Z6  ’>  2  3 

1 285 

LG  = 

:  Z5 

- 

2  *  '■ 

Z4  -  Z6  > 

1290 

LH  = 

w 

1 

CZ5  - 

Z4 )  2  3 

1295 

LI  = 

- 

*  ■:  zi 

-  Z5  > 

“7 

1  300 

GOS 

LIB 

6000 

1305 

IF 

M  - 

1 

GOT  0 

1315 

1310 

IF 

TL 

s 

=  AE 

IS  •;  100 

•4b 

•:'.  Z 1 

<N>>  •- 

Z 1  j 

GOTO 

30OG 

1315 

N  = 

H  + 

1 

1 320 

GOTO  7 

70 

2000 

60S 

;|JB 

6000 

69 


2001  IF  N  =  1  GOTO  2O03 

2002  IF  TL  >  =  PBS  <100  *  <Z1  -  R<  N >  ✓  Z1  >  GOTO  3000 

2003  N  =  N  +  1 
2005  GOTO  770 
2010  Z  =  21 

20 1 5  G0SU6  6400 
2020  2  =  24 
2025  GOSUB  6700 

2030  PI  =  PL 

2031  GOSUB  6600 

2032  IF  PO  =  1  GOTO  431 
2035  2=0 

2040  GOSUB  6700 
2045  P2  =  PL 

2050  GOSUB  4000 

2051  IF  PZ  =  1  GOTO  771 

2052  IF  M  =  1  GOTO  2059 

2055  IF  TL  >  =  PBS  <100  *  <Z4  -  C<M>)  /  24 >  GOTO  2065 

2059  M  =  M  +  1 

2060  GOTO  2020 

2065  LP  =  25  +  22  ✓  3 

2066  IF  H  >  84  GOTO  2075 

2070  LB  =  25  +  22  ✓  2 

2075  LC  =  25  +  22  -  2  *  27  ✓  3 

2065  GOSUB  6000 

2086  IF  N  *  1  GOTO  2086 

2087  IF  TL  >  =  PBS  <100  *  <21  -  P<N>>  /  21)  GOTO  3000 

2088  N  =  N  +  1 
2890  GOTO  770 

2331  IF  H  <  84  GOTO  2335 

2332  Fi  =  <22  *  61  *  K1  -  2  *  Cl  *  Kl  ^  .5)  *  <22  -  2  *  Cl  <61  *  K 1 

>  /  2 

2333  GOTO  2345 

2335  FI  »  61  *  22  2  *  Kl  2 

2340  F2  =  -  2  *  Cl  *  Kl  ^  .5  *  22 

2345  F3  =  G5  *  27  ••s  2  /  2 

23o0  FP  —  FI  +  F2  +  F3 

2351  IF  FP  >  0  GOTO  2365 

2352  PRINT  "**  WORN I NG  **  THE  NET  PCTIUE  FORCE  PGPINST  THE  WPLL  IN 
IL  #1  IS  NEGPTIUE. " 

2353  PRINT  :  PRINT 

2354  PRINT  "TO  RESTPRT  THE  PROGRPM  TVPE  <1>;  TO  EXERCISE  SOIL  41  PCTIUE 
PTIONS  TVPR  <  2  > . " 

2355  INPUT  PCs  IF  PC  -  1  bOTO  491:  GOTO  636 

2356  GOTO  636 


2370 

GOSUB  660O 

”?  T  ^ 

Z  =  Z4 

2360 

GOSUB  650O 

■  “TC-cr 

PI  =  PL 

70 


2386  60SUB  6380 

2333  IF  AD  =  1  GOTO  431 

2398  Z  =  0 

2395  60SU6  6500 

2400  P2  =  PL 

2405  60SUB  4000 

2406  IF  AZ  =  1  GOTO  771 

2407  IF  H  =  1  GOTO  2414 

2410  IF  TL  y  =  ABS  <100  *  <  Z4  -  C<H)}  /  Z4.'1  GOTO  2420 

2414  M  =  h  +  1 

2415  GOTO  2375 

2420  LA  =  25  +  22  /  3 

2421  IF  W  >  84  GOTO  2430 
2425  LB  =  Z5  +  Z2  ✓  2 

2430  LC  =  Z5  +  22  -  2  *  27  /  3 

2440  GOSUB  6000 

2441  IF  N  =  1  GOTO  2443 

2442  IF  TL  >  =  ABS  <100  *  <21  -  A<  N  > .')  /  21  >  GOTO  3000 

2443  N  =  N  +  1 
2445  GOTO  770 

3000  L$  =  CHRf  < 4>s  PRINT  L$i"PR#l" 

3808  PRINT  "*****•»****":  PRINT  "******* **•*" 

3010  PRINT  "INPUT  BATA" 

3011  PRINT  "****■***#**":  PRINT  "***•*•****•**" 

3012  PRINT  :  PRINT 

3020  PRINT  "WET  UNIT  WEIGHT ,  SOIL  #1  =  "G1"  PCF. " 

3021  PRINT 

3030  PRINT  "SATURATED  UNIT  WEIGHT,  SOIL  #1  ="62"  PCF." 

3031  PRINT 

3840  PRINT  "COHESION,  SOIL  #1  =  "C1"  PSF. " 

3041  PRINT 

3050  PRINT  "ANGLE  OF  INTERNAL  FRICTION,  SOIL  #1  =  "A1"  DEGREES." 

3051  PRINT 

3060  PRINT  "FRICTION  ANGLE,  SOIL  #1  ="B1'1  DEGREES." 

3061  PRINT 

3070  PRINT  "WET  UNIT  WEIGHT,  SOIL  #2  ="G3"  PCF." 

3071  PRINT 

3080  PRINT  "SATURATED  UNIT  WEIGHT,  SOIL  #2  ="G4"  PCF." 

3081  PRINT 

3090  PRINT  "COHESION,  SOIL  #2  ="C2"  PSF." 

3091  PRINT 

3100  PRINT  "ANGLE  OF  INTERNAL  FRICTION,  SOIL  #2  ="A2"  DEGREES." 

3101  PRINT 

3110  PRINT  "FRICTION  ANGLE,  SOIL  #2  ="B2"  DEGREES." 

3111  PRINT 

3120  PRINT  "WALL  HEIGHT  ="22"  FEET." 

3121  PRINT 

3130  PRINT  "  GROUND  WATER  DEPTH  ="23"  FEET." 

3131  PRINT 

31.33  PRINT  "ACTIVE  K,  SOIL  #1  ="Ki 


71 


3134 

3135 

3136 

3137 
3138* 
314Q 


PRINT 

PRINT  "PASSIVE  K,  SOIL  #2  ="K2 
PRINT 

PRINT  "ACTIVE  K.  SOIL  #2  =  "K3 
PRINT 

PRINT  "TOLERANCE  ="TL"  PERCENT." 


3141 

PRINT 

3150 

PRINT 

" SAFETV  FACTOR  ="SF". 

II 

3151 

3153 

PRINT 

:  PRINT 

PRINT 

"-*■*•**■*•*■■*•■*--*•■*■■*■":  PRINT 

"  *********** " 

3160 

PRINT 

"OUTPUT  DATA" 

3161 

PRINT 

"  *************  "  •  PR  I  NT 

" *********** " 

3162 

PRINT 

:  PRINT 

3165 

zz  =  zi 

*  100 

3166 

ZZ  =  I  NT  <ZZ> 

3167 

ZZ  =  zz 

/  100 

3176 

PRINT 

"REQUIRED  NALL  PENETRATION  ="ZZ"  FEET." 

3171 

PRINT 

:  PRINT 

3180 

PRINT 

L$i"PR#0" 

3180 

PRINT 

:  PRINT  :  PRINT  "DO  VOU  NANT  SUPPLEMENT ARV 

3200 

INPUT 

C*:X  =  ASC  <C*>:  IF 

X  <  84  GOTO  3520 

3210 

3218 

Lf  =  CHRS  <4>:  PRINT  L*j"PR#1" 

PRINT 

3219 

PRINT 

11  i# 

3220 

PRINT 

"SUPPLEMENTARV  DATA" 

3221 

PRINT 

"  ******************" 

3222 

PRINT 

:  PRINT 

3225 

PRINT 

"Zl  ="Z1"  FEET." 

3230 

PRINT 

"Z4  ="Z4 "  FEET." 

3240 

PRINT 

"Z5  ="Z5"  FEET." 

3250 

PRINT 

"Z6  ="Z6"  FEET." 

3260 

PRINT 

"Z7  ="Z7"  FT. " 

3270 

PRINT 

"Pt  ="P1"  PSF . " 

3280 

PRINT 

"P2  ="P2 "  PSF. " 

3290 

PRINT 

"P3  ="P3“  PSF." 

3300 

PRINT 

"P4  ="P4"  PSF. " 

3400 

PRINT 

"PR  ="PR"  PSF. " 

3405 

PRINT 

"FI  ="F1 

","LA  ="LA"FT. " 

3410 

PRINT 

"F2  ="F2"#. 

","LB  ="LB"FT." 

3415 

PRINT 

"F3  ="F3"#, 

","LC  ="LC"FT. " 

3420 

PRINT 

"F4  ="F4"#. 

","LD  ="LD"FT.  " 

3425 

PRINT 

"F5  ="F5"#. 

" # "LE  =”LE"FT.  " 

3430 

PRINT 

"F6  ="F6"#. 

"  .  "LF  ="LF"FT.  " 

3435 

PRINT 

"F7  ="F7"#. 

"."LG  ="LG"FT. " 

3440 

PRINT 

"F8  ="F8"#. 

","LH  ="LH"FT. " 

3445 

PRINT 

"F0  ="F9"#. 

"."LI  ="LI "FT. " 

3450 

PRINT 

"FO  = " F 0 " # . 

" . "L  J  ="L J"FT . " 

3460 

PRINT 

"FT  ="FT"#, 

" , "MT  ="MT"FT-#. " 

3500 

7505 

PRINT 

r»OME 

Li i "PR#0" 

72 


35 10  NEXT  R 

352U  PRINT  "THANK  VOU  FOR  USING  CONTHOLL. " 
3521.  PRINT 

3522  PRINT  "BVE-BVE" 

3523  END 

4000  Z5  =  Z1  -  < PR  *  CZ4  -  Zl>  x  <P1  -  PR)) 


4005 

IF  P2  >  =0  GOTO  4040 

4010 

F7  =  P2  *  Z4 

4015 

FS  =  <P1  -  P2>  *  Z4  x  2 

4020 

F3  =  PI  <25  -  Z4 ;•  x  2 

4025 

F0  *  PR  *  <21  -  Z5>  x  2 

4026 

LG  =  Z5  -  Z4  x  2 

4027 

LH  =  25  -  2  4  24  x  3 

402 S 

LI  =  2  *  < Z5  -  Z4 )  x  3 

4023 

LJ  =  -  2  *  <Z1  -  Z5>  x 

4030 

GOSUB  5000 

4035 

RETURN 

4040 

4.6  =  —  P2  4  Z4  x  <  PI  - 

P2) 

4045 

F7  =  P2  4  26  x  2 

4050 

F3  =  PI  4  <24  -  Z6 :>  x  2 

4055 

F3  =  PI  *  <25  -  Z4>  x  2 

4060 

FO  =  PR  *  <Z1  -  25  >  x  2 

406 1 

LG  =  Z5  -  Z6  x  3 

4062 

LH  =  25  —  <  2  4  24  +  26 ) 

/  -7 

4063 

LI  =  2  4  <25  -  24  >  x  3 

4o64 

LJ  =  -  2  4  <21  -  25)  x 

3 

4065  GO  SUE:  5000 
4070  RETURN 

5000  FT  =  FI  +  F2  +  F3  +  F4  +  F5  +  FS  +  F7  +  F8  +  F9  +  F0 

5605  C< H  >  =  Z4 

50 1 0  D<  M  >  =  FT 

5015  IF  H  >  =2  GOTO  5055 

5020  IF  Q<  M  >  <  0  GOTO  5040 

50i25  Z4  =  Z4  +  2 

5035  RETURN 

5040  Z4  =  Z4  -  2 

5030  RETURN 

5055  Z4  =  C<  H  -  1  >  -  D<  M  -  1  >  *  <  C<  H  >  -  C<  H  -  1  >  >  x  d<  H  > 

5056  IF  Z4  <  =  Z1  THEN  5065 

5060  HZ  =  1 

5065  RETURN 

6000  MT  =  FI  *  LH  +  F2  *  LB  +  F3  -  LC  +  F4  -  LD  +  F5  *  LE 

LG  +  FS  4  LH  +  Fa  *  LI  +  F0  4  LJ 

6001  PRINT  "I  MM  COMPUTING,  PLEHSE  BE  POT I ENT.":  FRINT  : 
6005  H< N  1  =  Z1 

6010  B<  N  >  =  MT 

6015  IF  N  >  =  2  GOTO  6055 

6O20  IF  B<  N  ■  <  0  GOTO  6040 

6025  Z1  =  Z1  +  2 

6035  RETURN 


-  0<  M  -  1  1  ■ 

+  F6  4  LF  f 
PRINT 
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6040  Z1  =  21  -  2 
6050  RETURN 

&U55  Zi  =  «<N  -  1  )  -  B<  N  -  1)  *  <  fl<  N  >  -  fi<N  -  1  >>  B<  H  ">  -  R<  N  -  L  v, 
6076  RETURN 

blU0  PL  =  —  64  —  65 )  +  2  *  K2  +  Z  *  C2  *  K2  ••'••  .5  +  65  * 

G5 >  *  2  ■*  K3  +  <61  •*•23  +  <  62  -  65  )  *  <  22  -  23  )  )  *  K3 

+  Z )  -  2  *  C2  *  K3  •"••  .  5 

6105  RETURN 

6200  PR  =  <<64  —  G5  >  *  Z  *•  K2  +  2  ■*  C2  *  K2  .5  +  <61  *  23  +  <62  —  G5 ) 

<  22  -  23  >  >  *  K2  +  65  *  <  22  -  23  +  2  > >  /  SF  +  2  +  C2  +  K3  •■••■  .5  -  <  64 

65  )  *  2  +  K3  —  65  +  Z 
6205  RETURN 

GbUU  PL  —  —  'i.  4!  *  L-2  +  k.2  ••*••  .5  +  63  +  <  Z3  —  22  )  +  K.2  +  65  +  <2  —  23  t  22 

<64  —  65  )  •*  <2  -  23  +  22 )  •*  K2 !>  •••''  SF  +  63  +  <23  -  22)  +  K3  -  2  *  C2 

K3  ■■■•■  •  5  +  65  *  <2  —  23  +  22 )  +  <  64  —  65 )  +  <2  —  23  +  22 )  +  k.3  •*■  61 

22  *  K.3 
6305  RETURN 

t>4fci0  PR  =  <63  +  <23  —  22)  +  K2  +  2  +  02  +  K2  •'••  .5  +  61  *  22  +  K2  +  <64  — 

5)  *  <  2  -  23  +  22)  *  K2  +  65  +  <  2  -  23  +  22))  •••"  SF  +  2  *  C2  +  k  3 

-  63  +  23  -  22 )  +  K3  -  <  64  -  65 )  +  <2  -  23  +  22 )  *  K3  -  65  +  '•  2  - 

3  +  22  ) 

6405  RETURN 

6“'Uy  F'L  =  —  <  63  ^  Z  +  K2  +  2  +  C2  +  K2  •■■■  .  5 )  SF  +  63  +  2  +  k3  -  2  —  6, 

K3  •••••  .5  +  61  *  22  +  K3  "  -------- 

6505  RETURN 

&600  PR  =  <  63  +  2  +  K2  +  2  +  C2  *  K2  ••'••  .5  +  6 1  +  22  +  K2 )  SF  +  2  -  62 
K3  ■-•  .5  -  63  +  2  *  K3 
66615  RtTURN 

:*7uu  PL  =  —  <63  *  2  +  K2  +  2*6-2  *  K2  •••••  .  5  1  SF  +  63  +  2  ■*•  k3  —  S  -  f_. 

k.3  .5  +  6 1  +  22  +  k.3 

6705  RETURN 

S800  IF  PI  <  0  GOTO  6835 
6805  PRINT  :  PRINT 

6810  PRINT  "**  MORNING  **  COMPUTED  MOLL  INSTOBILITV  <  SEE  USERS  MANUA 

It 

6  8  L  5  PR  I N  T  :  P R I N  T 

:*b26t  PRINT  "TO  RESTART  THE  PR06R0M  TVF'E  <  I.  ) »  TO  END  THE  PR06R0M  TVPE  ■  2 

6825  INPUT  AD:  IF  AD  =  1  GOTO  6835 
6830  END 
6835  •RETURN 


+  N 
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3.6  Variable  List  (CANTWALL  1) 


Input 


Q  =  #  of  runs 


G1  =  Wet  weight,  soil  #1 

G 2  =  Saturated  weight,  soil  #2 

Cl  =  Cohesion,  soil  #1 

Al  =  Phi  angle,  soil  #1 

Bl  =  Delta  angle,  soil  #1 

G3  =  Wet  weight,  soil  #2 

G4  =  Saturated  weight,  soil  #2 

C2  =  Cohesion,  soil  #2 

A2  =  Phi  angle,  soil  #2 

B2  =  Delta  angle,  soil  #2 

Z2  =  Wall  height 

Z3  =  Depth  to  GWT 

TL  =  Tolerance 

SF  =  Safety  factor 

K1  =  Active  earth  pressure  coefficient,  soil  #1 
K2  =  Passive  earth  pressure  coefficient,  soil  #2 
K3  =  Active  earth  pressure  coefficient,  soil  #2 
Flow  Control 
AZ  =  Z4  >  Zl 

X  =  Question  input 

W  =  Equivalent  active  force  in  soil  #1 
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AB  =  Restart  or  end  program 
AA  =  Restart  or  end  program 
AC  =  Restart  or  end  program 
AD  =  Restart  or  end  program 
Counters 
R  =  Run  number 
N  =  Moment  iterations 
M  =  Force  iterations 
Miscellaneous 


Z1 

= 

Wall  penetration  depth 

Z4 

= 

Pressure  depth  in  soil 

#2 

Z5 

= 

Pressure  depth  in  soil 

#2 

Z6 

= 

Pressure  depth  in  soil 

#2 

Z7 

= 

Tension  crack  depth  in 

soil  #1 

Z8 

= 

Tension  depth  in  soil 

#1 

Z 

= 

Depth 

G5 

= 

Unit  weight  of  water,  i 

62.4  pcf 

F1-F0 

= 

Forces 

LA-LJ 

= 

Lever  arms 

FA 

= 

E  forces  in  soil  #1 

FT 

= 

£  forces,  total 

MT 

= 

E  moments 

C  (M) 

= 

Z4 

D  (M) 

= 

FT 

A  (N) 

= 

Z1 

3.7  Problem  Verification 
3.7.1  Problem  #1 

HOW  MANY  PROBLEM  SETS  00  VGU  WANT  TO  RUM? 


WHAT  1$  THE  WET  WEIGHT  OF  SOIL  #1?  <PCF> 
?110 


WHAT  IS  THE  SATURATED  WEIGHT  OF  SOIL  #1?  <PCF) 
7122.4 


WHAT  IS  THE  COHESION  OF  SOIL  #1?  <PSF> 
?20O 


WHAT  IS  THE  ANGLE  OF  INTERNAL  FRICTION  FOR  SOIL  #1?  < DEGREES  > 
?30 


WHAT  IS  THE  FRICTION  ANGLE  BETWEEN  THE  WALL  AND  SOIL  #1?  (.'DEGREE 


WHAT  IS  THE  WET  WEIGHT  OF  SOIL  #2?  <PCF> 
?U0 


WHAT  IS  THE  SATURATED  WEIGHT  OF  SOIL  #2?  <PCF> 
7122.4 


WHAT  IS  THE  COHESION  OF  SOIL  #2?  CPSF> 
7600 


WHAT  IS  THE  ANGLE  OF  INTERNAL  FRICTION  FOR  SOIL  #2?  < DEGREES ) 
730 


WHAT  IS  THE  FRICTION  ANGLE  FOR  THE  WALL  AND  SOIL  #2?  < DEGREES) 
*15 


WHAT  IS  THE  HEIGHT  OF  THE  WALL?  < FEET > 

7  20 


WHAT  IS  THE  DEPTH  OF  THE  GROUND  WATER  TABLE7  (FEET:- 
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WHAT  IS  THE  TOLERANCE  FOR  DEPTH  CALCULATION  <  PERCENT ,  REC'D  .1  -  .01 >? 
?.  05 


HHAT  IS  THE  SAFETY  FACTOR 
?1 


DO  YOU  HANT  TENSION  CRACKS  WITH  WATER  PRESSURE  TO  BE  CONSIDERED  IN  SOIL  #1?’ 
OR  NO  > 

?N 


DO  YOU  WANT  TO  USE  AN  EQUIUALENT  ACTIUE  FORCE  IN  SOIL  #1?  < VES  OR  NO> 
?N 


DO  YOU  WANT  TO  INPUT  THE  COEFFICIENTS  OF  EARTH  PRESSURE  < YES >  OR  HAUE  THEN 
LATED  FOR  YOU  <NO>? 

?N 


INPUT  ASSUMED  DEPTH  OF  NALL  PENETRATION  < REFER  TO  USERS'’  MANUAL). 
?  1 00 


I  AM  COMPUTING,  PLEASE  BE  PATIENT. 


I  AH  COMPUTING,  PLEASE  BE  PATIENT. 


I  AM  COMPUTING,  PLEASE  BE  PATIENT. 


I  AM  COMPUTING,  PLEASE  BE  PATIENT. 


I  AM  COMPUTING,  PLEASE  BE  PATIENT. 


I  AM  COMPUTING,  PLEASE  BE  PATIENT. 


I  AM  COMPUTING,  PLEASE  BE  PATIENT. 


I  AM  COMPUTING,  PLEASE  BE  PATIENT. 


I  AM  COMPUTING,  PLEASE  BE  PATIENT 


I  AM  COMPUTING,  PLEASE  BE  PATIENT. 

I  AM  COMPUTING,  PLEASE  BE  PATIENT. 

I  AM  COMPUTING,  PLEASE  BE  PATIENT. 

I  AM  COMPUTING,  PLEASE  BE  PATIENT. 

I  AM  COMPUTING,  PLEASE  BE  PATIENT. 

^iiitititititititit 
it •£ it it it it it it it it 

INPUT  DATA 
***■**■*■***•»!• 

#*■**■■*#* 

MET  UNIT  WEIGHT.  SOIL  #1  =119  PCF. 

SATURATED  UNIT  WEIGHT,  SOIL  #1  =122.4  PCF. 
COHESION,  SOIL  #1  *200  PSF. 

ANGLE  OF  INTERNAL  FRICTION,  SOIL  #1  =30  DEGREES. 
FRICTION  ANGLE,  SOIL  #1  =15  DEGREES. 

WET  UNIT  WEIGHT,  SOIL  #2  =110  PCF. 

SATURATED  UNIT  WEIGHT,  SOIL  #2  =122.4  PCF. 
COHESION,  SOIL  #2  =600  PSF. 

ANGLE  OF  INTERNAL  FRICTION,  SOIL  #2  =30  DEGREES. 
FRICTION  ANGLE,  SOIL  #2  =15  DEGREES. 

WALL  HEIGHT  =20  FEET. 

GROUND  WATER  DEPTH  =10  FEET. 

ACTIUE  K,  SOIL  #1  =.301403678 


PASSIUE  K#  SOIL  #2  =4.97549393 
ACTIUE  K,  SOIL  #2  =.301403678 
TOLERAMCE  =.05  PERCENT. 

SAFETY  FACTOR  =1. 


■it-******-**** 

■*•■*■** 

OUTPUT  DATA 
*********** 


REQUIRED  HALL  PENETRATION 


DO  VOU  WANT  SUPPLEMENTARY 
?Y 

************** **** 

SUPPLEMENTARY  DATA 

•*• it *  Or it -fr-jr  •£  it & * it ■£ 


Z1  =5.21371744  FEET. 
Z4  =3.99707472  FEET. 
Z5  =4.23195843  FEET. 
Z6  =0  FEET. 

Z7  =0  FT. 

F‘l  =-3320.07579  PSF. 
P2  =-2199.1145  PSF. 
P3  =0  PSF. 

P4  =0  PSF. 

PR  =13880.0615  PSF. 
FI  =1657. 72023#. 

F2  =3315.44046#. 

F3  =904.211035#. 

F4  =-4392.01951#. 

F5  =3120#. 

F6  =0#. 

F7  =-3790.02499#. 

F8  =-2240.283#. 

F 9  =-389.915353#. 

FO  =6314.37162#. 


DATA? 


LA 

=  17 

.565291 

3FT. 

LB 

=9. 

2313584 

3FT. 

LC 

=  r  ■ 

5652917 

6FT. 

LD 

=  14 

.231358 

4FT. 

LE 

=  %  • 

5652917 

6FT . 

LF 

=24 

. 2319584FT . 

LG 

=2. 

2334210 

i7FT. 

LH 

=  1. 

5672418 

5FT. 

LI 

=  .  1 

56589 137FT. 

LJ 

=- . 

6546437 

’47FT 

So  I  L  *1 


fa  -  no  pc f  j  r^=\zz^pc(= 
%  c.  -  200  psF  j  ^  r  30°  ,  S 


\01L*Z  110  PCF  >  '£2< 

c  =  t>00  PSF  .  6  --  30“  ,  A  «  JS" 


4<3<  C  IaJ  (<*  -J)  [  |  f  l&M  (<t  +£)  StAj  (<J>  -  P 

l/s.o  (-  -J)  S 


ShU*  (^0^30°) 


£„o’  ^  (^0~>s)[\  +.  /  £iu  (SOt-is)  S./o  (30-0) 

/ S/aj  (‘iO'/s')  £iaj  (qo-^o) 


s.-oxA  Si/o  a*  4)  n  -  wiy^j 

L  /  SiaJ(~+-£)  Ci/j  (<*-*-(3) 


S./Wf0'3O 


^ /aJ  (‘j'Q-r/S’)  U  -  j  Sitj  (30-ris’)  Si  AO  (.^O  T  o) 
\/  s.im  (7o+/5")  Siaj 


££FER  To  FJ6-U/26 
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-V 


M 


a 


i. 


» 


(4) 


Fi  -  (o-urr)  z /z  ■  ,  30J4  (no  Pc?)  (io‘)  ^/z  -  IL  57,  7  * 

FZ '  ki  y  (4wr)(K-^£Jr)--.3oi 4  (  i  io)  (lojUd'-ic1)  -  35VKS~  ±  “ 

^3  -  k*  <{  (  H-(iUjt)'z/2  z,2oi4  (/Z2.4- £2.<t)  (2o'~/a‘)l‘/z 

-  *!o4.l  * 

F4  r  "2c  "A  (Y.CO  )  { .io\4  (2  o')  r  "4  3*7  2.0  d 

“  J\o  (*-&\AtY/2  --  (z?,4  ( 20  '-/o')  Vz  -  7/24* 


Composite  P/i£S5  w?/=  @£l-£‘.Zl'  ClF-t.o) 


EJMI  (pass we) 

-/V  -  )e^  -  (liz.4 -62.4)  S.zr  (4.174)  -  \SS’4.8  *77- 

-  Z  o  ( =  (?.  ( (o  o o  f> i f-')  \j 4.474  -  Z L>  lo  ■  3  *~  / p>  *- 
-^c  fy;  ~(U-6tJr)(lr)~l /?^  =  (n>‘(iZT.'1-4,Z.4)  *  (t/>  */o)  ;  //o'  : 

-  a 

/?/j0  -  (?i  *  tH'bujf))  y^a  -  fs.z/ ‘i-Zo'-ZJ’)  62.4  =  ■?«?<?./  *7-^ 

i'/’r  -  /  36  36./  T 
^-£Pr  £/t>e.  (active) 

~  ~  Y  (/y  1)  22 . 4  -  h  2 .4)  l  S’ '21')  ( .30  >-4 J  -  .  2  -' 

A  *  ^  c  /YY  -  Yoo)  /,30l4  -  6S8.S  H/Ft^ 

~  y ^  *  *  (0  l  ,  4  ( g",  21  ‘ )  -  3  2  s’,  1  e /pi-  “ 

--  ^39.5 


*2 -fr*  ^ 

136  34 . 1  + 

133  75.6  ^Ft*- 


ComPoS  iTf=  P/ZSSSu&e.  @  2-  -  24  -  4.0 


24  =  4.0 
1 


C  /  .r  £  6- 


FV?S5  i  v£ 


ten  v£ 


Lfz-F-T  5/bE  (pnmue.) 

^>r  -  (*4)  “  {lZl.4'62.4)  (4.0‘){4474)  I"±L_ 

-t?  c  -  '  2.  c.  f  fz^o  -  ~<4  (isOO^  {4-4  74  ~  ■  ^6  7  b  ■  3  /f-t 


6Z.<?  •  4,2)' 


-  'Z41.U  «b 


~  4  U  4 . 7  */?+■ 


L>  4  &  .  4 


1Zi(*n-r  $/&&  (rtcr/UE) 

(24)  k *.  c  ('27.4-6L  ■ 4 )  (4,0)  ,30i4  =  2Pl£l  ^ 

-  -p.c  />*_  1  ^>00)  ]/7ioi4  -  ~  4  £  £  .  £ 

--  £6^1- (?')+(» -6vr) 

~-  C'o'(<*°)  <■  (zo'-/a-)  (no)  -  r/;.<  ff/F^ 

7^.4 a  -  (?4  *-  [tt-bur))  f^o  -  (4.0  +  26- 10')  (olA  --  .  1 22 . 
--  744.S4  *lrt' 
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PI  =  £  ^ 

-  ~  3  3  20.2  Ft  ? 

Composite  ?f.essuize  @  €  -  o 


*a*;  v 

/MW 

c 

Pri*$n/e 

/)<zr> 

LEFT  <,jJ>E  (p*56/u£) 

*  -2  c  CZ^  3  (ioDo)  fiCpiF  -  -2.1,76.3, 

/Z/loMr  F> /D£  (/)  cr/  i/£ ) 

f&  .  -2c  /Z,  =  '2  (c,oo)fjZi4  =  -45-^.g 

=  (\o(\Z?.<1-6Z.<l)  t(zo'-to')  (lio)J.30l4  -  Siz.4  */Ft 

f>E°  ~  6“r)  faxo  r  (lo' -10')  kZ,  <i  -  AC±^/si 

2^  4  77.  £  *>**• 


Px.  "  ±^-1  -  2.  -f?r 

-  -2/48.7  7^- 


/AD-A129  271 

UNCLASSIFIEI 

COMPUTER 
RIR  FORCE 
D  K  EDDV 

APPLICATIONS  TO  GEOTECHNICAL  ENGINEERINGS) 
INST  OF  TECH  HRIGHT-PATTERSON  AFB  OH 

AUG  82  AFIT7CI7NR-82-86T 

F/G  13/2 

2/3  ^ 

NL 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

1 

M 

1 

■ 

^91 

flrzwj  (fi-EFF/Z  TO  E/^U/ze.  ) 


FotZCE 

LB^FK  ,n/zy* 

FI 

ZF  =  'h(io')i-  )o'  +4.  Z3  ’ 

-  n.t. 

FZ 

za  --  %(>o)  +  4.23' 

- 

F3 

LZ  s  7z  C‘0')  1-  4 .23' 

--  7, 

F< 

ZD  -  4  i^o)  1  4.Z3' 

=  /4 

FF 

LB  -  '/ 3  (id')  +  4.23  ' 

=  7. 

FI 

L&  =  7  1.(4')  +  (4.23-4’) 

=  z 

F% 

L  ti  ^7j  (4‘)+ (4.23  -4) 

=  / 

F4 

LX  *z/3  ( 4.Z3-4 )  v 

• 

*4 

LT  %  (f.U  -  4.? 3) 

• 

(*  LJ  IS  /JEZiflTl  V£  for  £W\  ^)  ) 


^  m  ^  <9  £S 

2W\  --  Fi  (la  V  F2.(ls)t  F3^Lc)  *  F4  £lb)  FBs(lb) 

+-  Fi  (  L6.)  ♦■  Fg  £lh)  «•  £li)  ^  F  <4>  ( L l) 

-  /65-7.7  (ll.St,)  i  33fS-.4  fc.zi)  *  7o4 .1.  (?.Fl) 

~  W 2.0(14.23)  +  3>2o(l.SL)  -  g7<K.g  (Z>&) 

-  ZZ4?  .  O  tl  -SC,)  -  33Q.sr(.'S’3)  4-  677l.3(-  -6S3) 

-  4  2.  S(,  Fir  -  dr  ~  O 

EEQuiee?)  dbpth  -  s.  2 1  1 


3.7.2  Problem  #2 


********** 

INPUT  DATA 

+**■*•***■*■•*•* 

MET  UNIT  HEIGHT,  SOIL  #1  =110  PCF. 

SATURATED  UNIT  HEIGHT,  SOIL  #1  =122.4  PCF. 
COHESION,  SOIL  #1  =0  PSF. 

ANGLE  OF  INTERNAL  FRICTION,  SOIL  #1  =32  DEGREES 
FRICTION  ANGLE,  SOIL  #1  =14  DEGREES. 

NET  UNIT  HEIGHT,  SOIL  #2  =100  PCF. 

SATURATED  UNIT  HEIGHT,  SOIL  #2  =112.4  PCF. 
COHESION,  SOIL  #2  =1000  PSF. 

ANGLE  OF  INTERNAL  FRICTION,  SOIL  #2  =23  DEGREES 
FRICTION  ANGLE,  SOIL  #2  =15  DEGREES. 

HALL  HEIGHT  =20  FEET. 

GROUND  HATER  DEPTH  =23  FEET. 

ACTIUE  K,  SOIL  #1  =.28006464 
PASS I UE  K,  SOIL  #2  =4.43311107 
ACTIUE  K,  SOIL  #2  =.325042267 
TOLERANCE  =.05  PERCENT. 

SAFE TV  FACTOR  =1. 


•*r  t*r  it  •fr  ir  ir  ■&  ib  i*  fr  "Sr  ■$*  * 

SUPPLEMENT PF:V  OhTn 


Z1  =6.15013255  FEET. 
Z4  =3.67646922  FEET. 
Z5  =4.32224336  FEET. 
Zb  =0  FEET. 

Z7  =0  FT. 

P 1  =—6047 . 833 1 3  PSF . 
P2  =-4659.82326  PSF. 
F'3  =-5907.2489  PSF. 
P4  =0  PSF. 

PR  =17140.5978  PSF. 


FI  =6161.42208#. 

Ln 

=10. 3839 15FT. 

F2  =0#. 

LB 

=  14. 3222484FT . 

F3  =0#. 

LC 

=84 . 3222484FT . 

F4  =—13979.4848#. 

LB 

=2. 82224836FT. 

F5  =-1371.13036#. 

LE 

=2 . 32224836FT . 

F6  =—47.5694097#. 

LF 

=.87 1268881 FT. 

F7  =-3336 . 07207# . 

LG 

=.3840 13751 FT. 

F3  =-1952.80034#. 
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3.7.3  Problem  #3 


INPUT  DATA 


MET  UNIT  HEIGHT  ,  SOIL  #1  =110  PCF. 

SATURATED  UNIT  HEIGHT  ,  SOIL  #1  =122.4  PCF. 
COHESION,  SOIL  #1  =0  PSF. 

ANGLE  OF  INTERNAL  FRICTION,  SOIL  #1  =32  DEGREES 
FRICTION  ANGLE,  SOIL  #1  =15  DEGREES. 

NET  UNIT  HEIGHT,  SOIL  #2  =100  PCF. 

SATURATED  UNIT  HEIGHT,  SOIL  #2  =112.4  PCF. 
COHESION,  SOIL  #2  =1000  PSF. 

ANGLE  OF  INTERNAL  FRICTION,  SOIL  #2  =28  DEGREE: 
FRICTION  ANGLE,  SOIL  #2  =15  DEGREES. 

WALL  HEIGHT  =20  FEET. 

GROUND  HATER  DEPTH  =60  FEET. 

ACT  I1-. IE  K,  SOIL  #1  =.279043023 
PASS I UE  K,  SOIL  #2  =4.43311107 
ACTIUE  K,  SOIL  #2  =.325042267 
TOLERANCE  =.05  PERCENT. 

SAFETV  FACTOR  =1. 


-*•*•*■■*•*■•*••*•***■>■■ 
OUTPUT  DATA 


REQUIRED  WALL  PENETRATION  =6.1  FEET. 


■+■  -*r  'fr + •*-  ★  ■*■  ★  -it-  4r  *  -fr  -r<r  ■*•  :*• 

SUPPLEMENTARY  DATA 


21  =6.16066258  FEET. 
Z4  =3.68574902  FEET. 


Z5  =4.31743173  FEET. 

Z6  =6.07054395  FEET. 

Z7  =0  FT. 

Pi  =-6196.54613  PSF. 

P2  =-4659.82826  PSF. 

P3  =-21292.1035  PSF. 

F'4  =0  PSF . 

PR  =17774.4774  PSF. 

FI  =6139.0565#. 

LA 

=10. 33403S4FT. 

F2  =0#. 

LB 

=14. 3174317FT. 

F3  =0#. 

LC 

=24 .317 43 1 7FT . 

F4  =0#. 

LD 

=0FT. 

F5  =0#. 

LE 

=0FT. 

F6  =0#. 

LF 

=0FT . 

F7  =-17221.5557#. 

LG 

=2. 46955722FT . 

F3  =-2839.66173#. 

LH 

=  1 . 85359905FT . 

F3  =-1926.1423#. 

LI 

=. 414455141FT. 

FO  =15348.3038#. 

LJ 

=-l. 18S84349FT 

FT  =3. 31469727E-06#. 

NT 

=-.66353391 1FT 
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CHAPTER  IV 


SLOPE  STABILITY 


4 . 1  Problem  Definition 

The  program,  BISH0P1,  calculates  the  stability  of  earth 
slopes.  The  Bishop  method  of  slices  is  used  (Bishop,  1955) . 
The  slope  is  assumed  to  fail  in  a  well  defined  circular 
arc.  The  earth  mass  is  assumed  to  remain  in  a  solid  state 
at  the  time  of  failure.  See  Figure  4.1  for  the  generalized 
geometry  of  the  problem.  A  reservoir  or  pool  may  be  specified 


Failure  Surface 


Reservoir 


Soil  #1 


^  Soil  #2 


£oil  #3 


Soil  #4 


Figure  4.1.  General  Problem  Definition 


103 


c_ 

u 


r.-. 

ft 


r-; 

t  - 
v.- 


T* 


•  r 

V 

li 


I 

£ 


A' 


-T. 


;;; 


tj 


4 . 2  Background  Theory 

4.2.1  General 

A  landslide  is  the  failure  of  a  soil  mass  beneath  a 
slope  [1] .  A  slide  involves  the  downward  and  outward  move¬ 
ment  of  the  soil  mass.  Slides  occur  at  various  rates.  Some 
slides  are  sudden  and  provide  no  warning;  conversely,  other 
slides  fail  slowly  producing  cracks  and  other  visible  signs 
of  impending  failure.  The  safety  of  the  earth  mass  against 
failure  is  termed  its  stability  [2] .  The  stability  of  earth 
masses  must  be  considered  whenever  the  possible  failure 
of  a  slope  may  damage  a  structure  or  cause  harm  to  individuals. 

Slope  failure  occurs  when  the  shear  strength  of  the 
soil  is  exceeded  by  the  shear  stresses  distributed  over 
a  finite  continuous  surface.  Among  the  major  factors  which 
influence  slope  stability  are:  failure  plane  geometry,  non¬ 
homogeneity  of  soil  layers,  tension  cracks,  dynamic  loading 
or  earthquakes,  and  seepage  flow  [4]. 

4.2.2  Limit  Equilibrium  Analysis 

The  limit  equilibrium  analysis  of  slope  stability  is 
composed  of  several  methods  varying  in  complexity  and  applic¬ 
ability  to  particular  conditions.  The  ultimate  solution 
of  a  limit  equilibrium  analysis  is  the  determination  of 
a  factor  of  safety.  "The  factor  of  safety  is  that  factor 
by  which  the  shear  strength  parameters  may  be  reduced  in 
order  to  bring  the  slope  into  a  state  of  limiting  equilibrium 


4 


*  - 


along  a  given  slip  surface”  [5] .  This  definition  implies 
a  uniform  state  of  stresses  along  a  given  failure  plane. 

Classes  of  analysis  within  the  limit  equilibrium  analysi 
include  the  Culmann  method  (straight  line  failure  plane) , 
Bishop's  method  (circular  arc  failure  plane),  the  logarithmic 
failure  plane  method,  and  the  irregular  failure  plane  method 
(Janbu,  1954) .  Each  method  introduces  an  additional  degree 
of  complexity;  subsequently,  the  resulting  factor  of  safety 
is  hoped  to  reflect  the  actual  conditions  better  than  its 
predecessor  [4] . 

4.2.3  Bishop's  Method 

The  slice  method,  as  originally  proposed  by  Fellenius 
(1927),  is  the  basis  of  Bishop's  method  of  slope  stability 
analysis.  Bishop's  method  assumes  a  circular  failure  surface 
Circles  centers  and  radii  are  varied  until  the  minimum  factor 
of  safety  is  calculated.  Figure  4.2  illustrates  the  safety 
factor  contours  created  by  iteratively  solving  the  stability 
analysis.  The  circle  center  is  the  minimum  factor  of  safety 
of  the  slope. 


Figure  4.2.  Slope  Stability  Safety  Factor  Contours  [2], 


105 


For  each  slice,  the  forces  are  evaluated  according 
to  the  limit  equilibrium  of  the  slice.  Figure  4.3  illustrates 
the  forces  considered  to  act  on  a  slice. 


Figure  4.3.  Forces  Acting  on  a  Single  Slice  [1]. 

where , 

b  =  slice  width 

a  =  slice  inclination 

P  =  normal  force  on  slice  bottom 

W  =  slice  weight 

S  =  shear  force  on  slice  bottom 

E  -  normal  force  on  slice  side 
n 

T  =  shear  force  on  slice  side 
n 

The  equilibrium  of  the  entire  soil  mass  is  evaluated 
by  summing  the  forces  acting  on  all  the  slices.  Because 
the  factor  of  safety  appears  on  both  sides  of  the  limiting 


equilibrium  equation,  iterations  are  required  for  each  failure 
surface.  The  equation  is  as  follows: 


p  _  z  (Cb  +  Ptan<j) )  sec  a 
.  tan  4>  tan  a 

_ _ F _ 

ZW  sin  a 

where , 

<j>  =  angle  of  internal  friction 
C  =  cohesion 
a  =  slice  inclination 
F  =  factor  of  safety 

The  user  can  find  the  derivation  of  this  equation  in 
reference  5,  p.  161. 

Earthquake  forces  treated  as  equivalent  static  forces 
are  added  as  additional  driving  forces.  The  force,  F,  is 
calculated  by  accelerating  the  slice  mass. 

F  =  ma 


therefore , 


where , 

F  =  earthquake  force 
m  =  slice  mass 
W  =  slice  weight 
g  =  acceleration  of  gravity 
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k  =  seismic  coefficient 
eq 

The  earthquake  force  may  reduce  the  normal  force  acting 
on  the  bottom  of  the  slice  [3]  . 

4 . 3  Program  Use  and  Limitations 

4.3.1  General 

BISHOP1  was  translated  into  Applesoft  Basic  Language 
by  the  author.  The  original  slope  stability  program  was 
programmed  in  TRS-80  Basic  language  and  was  presented  in 
reference  3  by  Cross.  Since  the  author  did  not  program 
the  software,  the  programming  rationale  for  selection  of 
particular  program  flow  paths  is  not  presented.  During 
translation  of  the  program,  all  print  and  input  statements 
were  modified;  additionally,  the  program  was  altered  to 
calculate  a  rapid  drawdown  without  repetitive  data  input. 

The  program  is  user  oriented.  All  input  is  prompted 
by  brief,  concise  statements  and  questions.  As  with  any 
program,  the  user  should  be  familiar  with  the  associated 
theory,  required  input  data,  and  the  sequence  of  input. 

The  user  is  urged  to  tabulate  and  check  the  input  data  prior 
to  running  the  program. 

4.3.2  Data  Input 

As  depicted  in  Figure  4.4,  the  geometry  of  the  physical 
conditions  is  specified  by  a  series  of  point  and  line  numbers. 
The  user  should  graphically  display  the  problem  using  an 
appropriate  scale  such  that  the  entire  problem  is  in  the 
first  quadrant.  Points  on  the  top  soil  line  must  be  numbered 


A  a'aV.  V_  v  ■ 


(ZD" Point  # 

/\~  Line  # 


Horizontal  Distance 


Figure  4.4.  Input  Parameters  for  Slope  Stability  Analysis. 


from  left  to  right.  The  user  must  insure  that  the  extreme 
left  points  are  outside  of  the  specified  failure  circle; 
similarly,  the  extreme  right  points  must  be  outside  the 
failure  circle.  Below  the  top  soil  line,  points  may  be 
numbered  at  random,  but  the  interpretation  of  output  data 
is  simplified  by  maintaining  a  lef t-to-right  numbering 
order.  Point  numbers  are  specified  by  their  x  and  y  coor¬ 
dinates  . 

Lines  are  specified  by  the  left  and  right  point  numbers 


Vertical  lines  are  not  permitted.  If  a  vertical  line  is 


required,  the  user  may  slope  the  line  by  an  insignificant 
amount;  i.e.,  .01.  Interpretation  of  the  output  data  is 
simplified  by  maintaining  a  consistent  lef t-to-right 
numbering  order,  but  it  is  not  required. 

Soil  types  are  specified  by  number.  While  inputting 
line  data,  the  user  must  specify  which  soil  is  beneath  the 
particular  line.  Soil  physical  properties  are  input  as 
the  angle  of  internal  friction,  cohesion,  the  unit  weight, 
and  an  indication  if  the  soil  is  saturated  (0  =  yes,  1  =  no). 
Saturation  implies  the  soil  is  below  the  groundwater  table; 
this  is  the  method  used  to  specify  groundwater  table  location 

The  failure  circle  is  specified  by  inputting  the  x 
and  y  coordinates  of  the  circle  center  and  the  radius. 

The  pool  or  reservoir  is  specified  by  the  pool  surface 
elevation  and  the  extreme  left  and  right  x-coordinates  of 
the  line  defining  the  pool  water  level.  The  piezometric 
surface  within  the  soil  mass  is  identified  by  indicating 
a  saturated  soil  while  inputting  the  soil  data;  therefore, 
the  phreatic  surface  is  a  series  of  lines  with  a  wet  soil 
above  and  a  saturated  soil  beneath. 

Any  units  may  be  used  but  must  be  consistent;  i.e., 
feet,  pounds,  and  degrees. 

The  sequence  of  input  is  as  follows: 

a)  problem  heading 

b)  pool  elevation 
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c)  extreme  left  x-coordinate  of  pool 

d)  extreme  right  x-coordinate  of  pool 

e)  water  unit  weight 

f)  earthquake  seismic  coefficient,  k 

g)  total  number  of  points  to  be  specified  for  all 
lines  (see  Section  4.3.5,  Limitations).  Note  that 
the  pool  water  level  does  not  have  line  points 
associated  with  it. 

h)  point  #1,  x-coordinate 

i)  point  #1,  y-coordinate 

j)  point  #i,  x-coordinate 

k)  point  #i,  y-coordinate 

l)  number  of  lines  (see  Section  4.3.5) 

m)  left  point,  line  #1 

n)  right  point,  line  #1 

o)  soil  beneath  line  #1 

p)  left  point,  line  #i 

q)  right  point,  line  #i 

r)  soil  beneath  line  #i 

s)  number  of  soil  types  (see  Section  4.3.5) 

t)  unit  weight,  soil  #1 

u)  cohesion,  soil  #1 

v)  phi  angle,  soil  #1 

w)  is  soil  saturated,  0  =  yes,  1  =  no 

x)  unit  weight,  soil  #i 

y)  cohesion,  soil  #i 


z)  phi  angle,  soil  #i 
aa)  is  soil  saturated,  0  =  yes,  1  =  no 

bb)  minimum  number  of  slices  to  be  used  (see  Section  4.3.3) 

cc)  x-coordinate  of  circle  center 

dd)  y-coordinate  of  circle  center 

ee)  failure  circle  radius 

ff)  drawdown  option  (see  Section  4.3.3) 

See  Section  4.5  for  a  printout  of  the  above  statements 
and  associated  input  for  an  example  problem.  As  previously 
mentioned,  the  user  should  create  a  table  of  input  parameters 
prior  to  running  BISH0P1.  Example  Problem  #1  has  approxi¬ 
mately  one  hundred  input  variables. 

4.3.3  Options 

BISHOPl  can  calculate  a  rapid  drawdown  condition  without 
repetitive  input  of  the  problems  physical  geometry.  This 
option  is  prompted  by  a  question  after  the  initial  input 
of  data  and  again  following  the  output  of  calculated  data. 

If  the  user  desires  to  calculate  the  safety  factor  for  only 
a  drawdown  condition,  a  pool  elevation  of  C  should  be  input. 

The  rapid  drawdown  option  sets  the  pool  elevation  equal 
to  zero  and  does  not  affect  the  phreatic  surface. 

In  most  design  studies,  both  the  steady  state  and  rapid 
drawdown  conditions  are  analyzed.  To  analyze  both  cases, 
first  calculate  the  steady  state  safety  factor  by  specifying 
a  pool  elevation.  At  the  completion  of  this  analysis,  the 
computer  asks  if  a  rapid  drawdown  analysis  is  desired.  At 


that  time,  specify  a  rapid  drawdown  condition  and  calculate 
the  slope  stability.  This  sequence  provides  the  most  infor¬ 
mation  with  the  least  required  user  input. 

A  minimum  number  of  slices  can  be  specified  (the  computer 
may  use  more  slices  in  some  instances).  The  original  program 
automatically  sets  the  minimum  number  of  slices  at  10  [3]. 

The  program  was  modified  to  allow  the  user  to  specify  a 
minimum  number.  The  maximum  number  is  set  at  25  slices 
due  to  the  memory  constraints  of  the  hardware.  The  author 
has  successfully  used  25  slices,  but  in  approximately  50% 
of  the  cases,  an  "Out  of  Memory"  error  was  encountered. 

The  number  of  slices  specified  is  contingent  on  the  number 
of  points,  lines,  and  soil  types  input;  the  user  must  use 
trial  and  error  to  best  suit  the  minimum  number  of  slices 
to  the  physical  geometry  of  the  particular  problem.  The 
author  never  encountered  an  "Out  of  Memory"  error  when  10 
slices  were  specified.  Fewer  than  10  slices  should  not 
be  used  for  most  problems. 

The  key  advantage  of  BISHOP1  is  the  user  option  to 
specify  various  failure  circles  without  re-entering  physical 
geometry  data.  At  the  end  of  each  run,  the  user  is  prompted 
to  define  another  failure  circle.  After  defining  the  new 
circle,  the  program  repeats  the  stability  analysis  using 
the  physical  geometry  data  previously  input.  The  rapid 
drawdown  option  is  used  concurrently  with  the  above  option. 

As  depicted  in  Figure  3.2,  the  user  can  manually  "search" 


for  the  minimum  factor  of  safety  by  establishing  contours 
based  on  successive  runs. 

4.3.4  Output 

Output  of  calculated  data  is  composed  of  three  general 
types.  The  first  is  automatically  printed  at  the  end  of 
each  run  and  consists  of  the  calculated  factor  of  safety, 
circle  definition,  and  earthquake  loading  factor.  The  second 
section  of  output  is  at  the  option  of  the  user.  When  speci¬ 
fied,  a  "formal  printout"  of  all  input  data  will  be  printed. 
When  successive  runs  are  made  while  searching  out  the  minimum 
factor  of  safety,  the  "formal  printout"  should  not  be  speci¬ 
fied  until  after  the  last  run.  The  "formal  printout"  is 
used  to  document  the  output  data,  thus  enabling  the  user 
to  identify  a  particular  series  of  circle  definitions  and 
factors  of  safety  to  a  particular  problem. 

The  third  option  is  a  list  of  variables  calculated 
during  a  particular  stability  analysis.  Each  diagnostic 
list  will  be  different  for  each  circle  definition  or  physical 
geometry.  The  diagnostic  option  creates  a  list  of  slices 
and  the  corresponding  weights,  inclinations,  cohesions, 
widths,  effective  weights,  phi  angles,  and  lever  arms  (X) 
used  in  the  stability  analysis.  Additionally,  the  factor 
of  safety  iterations  are  listed  and  the  net  force  required 
to  provide  a  factor  of  safety  equal  to  one.  See  Section 
4.5  for  a  printout  of  the  diagnostic  list. 
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4.3.5  Limitations 

The  number  of  points,  lines,  soil  types,  and  number 
of  slices  are  limited.  To  avoid  program  errors,  the  user 
must  not  specify  more  variables  than  allowed  by  Table  4.1. 


Table  4.1.  Input  Limitations. 


Parameter 

Maximum 

Points 

20 

Lines 

20 

Soil  types 

5 

Slices 

25* 

*Refer  to  Section  4.3.3. 

The  physical  geometry  of  the  problem  must  be  specified 
completely  within  the  first  quadrant.  No  negative  values 
of  x  and  y  coordinates  are  allowed.  No  vertical  lines  may 
be  specified.  As  previously  discussed,  offsetting  the 
x-coordinate  of  a  vertical  line  by  an  insignificant  amount 
will  suffice  to  avoid  this  limitation. 

A  slope  must  be  specified  to  move  up  and  out  away  from 
the  coordinate  origin  as  depicted  in  Figure  4.1.  Erroneous 
slope  stabilities  were  calculated  when  this  rule  was  not 
followed.  By  rotating  the  scope  geometry  180°  such  that 
the  embankment  sloped  up  and  out  away  from  the  origin,  a 


correct  factor  of  safety  was  obtained.  This  limitation 
was  not  mentioned  by  Cross  in  his  presentation  [3] . 

A  limitation  of  BISH0P1  is  the  lack  of  a  search  routine 
which  would  yield  a  minimum  factor  of  safety  without  user 
intervention.  Main  frame  computing  facilities  have  programs 
available  to  perform  this  task;  for  example,  ICES  LEASE  I. 

Once  again,  this  ability  is  dependent  upon  memory  requirements 
beyond  the  capacity  of  most  personal  computers. 

4.3.6  Warnings 

In  Section  4.3.2,  the  user  was  advised  to  ensure  that 
the  extreme  left  and  right  points  were  outside  the  failure 
circle.  If  the  user  ignored  this,  the  warning,  "Circle 
exceeds  top  line  end  points"  will  be  printed.  The  user 
will  then  be  prompted  to  input  another  failure  circle.  If 
the  failure  circle  previously  specified  must  be  used,  the 
point  x-coordinates  must  be  altered  such  that  the  failure 
circle  is  within  the  extreme  left  or  right  points.  The 
user  must  restart  the  program  to  redefine  the  problem. 

When  the  user  defines  a  circle  with  a  radius  not  large 
enough  to  intercept  the  slope  the  warning,  "Circle  does 
not  intercept  slope"  will  be  printed.  The  user  will  be 
prompted  to  define  a  new  failure  circle.  Input  of  physical 
geometry  is  not  required.  See  Section  4.5  for  an  example 
of  this  warning. 

The  final  warning  advises  the  user  that  more  than  ten 
iterations  are  required  to  calculate  the  factor  of  safety. 


No  user  intervention  is  required  as  the  program  will  continue 
to  iterate  the  factor  of  safety.  If  the  warning  continues 
to  appear,  the  preset  iteration  tolerance  may  be  changed 
by  altering  line  3860.  The  preset  tolerance  is  .005.  If 
the  user  is  experiencing  repetitive  warnings,  try  changing 
.005  to  .01  to  avoid  this  problem.  More  than  ten  iterations 
were  never  required  to  calculate  the  factor  of  safety  during 
verification  of  this  program. 
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5  SPEEO=  ISO 

15  PRINT  "^SIMPLIFIED  BISHOP  SLOPE** " 

20  PRINT  "**  STABILITY  ANALYSIS  **" 

25  PRINT  "***************************" 

27  PRINT  :  PRINT  :  PRINT  :  PRINT  : 

30  PRINT  “  DANA  K.  EDDY,  578-80-8378" 

35  PRINT  "  GEORGIA  INSTITUTE  OF  TECHNOLOGY* " 

40  PRINT  "  DEPT.  OF  CIUIL  ENGINEERING" 

45  PRINT  "  GEOTECHNICAL  ENGINEERING  DIU. " 

47  PRINT  :  PRINT  :  PRINT 

50  PRINT  "  PROGRAM  DATE:  JUNE  1983" 

55  PRINT  "  SVSTEM  HARDWARE:  APPLE  II  PLUS,64K" 

80  PRINT  "  SVSTEM  SOFTWARE:  DOS  3.3,  APPLESOFT" 

62  PRINT  :  PRINT  :  PRINT 

65  PRINT  "  THIS  PROGRAM  CALCULATES  THE  FACTOR  OF  SAFETY'  OF  AN  EARTH  SLOP 
E  AGAINST  A  CIRCULAR  FAILURE.  THE  SIMPLIFIED  BISHOP  SLOPE  STABILITY 
ANALYSIS  IS  USED." 

67  PRINT  :  PRINT  :  PRINT 
70  SPEE0=  255 

140  DIM  P( 20,2  ),L< 20,3  ),S2< 5,4),A< 50  >,F< 50, 14  ),Z< 50,8 > 

155  PI  =  3.14159 
160-  J6  =  0 

170  REM  *** INPUT  PROGRAM  VARIABLES*** 

180  PRINT  "PROBLEM  HEADING" 

190  INPUT  H$ 

195  PRINT 

2O0  PRINT  "SUBMERGENCE  ELEVATION  (0  IF  NO  SUBMERGENCE  >" 

210  INPUT  SO 
215  PRINT 

220  PRINT  "FROM  X-CGORDINATE" 

230  INPUT  S6 
235  PRINT 

240  PRINT  "TO  ^-COORDINATE" 

250  INPUT  S7 
255  PRINT 

260  PRINT  "WATER  UNIT  WEIGHT" 

270  INPUT  W0 
275  PRINT 

280  PRINT  "EARTHQUAKE  LOADING  FACTOR" 

290  INPUT  El 

295  PRINT  :  PRINT 

300  PRINT  "NUMBER  OF  POINTS" 

310  INPUT  PI 

315  PRINT 

320  FOR  I  =  1  TO  PI 

330  PRINT  "POINT  #"I 

340  PRINT  "X-COOROINATE" 

350  INPUT  PC  1,1 > 

355  PRINT 
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360 

PRINT  "Y-C00RDINATE" 

•  ' 

370 

INPUT  P< I #2 ) 

375 

PRINT 

380 

NEXT  I 

385* 

PRINT 

390 

PRINT  "NUMBER  OF  LINES" 

400 

INPUT  LI 

405 

PRINT 

410 

FOR  I  =  1  TO  LI 

420 

PRINT  "LINE  #"I 

425 

PRINT 

430 

440 

445 

450 

460 

PRINT  "LEFT  POINT" 

INPUT  L< I » 1 > 

PRINT 

PRINT  "RIGHT  POINT" 

INPUT  L< I #2  ) 

wm 

465 

PRINT 

470 

PRINT  "SOIL  BENEATH  LINE  #"I 

teem 

480 

INPUT  L<I#3> 

485 

PRINT 

490 

NEXT  I 

495 

PRINT  :  PRINT 

500 

PRINT  "NUMBER  OF  SOIL  TYPES" 

1**4 

510 

INPUT  SI 

w  --tM 
KrtW 

515 

PRINT 

*.  4  j 

520 

FOR  I  =  1  TO  SI 

530 

PRINT  "SOIL  #"I 

•  ‘  -t' 

c--?cr 

J-jJ 

PRINT 

540 

PRINT  "UNIT  HEIGHT" 

-  \ 

550 

INPUT  S2< I  #  1  ) 

555 

PRINT 

\ 

% 

560 

PRINT  "COHESION" 

570 

INPUT  S2< I # 2  > 

kV 
'  s  \ 

575 

PRINT 

580 

PRINT  "PHI  ANGLE" 

,V  *S 

590 

INPUT  S2< I  #3) 

595 

PRINT 

“A 

600 

PRINT  "IS  SOIL  SATURATED#  0=VES, 

-NO  ?" 

>> 

A 

610 

INPUT  S2< I # 4  > 

615 

PRINT 

V-Vv 

620 

NEXT  I 

.*•  A 

625 

PRINT  :  PRINT 

627 

PRINT  "SPECIFY  THE  MINIMUM  NUMBER 

OF  SLICES  TO  BE  USED  <  MAX  =  25  ,'." 

:  vi 

628 

INPUT  S3 

629 

PRINT  :  PRINT 

630 

REM 

640 

F3  =  0 

650 

PRINT  "FAILURE  CIRCLE  DEFINITION" 

+-0I 

655 

PRINT  :  PRINT 

-  -  rt 

PRINT  "X-COORDINATE  OF  CENTER" 

\  A ,  > 

s'.\ 


67©  INPUT  X 
675  PRINT 

680  PRINT  " V-C0ORD I NATE  OF  CENTER" 

690  INPUT  V 
695  .  PRINT 

700  PRINT  "CIRCLE  RADIUS" 

710  INPUT  R 

712  PRINT  :  PRINT 

713  PRINT  "FOR  A  RAPID  DRAWDOWN  TVPE  <0>*  FOR  SUBMERGENCE  ELEUAT I ON  hS  P 
EUIOUSLV  SPECIFIED  TVPE  Cl>" 

714  INPUT  ZZ:  IF  ZZ  =  1  THEN  716 

715  SO  =  0 

716  PRINT  :  PRINT 

720  REM  **CHECK,  CIRCLE  EXCEEDS  TO  LINE  END  POINTS** 

730  U1  =  PI 

740  FOR  I  =  2  TO  PI 

750  IF  P<  1,1)  <  PCI  -  1,1)  AND  U1  =  PI  THEN  770 
760  GOTO  78 U 
770  U 1  =  I  -  1 
780  NEXT  I 

790  J1  =  R  *  R  ~  CPC  1*2)  ~~  V)  2 
80O  J2  =  R  *  R  -  C  PC  U 1 ,2 )  -  V)  *  2 
310  IF  J1  <  =  0  THEN  330 

820.  IF  J1  >  0  AND  PC  1,1)  >  X  -  SQR  CJ1)  THEN  860 
830  IF  J2  <  =  0  THEN  850 

840  IF  J2  >  0  AND  PC U 1,1)  <  X  +  SQR  C J2  >  THEN  860 
350  GOTO  380 

660  PRINT  “*  CIRCLE  EXCEEDS  TOP  LINE  END  POINTS  *" 

370  GOTO  4380 

380  REM  **DEFINE  INTERSECTION  OF  CIRCLE  WITH  LINES** 

890  FOR  I  =  1  TO  LI 
900  XI  =  PC  LC I  *  1  )  *  1 > 

910  VI  =  PCLCI,1),2) 

920  X2  =  PCLCI*2)*1) 

930  V2  =  PC  LC I *2 )*2  ) 

340  IF  X2  =  XI  THEN  960 

950  GOTO  970 

360  S  =  9.99E  +  10 

370  IF  X2  <  >  XI  THEN  990 

980  GOT  0  1 000 

990  S  =  C  V2  —  VI J  /  C  X2  —  X 1 > 

1000  IF  ABS  CS)  <  1.0E  -  5  THEN  1150 
1010  Cl  =  XI  -  VI  S 

1020  C2  =  1  /  S  A  2+1 

1030  C3  =  2  *  Cl  /  S  -  2  *  X  /  S  -  2  *  V 

104O  C4  =  Cl  A  2  -  2  *  X  *  Cl  +  X  2  +  V  2  -  R  2 

1950  C5  =  C3  2  -  4  *  C2  *  C4 

I860  IF  C5  <  0  THEN  1080 

1070  GOTO  1030 

1080  ZC I , 1 >  =  0 
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109©  IF  C5  <  0  THEN  1630 

1100  Q1  =  <  -  C3  +  SQR  <C5>>  ✓  <2  *  C2:> 

1110  Q2  =  <  -  C3  -  SQR  CC5>>  ✓  <2  *  C2> 

1120  Q3  =  Q1  s  S  +  Cl 
1130  Q4  =  02  S  +  Cl 
1140  GOTO  1240 

1150  C5  =  R  ''  2  -  <V  -  VI  >  a  2 
1160  IF  C5  <  0  THEN  1180 
1170  GOTO  1190 
1180  Z< 1*1 >  =  © 

1190  IF  C5  <  0  THEN  1630 
1200  03  =  X  +  SQR  < C5  > 

1210  Q4  =  X  -  SQR  <C5> 

1220  Q1  =  VI 

1230  Q2  =  VI 

1240  Jl  =  O 

1250  J2  =  0 

1260  IF  083  <S>  <  =  9.99E  +  9  AND  Q3  =  >  XI  AND  Q3  <  =  X2  THEN  1 

1270  GOTO  1290 
1280  Jl  =  1 

1290  IF  AB3  <S>  <  =  9.99E  +  9  AND  Q4  =  >  XI  AND  Q4  <  =  X2  THEN  1 

1300  GOTO  1320 
1310  J2  =  1 

1320  IF  S  <  -  9.99E  +  9  AND  Q1  >  =  V2  AND  Q1  <  *  VI  THEN  1340 

1330  GOTO  1350 
1340  Jl  =  1 

1350  IF  3  <  -  9.99E  +  9  AND  Q2  >  =  V2  AND  02  <  =  VI  THEN  1370 

1360  GOTO  1380 
1370  J2  =  1 

1380  IF  S  >  9.99E  +  9  AND  01  >  =  VI  AND  01  <  =  V2  THEN  1400 

1390  GOTO  1410 
1400  Jl  *  1 

1410  IF  S  >  9.99E  +  9  AND  02  >  =  VI  AND  02  <  -  V2  THEN  1430 

1420  GOTO  1440 

1430  J2  =  1 

1440  Z<I,1>  =  Jl  +  J2 

1450  IF  Jl  =  1  THEN  1470 

1460  GOTO  1480 

1470  Z<I*2>  =  03 

1480  IF  Jl  =  1  THEN  1500 

1490  GOTO  1510 

1500  Z<I»3)  =  01 

1510  IF  Jl  =  0  AND  J2  =  1  THEN  1530 
1520  GOTO  1540 
1530  Z<  I  »2  ':>  =  04 

1540  IF  Jl  =  0  AND  J2  =  1  THEN  1560 
1550  GOTO  1570 
1560  Z<I*3>  =  02 

1570  IF  Jl  =  1  AND  J2  =  1  THEN  1590 
1580  GOTO  160O 
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1 590 

I  ,  4  .•*  =  Q4 

1  600 

IF  J1  =  1  PND  J2  = 

1  THEN 

1610 

GOTO  1630 

1620 

ZC  1 ,5>  =  Q2 

1630 

NEXT  I 

1640 

X4  =  0 

1650 

X5  =  9.99E  +  20 

1 660 

11  =  1 

1670 

FOR  I  =  1  TO  LI 

1680 

IF  ZC 1,1 >  =  >1  THEN  1700 

1690 

GOTO  1710 

1  700 

pc  ii  :>  =  zc  i,2) 

1710 

IF  ZCI,1>  =  >1  THEN  1730 

1 720 

GOTO  1740 

1730 

11  =  11  +  1 

1740 

IF  ZC  I  ,  1  >  -  THEN 

1760 

1750 

GOTO  1770 

1 760 

PC  11  )  ■  ZC  1,4) 

1770 

IF  ZC 1,1  )  =  2  THEN 

1790 

1 780 

GOTO  1800 

1790 

11  =  11  +  1 

1  300 

NEXT  I 

1310  IF  II  =  1  THEN  1830 
1820  GOTO  1840 

1830  PRINT  "CIRCLE  DOES  NOT  INTERCEPT  SLOPE" 

1840  IE  II  =  1  THEN  4380 

1350  REM  **SET  UP  OF  SLICE  PP.RPV** 

1880  FOR  I  =  1  TO  II  -  1 
1370  IF  PC I >  >  X4  THEN  1390 
1880  GOTO  1900 
1890  X4  =  PCI) 

1900  IF  PC I >  <  X5  THEN  1920 
1910  GOTO  1930 
1920  X5  =  PC  I  ) 

1 930  NEXT  I 

1940  FOR  I  =  1  TO  Pi 

1950  IF  PC  1,1)  <  X4  PNO  PC  1,1)  >  X5  THEN  1970 
I960  GOTO  1980 
i  970  PC  1 1  >  -  PC  1 , 1 

1980  IF  PC  1,1)  <  X4  PND  PC  1,1 >  >  X5  THEN  2000 

1990  GOTO  2010 

2000  II  =  II  +  1 

2010  NEXT  I 

2020  II  =  II  -  1 

2030  FOR  I  =  1  TO  1 1 

2040  FOR  J  =  1  TO  II  -  1 

2050  IF  PC  J  +  1  >  >  PCJ>  THEN  2090 

2060  J1  =  PC  J  +  1> 

2070  PC  J  +  1  )  =  PC  J  > 

2080  PC  J  )  =  J1 


*  _si 

•j 


j 

j 


j 


£090 

2100 

2110 

2120 

2130' 

2140 

2150 

2160 

2170 

2130 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2230 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

2370 

2380 

2330 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 


NEXT  J 
NEXT  I 
U1  =  0 

FOR  I  =  1  TO  II  -  1 
IF  A<  I  >  <  A<  I  +  1  >  THEN  2150 
GOTO  2190 
U1  =111  +  1 

IF  A<I>  <  A<I  +  1>  THEN  2180 
GOTO  2130 
AOJ1  >  =  A<  I> 

NEXT  I 
111  =  U1  +  1 
A<U1>  =  R<I1> 

II  =  U1 

REM  **OEFINE  SLICE  BOUNDARIES** 

Q1  =  A<  II  >  -  A<  1  > 

Q2  =  Q1  ✓  S3 
U1  =  II 

FOR  I  =  1  TO  U1  -  1 
03  =  A< I  +  1  )  -  A<  I  > 

04  =  INT  033  x  02 >  +  1 
Cl  =  03  £  04 
C2  =  A<  I  > 

FOR  J  =  1  TO  04 
IF  J  <  04  THEN  2350 
GOTO  2360 
11  =  II  +  1 
IF  J  <  04  THEN  2380 
GOTO  2390 
A<  II >  =  C2  +  Cl 
IF  J  <  04  THEN  2410 
GOTO  2420 
C2  =  C2  +  Cl 
NEXT  J 
NEXT  I 

FOR  I  =  1  TO  II 
FOR  J  =  1  TO  II  -  1 
IF  AOJ  +  1>  >  A<  J )  THEN  2500 
J 1  =  A<  J  +  1  > 
fi<  J  +  1  >  =  R<  J  > 

A<  J  >  =  J 1 
NEXT  J 
NEXT  I 

REM  ++DEFINE  SOIL  PARAMETERS  FOR  EACH  SLICE 
FI  =  II  -  1 
FOR  I  =  1  TO  FI 
F1’  I  »4  =  A'.'  I  +  1  >  —  A<  I  > 

X6  =  F<  I  #  4  > 

F>;  1 ,7  ">  -  <  A(  I  +  l  )  +  H<  I  >  >  /  2 

X3  =  F< I #7 ) 
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259Q  VI  =  V  -  SQR  <R  ^  2  -  < P<  I  >  -  X)  2) 

2600  V2  =  V  -  SQR  <  R  ~  2  -  <  P<I  +  1  >  -  X  >  ~  2  > 

2610  05  =  PTN  C  PBS  <V2  -  VI  >  ✓  F<I,4>> 

2620  .  IF  V2  <  VI  THEN  2640 

2630  GOTO  2650 

2640  P5  »  -  05 

2650  F(I,2)  =  05 

2660  IF  05  =  0  THEN  2680 

2670  GOTO  2690 

2680  F<  1 ,2  >  =  1.0E  -  5 

2690  V3  =  V  -  SQR  <  R  2  -  <  X3  -  X )  A  2  > 

2700  14  =  0 

2710  FOR  J  =  1  TO  LI 

2720  L5  =  L<  J  ,  1 ) 

2730  L6  =  L<J,2) 

2740  IF  P<L5,2)  <  =  V3  PND  P<L6,2>  <  =  V3  THEN  2840 

2750  IF  P<  L5,  1  )  <  X3  0ND  P< L6,l  >  <  X3  THEN  2840 

2760  IF  P<L5,1)  >  X3  0ND  P<L6,1>  >  X3  THEN  2840 

2770  V6  =  P< L5*2 )  +  <P<L5,2>  -  P<L6,2)>  ✓  <P<L5,1>  -  P<L6,1>>  • 

,  1  )  > 

2780  IF  V6  <  =  V3  THEN  2840 

2 { 90  14  =  14  +  1 
2800.  Z<I4,1>  =  V6 
2810  2<  14,2  )  -  L<J#3> 

2320  W  =  w 
2330  E  =  0 
2340  NEXT  J 

2850  IF  14  =  1  THEN  2970 

2860  FOR  J  =  1  TO  14 

2370  FOR  J1  =  1  TO  14  -  1 

2880  IF  Z<J1,1>  =  >  Z<J1  +  1,1)  THEN  2950 

2330  L5  =  01 , 1  ) 

2900  L6  =  Z<  J1 ,2  > 

2910  Z<  J1 , 1  )  =  Z<  J1  +  M) 

2920  Z<  J 1 » 2  )  =  Z<J1  +  1  #2  > 

2930  Z<J1  +■  1,1)  =  L5 

2940  Z<J1  +  1,2)  =  LB 

2950  NEXT  J1 

2360  NEXT  J 

2970  14  =  14  +  1 

2980  Z< 14, 1 )  =  V3 

2990  FOR  J1  =  1  TO  14  -  1 

3O00  IF  I  =  1  0ND  J1  =  1  0ND  X3  =  >  S6  THEN  3020 

3O10  GOTO  3030 
3O20  16  =  S0  -  VI 

3030  IF  I  =  FI  PND  J1  =  1  PND  X3  =  >  86  PND  X3  <  =  S7  THEN 

304O  GOTO  3060 
3050  JS  =  S0  -  V2 

3060  M  =  N  +  (  Z<  J 1 , 1  )  -  Z(  J1  +  1,1  )>  ♦  X6  *  S2(Z<  J1 ,2),1  > 

307O  IF  Z<  Jl , 1 >  <  SO  PND  X3  =  >  S6  PHD  X3  <  =  S7  THEN  303O 


i’-  -■>  a  - 
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308O  GOTO  3100 

3090  W  =  N  +  (30  -  Z<  J 1 , 1  >  >  *  K6  *  WO 
3100  IF  S2<  Z(  J 1 ,  2  >  ,  4 >  >  0.95  THEN  3120 
3110  GOTO  3130 
3120  E4  =  S2(Z<  J1,2>,1> 

3130  IF  S2(Z<  Jl,2>,4>  <  0.95  THEN  3150 

3140  GOTO  3160 

3150  E4  =  S2( Z< J1 #2 >, 1 >  -  WO 

3160  E  =  E  +  <Z(Jl,i:>  -  Z(J1  +  1,1  >>  *  X6  *  E4 
317®  NEXT  J1 
3180  F( I , 1  )  =  W 
3190  F< I  ,5  >  =  E 

320O  F<  1,3  >  =  S2<  Z<  1 4  -  1,2),2> 

3210  F<  1 ,6  >  =  2  *  PI  *  < S2< Z(  1 4  -  1,2>,3>  /  360; 

3220  NEXT  I 

3230  IF  F9  =  O  THEN  3360 

3235  L*  =  CHRX  (4):  PRINT  LX;"PR#1'* 

3240  PRINT  "SLICE  WEIGHT  INCLINATION  COHESION  WIDTH  EFF  WEIGH 

T  PHI  X" 

3241  PRINT 

3280  0  =  360  /  (2  *  PI  •' 

3290  FOR  I  =  1  TO  FI 

3300  POKE  36,3:  PRINT  Ii:  POKE  36,7:  PRINT  F< 1,1 POKE  36,19:  PRINT  F< I 
,2)  *  0;:  POKE  36,34:  PRINT  F(I,3>;:  POKE  36,41:  PRINT  F<I,4)s:  POKE 
36,53:  PRINT  F(I,5>;:  POKE  36,65:  PRINT  F(I,6>  *  0;:  POKE  36,69:  PRINT 


f<  i  ,7  :> 

3340 

NEXT  1 

3350 

PRINT  Lii "PR#0" 

3360 

D  =  O 

3370 

FOR  I  =  1  TO  FI 

3 380 

D  =  0  +  Ft.  1,1  .;  *  SIN 

(  PBS  (  F(  1,2.;  >  .'■<  *  (  F(  1 ,2  )  / 

3390 

D  =  D  +  El  *  F<  1,1)  * 

COS  (  PBS  <F(  I,2>>> 

3400 

NEXT  I 

3410 

IF  16  >  O  THEN  3430 

3420 

GOTO  3440 

3430 

17  —  WO  *  IS*  16  *  ( R 

-  16  /  3 ( 2  *  R  > 

3440 

IF  16  >  O  THEN  3460 

3450 

GOTO  3470 

3460 

D  =  D  -  SGN  <D>  *  17 

3470 

IF  16  >  O  PND  F9  =  1 

THEN  3485 

3480 

GOTO  3510 

3485 

LX  =  CHRX  (4):  PRINT 

LX ; "PR#1 " 

3487 

PRINT  :  PRINT 

3490 

PRINT  "DRIUING  FORCE 

COUNTER  BALANCE  OF  "17"#." 

3  *4  h  *5 

PRINT  :  PRINT 

7300 

PRINT  LX i "PR#0" 

351 0 

IF  J6  0  THEN  3530 

35J.0 

GOTO  3540 

}  S  "70 

17  =  NO  ♦  JP  ♦  JG  ♦  *.  F 

!  -  J6  3  >  •■■■  <  2  *  R 

3  j40 

IF  J6  O  THEN  3560 

550  GOTO  3570 
560  D  =  0  +  SGN  <10}  17 

570  IF  J6  >  0  AND  F9  =  1  THEN  3585 
580  GOTO  3610 

585  Li  =  CHRi  < 4 >:  PRINT  Li^PR#!" 

587  PRINT  :  PRINT 

590  PRINT  "0RIUIN6  FORCE  INCREASE  OF  "17"#." 

595  PRINT  :  PRINT 
800  PRINT  Li; "PR#0" 

BIO  REM  **ITERATIUE  SOLUTION  FOR  FACTOR  OF  SAFETY** 

620  FO  =  1 

630  R4  =  0 

640  16  =  0 

650  FOR  I  =  1  TO  FI 

660  R1  =  F'l  I  »3  >  *  F<  1*4}  +  F<I#5}  *  TAN  <F<I#8}) 

670  R2  =  1  COS  i.  ABS  'I  F  >!  1*2  )  }  3 

680  R3  =  1  +  TAN  CF<I,6  3>  *  TAN  <  F<  1*233  F0 


690 

F:4  =  R4  +  R1  *  •: 

:R2  ✓  R3> 

700 

NEXT  I 

710 

F2  =  R4  ✓  D 

720 

16=  16+  1 

730 

IF  F9  =  1  THEN 

3750 

740 

hijTO  3820 

750. 

IF  16  =  1  THEN 

3765 

760 

GOTO  3775 

765 

Li  =  CHRi  <4  3: 

PRINT  Li; "PR#!" 

770 

PR I NT  "I TERAT I ON " * " INITIAL" * " CALCULATED " 

775 

Li  =  CHRi  (4): 

PRINT  Li; "PR#1 " 

780 

PRINT  TAB<  5 3I6*F0*F2 

300 

PRINT  Li; "PR#Q 

II 

820 

IF  16  >  10  THEN  3846 

830 

GOTO  3850 

840  PRINT  "WILL  HOT  CLOSE" 

8:50  IF  16  >  10  THEN  3970 

■860  IF  ABS  C  ABS  < F0 >  -  ABS  <  F2  >  >  <  0.005  THEN  3910 

870  F0  =  ABS  <F2  3 

1880  R4  =  0 

'890  GOTO  3650 

:895  PRINT  :  PRINT 

;910  Li  =  CHRi  ( 4  3 :  PRINT  Li;  "PR  #1":  PRINT 

:915  PRINT  Hi 

■'316  PRINT  " ****•»■***•*•**•*■•*•*••*•*#•*••*■*••*■•*•*••«•**•*•»•*•*■•*■•*••*• " 

:920  PRINT  "FACTOR  OF  SAFETY  =  "F2;"  AT  X  =  "N;"  Y  =" 
>935  PRINT 

’940  PRINT  "EARTHQUAKE  LOADING  FACTOR  ="E1 
5955  PRINT  :  PRINT 
■360  PRINT  Li; "PR#0" 

5965  HOME 

>970  PRINT  "DO  YOU  WANT  A  FORMAL  PRINTOUT  V  OR  N  3" " 
>980  INPUT  Ai 


PR  #1":  PRINT  :  PRINT 

****•*••#•*#•»  •*•*•*■*■*•*  " 


PRINT  "FACTOR  OF  SAFETY  =  "F2;"  AT  X  ="X;"  Y  ="Y;"  RADI 
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3*385  PRINT  :  PRINT 

3930  IF  AT  =  "N"  THEN  4326 

4905  LT  =  CHRT  <  4  >:  PRINT  LT;"PR#1" 

4007  PRINT  HT 
4003  PRINT  :  PRINT 

4010  PRINT  "HATER  UNIT  HEIGHT  =“W0;“  EARTHQUAKE  LOADING  FACTOR  =  "E1 
4020  PRINT  "SUBMERGENCE  ELEU  ^'SQ;**  FROM  X="SS;"  TO  X="S7 
4030  PRINT 

4040  PR  I  NT  "  PO I  NT  "  *  "  X-CQORD I  NATE  "  *  "  V-COORD I  NOTE 11 

4050  FOR  I  =  1  TO  PI 

4060  PRINT  I, PC  1,1 >,P< I,2> 

4070  NEXT  I 

4080  PRINT  :  PRINT 

4030  PRINT  "LINE";"  LEFT  PT";"  RIGHT  PT" ;"  SOIL" 

4100  FOR  I  =  1  TO  LI 

4110  POKE  36*2:  PRINT  I;:  POKE  36*10:  PRINT  LCI*1>;:  POKE  36*13:  PRINT  L 
1*2!);:  POKE  36*27:  PRINT  LCI*3> 

4120  NEXT  I 
4130  PRINT  :  PRINT 

4140  PRINT  "SOIL";"  UNIT  HEIGHT";"  COHESION";"  PHI";"  SAT'D" 

4150  FOR  I  =  1  TO  SI 

4160-  POKE  36*2:  PRINT  I;:  POKE  36*10:  PRINT  S2CI*1>;:  POKE  36*21:  PRINT 
2C 1 ,2  >*:  POKE  36*30:  PRINT  S2C 1 ,3 >; :  POKE  36*37:  PRINT  S2CI*4> 

4170  NEXT  I 
4180  PRINT  :  PRINT 

4130  PRINT  "CIRCLE:  X="X;"  V="V;"  RADIUS  ="R;"  FACTOR  OF  SAFETV  ="F2 

42O0  PRINT  :  PRINT 
4210  PRINT  LT; "PR#0" 

4215  HOME 

4320  PRINT  "DO  VOIJ  HANT  A  DIAGNOSTIC  RUN  <  V  OR  NT?" 

4330  INPUT  AT 

a732  PRINT  :  PRINT 

4340  IF  AT  =  "V"  THEN  4360 

4350  GOTO  4370 

4360  F3  =  1 

4370  IF  AT  =  "V"  THEN  720 

4380  PRINT  "DO  VOU  HANT  TO  DEFINE  ANOTHER  FAILURE  CIRCLE  CV  OR  N>?" 

4330  INPUT  AT 

4331  PRINT  :  PRINT 

4400  IF  AT  =  "V"  THEN  630 

4410  PRINT  "BVE-BVE" 
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4 . 5  Program  venticatxon 
4.5.1  Problem  #1 


**SIMPLIFIED  BISHOP  SLOPE-** 
**  STABILITY  analysis  ** 


i  *-*"*■■*•*■■**■*■•*■*•■*• ■> 


DANA  K.  EODV,  578-80-8373 
GEORGIA  INSTITUTE  OF  TECHNOLOGV 
DEPT.  OF  CIVIL  ENGINEERING 
GEOTECHNICAL  ENGINEERING  DIO. 


PROGRAM  DATE:  JUNE  1983 

SYSTEM  HARDWARE:  APPLE  II  PLUS, 64 K 

SYSTEM  SOFTWARE:  DOS  3.3,  APPLESOFT 


THIS  PROGRAM  CALCULATES  THE  FACTOR  OF  SAFETY  OF  AH  EARTH  SLOPE  AGAINST  A 
LAR  FAILURE.  THE  SIMPLIFIED  BISHOP  SLOPE  STABILITY  ANALYSIS  IS  USED. 


PROBLEM  HEADING 
?AN  EXAMPLE  PROBLEM 

SUBMERGENCE  ELEVATION  <0  IF  NO  SUBMERGENCE > 
?  1 040 

FROM  X-COORDINATE 

?8 

TO  X-COORDINATE 
7132 

WATER  UNIT  WEIGHT 
762.4 

EARTHQUAKE  LOADING  FACTOR 

70 


NUMBER  OF  POINTS 

713 

POINT  #1 

X-COORDINATE 

'0 

^-COORDINATE 

71020 


POINT  #2 
^-COORDINATE 
7  1  O0 


V-C00RD I NOTE 
?1020 

POINT  #3 
X-COORDINATE 

7132 

^-COORDINATE 

71040 

POINT  #4 

X-COORDINATE 

7162.5 

V-C00RDINATE 

71045 


POINT  #5 

X-COORDINATE 

7172.5 


V-C00RDINATE 

71045 

POINT  #6 
X-COORDINATE 


t'-C00ROINATE 

71020 

POINT  #7 
X-COORDINATE 

7  O 

t'-COOROINATE 

71000 


POINT  #3 
X-COORDINATE 

7  300 


V'-COORDINATE 
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POINT  #9 
X-COORDINATE 

7138 

V-COOROINATE 

?1035 

POINT  #10 
^-COORDINATE 

7 1 50 

V-C00RD I NATE 
71030 

POINT  #11 

X-COORDINATE 

7161 

V-C00ROINATE 

71027.5 

POINT  #1£ 

X-COORDINATE 

7168 

^-COORDINATE 

71025 

POINT  #13 

X-COORDINATE 

7171 

V-C00RD I NATE 
71020 


NUMBER  OF  LINES 
713 

LINE  #l 

left  point 

T  i 

RIGHT  POINT 


SOIL  BtNtATH  LINE  #1 


LINE  #2 


LEFT  POINT 


RIGHT  POINT 

73 

SOIL  BENEATH  LINE  #2 

LINE  #3 

LEFT  POINT 
73 

RIGHT  POINT 
?4 

SOIL  BENEATH  LINE  #3 
?1 

LINE  #4 

LEFT  POINT 
?4 

RIGHT  POINT 
?5 

SOIL  BENEATH  LINE  #4 
71 

LINE  #5 

LEFT  POINT 
?5 

RIGHT  POINT 
76 

SOIL  BENEATH  LINE  #5 
71 

LINE  #6 

LEFT  POINT 

73 

RIGHT  POINT 
79 

SOIL  BENEATH  LINE  #6 


LINE  #7 


LEFT  POINT 

79 

RIGHT  POINT 
'TO 

SOIL  BENEATH  LINE  #7 


LINE  #3 

LEFT  POINT 

?  1 0 

RIGHT  POINT 
71 1 

SOIL  BENEATH  LINE  #3 


LINE  #9 

LEFT  POINT 
711 

RIGHT  POINT 
?12 

SOIL  BENEATH  LINE  #3 


LINE  #10 

LEFT  POINT 

?12 

RIGHT  POINT 
713 

SOIL  BENEATH  LINE  #10 


LINE  #11 
LEFT  POINT 


RIGHT  POINT 
713 


SOIL  BENEATH  LINE  #11 


LINE  #12 


LEFT  POINT 

713 

RIGHT  POINT 
76 

SOIL  BENEATH  LINE  #12 

73 

LINE  #13 

LEFT  POINT 

!  »** 

RIGHT  POINT 

;*  y 

SOIL  BENEATH  LINE  #13 

74 


NUMBER  OF  SOIL  TVPES 
74 

SOIL  #1 

UNIT  WEIGHT 
■•115 

COHESION 

'200 

PHI  ANGLE 
726 

IS  SOIL  SATURATED,  0=VES,  1=N0 
*1 

SOIL  #2 

UNIT  WEIGHT 
7 11 7 

COHESION 

'200 

PHI  ANGLE 
726 


IS  SOIL  SATURATED,  0=VES»  1 =Nu 
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SOIL  #3 

UNIT  WEIGHT 

7120 

COHESION 

7100 


PHI  ANGLE 


IS  SOIL  SATURATED,  0=VES,  l=NO  7 

70 

SOIL  #4 

UNIT  WEIGHT 
7150 

COHESION 
7 1 0000 

PHI  ANGLE 

745 

IS  SOIL  SATURATED#  G*VES#  l=NO  7 
70 


SPECIFV  THE  MINIMUM  NUMBER  OF  SLICES  TO  BE  USED  < MAX  =  25). 
710 


FAILURE  CIRCLE  DEFINITION 


X-COOROINATE  OF  CENTER 

7 1 03 .  7 

V-CQORDINATE  OF  CENTER 
7 1 100.3 

CIRCLE  RADIUS 
731.0 


FOR  A  RAPID  DRAWDOWN  TVPE  <  0  >#  FOR  SUBMERGENCE  ELEUATION  AS  PREVIOUSLY  SPECIFIED,..^ 

tvpe  <  i )  ' 

7  j. 


i 
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AN  EXAMPLE  PROBLEM 

it ■£ ■£  ■£ ■*■ -fc * •*■  •*■£•*“ •£ it -£ •#• #■ $■• ■£ & ■£ ■$ ■£  •£ 

FACTOR  OF  SAFETY  =  2.64916612  AT  X  =108.7  V  =1199.3  RADIUS  =81 
EARTHQUAKE  LOADING  FACTOR  =9 


DO  YOU  WANT  A  FORMAL  PRINTOUT  < Y  OR  N>? 
?V 


AN  EXAMPLE  PROBLEM 


HATER  UNIT  HEIGHT  =82.4  EARTHQUAKE  LOADING  FACTOR  =9 
SUBMERGENCE  ELEU  =1949  FROM  X=9  TO  X=132 


POINT 

X-COORDINATE 

V-COORO I NOTE 

1 

9 

1929 

J 

*— 

100 

1929 

~7 

132 

1949 

4 

162.5 

1945 

5 

172.5 

1945 

6 

235 

1929 

1* 

9 

1099 

o 

399 

1900 

3 

138 

1035 

19 

159 

1939 

11 

161 

1027.5 

12 

163 

1925 

13 

171 

1020 

L  i  NE 

LEFT  PT 

RIGHT  PT 

SOIL 

1 

1 

2 

3 

a! 

c. 

•«!» 

ul 

3 

4 

1 

4 

4 

5 

1 

er 

5 

K 

1 

•— 1 

rJ 

-■ 

-7 

;Z( 

•  J 

A- 

“7 

9 

19 

8 

10 

11 

cL 

9 

11 

12 

19 

•  12 

13 

11 

13 

"7 

13 

1.2 

s 

-• 

13 

s 

4 
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SOIL  UNIT  HEIGHT  COHESION 

PHI 

SAT'D 

1 

115  200 

26 

1 

2 

117  200 

26 

0 

1  120  100 

ilo 

0 

4 

150  10000 

45 

0 

CIRCLE: 

X= 108.7  V= 1 1 00 . 8  RADIUS 

=31  FACTOR  OF  SAFETV  =2.64016612 

DO  VOU 

?V 

WANT  A  DIAGNOSTIC  RUN 

<  V  OF 

t  N  >? 

SLICE 

HEIGHT 

INCLINATION 

COHESION 

HIDTH 

EFF  HEIGHT 

PHI 

1 

3419.83601 

-4.96531439 

200 

2.64388013 

137.423009 

26 

101. 689478 

C. 

15437. 8547 

-2. Oi 358368 

100 

5 . 68853236 

1135.39153 

28 

1 05 . 855709 

3 

16542. 1386 

2.01358368 

100 

5.63858236 

2290. 17546 

O 

111.54429 1. 

4 

10617.9847 

6.11577104 

200 

5. 37047255 

3337.01816 

26 

117.323819 

5 

11212.2684 

10.3152172 

200 

5.87047255 

429 1 . 56704 

26 

183. 19489 1 

6 

11493.2192 

14.5713242 

200 

5.37047255 

5052. 69402 

26 

189. 064764 

11408.9623 

18.9619938 

200 

6 

6425. 16599 

26 

135 

o 

9434.07119 

23.3314633 

200 

5. 43929873 

6943. 37233 

26 

140. 719648 

9 

8372.33118 

27.5671964 

200 

5. 43929673 

7433.64569 

26 

146. 158945 

10 

1556.34333 

30. 1967314 

200 

1. 12140656 

1556.34339 

26 

149. 439297 

11 

6689.7098 

32.9750275 

200 

5.5 

6689. 7098 

26 

152.75 

12 

4794.2082 

37. 7555341 

200 

5.5 

4794.2082 

26 

153.25 

13 

912.06297 

40.9187278 

200 

1.5 

912.06297 

26 

161.75 

14 

1399.49512 

44.0393237 

200 

4.91422302 

1399.49512 

26 

164.957112 

DRIVING  FORCE  COUNTER  BALANCE  OF  11202.7345#. 


ITERATION 

INITIAL 

CALCULATED 

1 

1 

2.3702 

899 

2. 3702899 

8.61 83 

3446 

3 

2. 61883446 

2 . 6387 

2113 

4 

2.63872113 

2 . 640 16612 
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AN  EXAMPLE  PROBLEM 

FACTOR  OF  SAFETY  =  2.64016612  AT  X  =108.7  Y  =1100.8  RADIUS  =81 
EARTHQUAKE  LOADING  FACTOR  =0 


DO  YOU  NANI  A  FORMAL  PRINTOUT  <Y  OR  N )? 
?N 


DO  YOU  WANT  A  DIAGNOSTIC  RUN  <Y  OR  N  >? 
?N 


DO  YOU  WANT  TO  DEFINE  ANOTHER  FAILURE  CIRCLE  <Y  OR  NX? 
'Y 


FAILURE  CIRCLE  DEFINITION 


X-CGORDINATE  OF  CENTER 
? 103. 7 

V-COORDINATE  OF  CENTER 
?  1 1 00 . 8 

CIRCLE  RADIUS 

•'  o  1 .  O 


FOR  A  RAPID  DRAWDOWN  TYPE  <0),  FOR  SUBMERGENCE  ELEUATION  AS  PREUIOIJSLY  SPECIFIE 
TYPE  < 1 > 

?0 


AN  EXAMPLE  PROBLEM 
*****■****#*«*■*•**•**•» 

FACTOR  OF  SAFETY  =  1.45417058  AT  X  =103.7  V  =1100.8  RADIUS  =81 

EARTHQUAKE  LOADING  FACTOR  =0 


DO  YOU  WANT  A  FORMAL  PRINTOUT  < V  OR  NY? 
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PH  EXAMPLE  PROBLEM 


HATER  UNIT  HEIGHT  =62.4  EARTHQUAKE  LOADING  FACTOR  =0 
SUBMERGENCE  ELEU  =0  FROM  X=0  TO  X=132 


POINT 

X-COORDINATE 

V-COORDINATE 

1 

0 

1020 

2 

100 

1 020 

132 

1040 

4 

162.5 

1045 

Cj 

172.5 

1045 

K 

235 

1020 

»' 

0 

1000 

O 

300 

1000 

9 

138 

1035 

10 

150 

1030 

11 

161 

1027.5 

12 

168 

1025 

13 

171 

1020 

LINE 

LEFT  FT 

RIGHT  PT 

SOIL 

1 

1 

o 

3 

2 

il* 

3 

3 

3 

4 

1 

4 

4 

5 

1 

5 

5 

6 

1 

G 

•  *» 

9 

cL 

!' 

9 

10 

£. 

c‘ 

10 

11 

O 

8 

11 

12 

10 

12 

13 

j 

11 

cL 

13 

*7 

12 

13 

6 

3 

13 

■7 

i 

8 

4 

SOIL 

UNIT  HEIGHT 

COHESION 

PHI  SAT'D 

1 

115 

200 

26  1 

•2. 

117 

200 

26  O 

3 

120 

1  00 

28  0 

4 

150 

1  0000 

45  O 

CIRCLE:  X=108.7 


V=1 100.8  RADIUS  =81  FACTOR  OF  SAFETV  =1 


.45417058 
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DO  VOU  WANT  A  DIAGNOSTIC  RUN  <  V  OR  N 
?V 


CE 

WEIGHT 

INCLINATION 

COHESION 

WIDTH 

EFF  WEIGHT 

PHI 

1 

294.477875 

-4.96581439 

200 

2.64388013 

137.423009 

26 

1  O  1 . 66947 •: : 

•p 

2538. 28963 

-2. 01353368 

100 

5.68858236 

1185.89153 

23 

105.855709 

3 

4904.59238 

2.01358368 

100 

5. 68858236 

2290. 17546 

iiO 

111. 54429  L 

4 

7257. 89607 

6.11577104 

200 

5.37047255 

3387. 01816 

26 

117.323819 

5 

9196.21509 

10.3152172 

200 

5.37047255 

4291.56704 

26 

123. 194291 

6 

1O827.2015 

14.5713242 

200 

5.37047255 

5052.69402 

26 

1 29 . 064764 

—f 

1 

11408.9623 

18.9619938 

200 

6 

6425. 16599 

2t- 

135 

3 

9484.07119 

23.3014683 

200 

5. 43929673 

6943. 87233 

26 

140.719648 

8 

8372.38118 

27.5671964 

200 

5.43929673 

7433.64569 

26 

146.  158945 

10 

1556.34339 

30. 1967314 

200 

1 . 12140656 

1556.3 J 339 

26 

149. 439297 

11 

6689. 7093 

32.9750275 

200 

5.5 

6689.7098 

26 

1  r  '  c_- 

1  JiL  .  1  J 

12 

4794.2082 

37.7555341 

200 

5.5 

4794.2082 

26 

158.25 

13 

912.06297 

40.9187278 

200 

1.5 

912. 06297 

26 

161.75 

14 

1399.49512 

44.0393237 

200 

4.91422302 

1399.49512 

26 

1 64 . 957 1 1 2 

NATION  INITIAL 

CALCULATED 

1 

1 

1 . 380869 

36 

2 

1. 

38086936 

1.445166 

03 

•-* 

1. 

44516603 

1 . 453207 

06 

4 

1. 

45320706 

1.454170 

58 

AN  EXAMPLE  PROBLEM 


FACTOR  OF  SAFETY  =  1.45417058  AT 


= 1 93 .7  V  = 1 1 00 . 8  RAD I US  =8 1 


EARTHQUAKE  LOADING  FACTOR  =0 


DO  VOU  WANT  A  FORMAL  PRINTOUT  <Y  OR  N  >'? 
7  N 


DO  VOU  WANT  A  DIAGNOSTIC  RUN  C V  OR  N>? 


139 


DO  YOU  WANT  TO  DEFINE  ANOTHER  FAILURE  CIRCLE  <V  OR  N)? 


FAILURE  CIRCLE  DEFINITION 


X-COORD I NATE  OF  CENTER 

?  1 03 .  7 

V-OOOROINATE  OF  CENTER 
?1 100.3 

CIRCLE  RADIUS 
*40 


FOR  A  RAPID  DRAWDOWN  TYPE  <0>,  FOR  SUBMERGENCE  ELEVATION  AS  PREVIOUSLY  SPECIFIED 
TYPE  < 1 > 

U 


CIRCLE  DOES  NOT  INTERCEPT  SLOPE 

DO  YOU  WANT  TO  DEFINE  ANOTHER  FAILURE  CIRCLE  <Y  OR  N >? 
*N 


BYE-BYE 


/OO  PSF 
120  PoF 


LinE. 


LEFT  PT 


/Ilia Mr  PT 


0>O/ L.  J3E/J T ti 
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/ 

/ 

2. 

3 

2 

2. 

3 

Z. 

3 

3 

4 

1 

4 

A 

3" 

1 

S 

(* 

1 

<£ 

3 

4 

2 

7 

4 

/<3 

2 

© 

o 

20 

/  / 

2 

4 

/  / 

/  2. 

2 

/o 

>2 

>3 

2 

// 

Z. 

/  7 

3 

12. 

1  3 

(p 

J 

1 3 

7 

3 

4 

So/ 1 

C  (P5f) 

& 

S/4TurZAT£Z)  (t-Ms) 

/ 

//4 

zoo 

ztc 

\ 

2 

HI 

zoo 

Z<o 

0 

3 

120 

/  00 

zz 

0 

4  fcac*) 

/  5"0 

J  0,000 

AS 

0 

SA/ou/ZE  CirzciE  CsAjrgfc 

X*  1  cs.i 1  y  =  )\co.  8 '  £ai>ius  -  &) ' 

Porz.  Hul  poo  u  '  fif^noAj  =  /<?4<2 

P7z<?vvi  X- O'  X  =  1 3Z' 


OC S  CYBEfc  730 
0 1  lC£5  LH’ASt  X 

Bishop  j_ 


3P  -  2.7  2, 


$P  =  Z.64 


4  =  Z.4% 


Po/^  PaPiD  'bft  ftUi'fcauJAj  : 


OCS  CYS££  ^30 

TrTices  LEASE  r 

3/5  HOP  i. 


SF  =  /,42- 
6F-  =  /•  4  5” 


4.5.2  Problem  #2 


Example  Problem  Presented  by  Cross  [3] . 


Coha- 

Satu¬ 

Point 

X 

Y 

UlM 

La  ft 

Right 

Soil 

Soil 

a 

cion 

Phi 

rated 

1 

0 

100 

1 

1 

2 

3 

1 

127 

2000 

20 

No 

2 

300 

100 

2 

2 

3 

3 

2 

130 

1000 

33 

No 

3 

400 

150 

3 

3 

4 

2 

3 

130 

1000 

33 

Yes 

4 

500 

200 

4 

4 

5 

1 

5 

800 

250 

5 

5 

6 

1 

6 

1000 

250 

8 

4 

7 

2 

7 

1000 

200 

7 

3 

8 

3 

8  1000  150 


FACTOR  of  SAFETY-  1 .96  AT  X-  400  Y-  300  R-  430 
EARTHQUAKE-  0.03 

DO  YOU  UISH  A  FORMA*.  PRINTOUT  (Y  OR  N > T 
SAMPLE  SLOPE  STABILITY  PROBLEM 
WATER  UNIT  WEIGHT-  42.40  EARTHOUAKE-O. 03 
SUBMERGENCE  AT  130.00  FROM  0.0  TO  400.0 


POINT 

X-ORD 

Y-ORD 

1 

0.00 

100.00 

2 

300 . 00 

100.00 

3 

400.00 

130.00 

4 

500 . 00 

200.00 

5 

600.00 

230.00 

6 

1000.00 

230.00 

7 

1000.00 

200.00 

8 

1000.00 

130.00 

LINE  LEFT 

RIGHT 

SOIL 

1  1 

2 

3 

2  2 

3 

3 

3  3 

4 

2 

4  4 

3 

1 

5  3 

6 

1 

6  4 

7 

2 

7  3 

8 

3 

SOIL 

UNIT  WEIGHT 

COHESION 

PHI 

SATURATED 

1 

127 

2000 

20 

1 

2 

130 

1000 

33 

1 

3 

130 

1000 

33 

0 

CIRCLE  X-ORO  Y-ORD  RADIUS  FACTOR  OF  SAFETY 

400.0  300.0  430.0  1 . 9 A 


DO  YOU  UISH  A  DIAGNOSTIC  RUN  CY  OR  N)T 


SLICE 

WEIGHT 

INCLINATION 

COHESION  WIDTH  EFF  UEIGTH 

1 

252336.4 

-23.6 

1000 

53. 1 

45203.8 

J 

305318.3 

-16.4 

1000 

33. 1 

114232.4 

3 

402330.9 

-9.6 

1000 

30.0 

109976.2 

4 

603232.6 

-3.2 

1000 

30.0 

293401 .0 

3 

726732.6 

3.2 

1000 

30.0 

416901 .0 

6 

052030.9 

9.6 

1000 

30.0 

360476.2 

7 

939262 . 7 

16.2 

1000 

50.0 

602516. 6 

0 

902802.6 

22.9 

1000 

30.0 

701399.0 

9 

759268.3 

27.4 

1000 

41.4 

650306.4 

10 

617173.0 

33.7 

1000 

41.4 

376418.3 

11 

310439.1 

43.6 

1000 

32.6 

310439. 1 

12 

131972.9 

32  •  2 

2000 

38.0 

131972.9 

DRIVING  FORCE  COURIER  BALANCE 

or 

73111.11 

ITERATION 

INITIAL 

CALCULATED 

l 

1 . 0000 

1.8106 

2 

1.8106 

1.9435 

3 

l .9433 

1 . 95  70 

4 

1.9370 

1.9383 

FACTOR  OF  SAFETY-  1.96  AT  X-  400  Y«  500  R-  430 
EARTHQUAKE*  0.03 


Example  Problem  Calculated  by  BISH0P1 


(BISHOP1  is  the  slope  stability  program  presented  by  Cross  [3] 
translated  into  Applesoft) 


SAMPLE  SLOPE  STABILITY  PROBLEM 

***•#  •*•*■*■■*•■*■■*■+■!►•■*■*■  nr  H^**^**^-#*-*- -S' 

FACTOR  OF  SAFETV  =  1 . 35830517  AT  X  =400  V  =590  RADIUS  =450 
EARTHQUAKE  LOADING  FACTOR  =.05 

SAMPLE  SLOPE  ST AS I L I TV  PROBLEM 


WATER  UNIT  HEIGHT  =B2.4  EARTHQUAKE  LOADING  FACTOR 
SUBMERGENCE  ELEU  =150  FROM  K=0  TO  K=400 


POINT 

1 


M-COGROINATE 

Y-CC 

0 

i  00 

300 

100 

400 

150 

500 

200 

S00 

250 

1  000 

250 

i.  ODD 

200 

1  000 

150 

ME  LEFT  F'T  RIGHT  FT 
1  1  2 


UNIT  HEIGHT  COHESION 
127  2000 

130  1000 

1 3ti  1  000 


:«H  i  J 
1 
1 

n 


ClRcLt:  =4P0  V=500  RADIUS  =450  FACTOR  OF  SAFETV  =1 . 05S30517 
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CE 

NtlGHT 

I NCL I NAT I ON 

COHESION 

HIDTH 

EFF  HEIGHT 

PHI 

■;;* 

1 

£52536.452 

-23.5767355 

1000 

53.0776408 

45205 . 7908 

33 

220.36353: 

4 

385318. 321 

-16.363501 

1 000 

53.0776407 

1 14252. 363 

273.46116 

3 

482338 . 898 

-9. 60948742 

1000 

50 

189976. 227 

”7  “7 

4 

603232.626 

-3. 18968777 

1  ooo 

50 

293400 . 965 

33 

375 

5 

726732. 626 

3.  1 8968  f  i''  f' 

1  ooo 

50 

4 1 6900 . 965 

33 

425 

Ei 

852838 . 898 

9. 60943742 

1  ooo 

50 

560476 . 227 

"7  "7 

475 

1 

339262.679 

16. 1554181 

1000 

50 

6825 1 6 . 593 

33 

cr  c: 

,__l  ^  ( 

8 

982882.646 

22.9295296 

1000 

50 

78 1 398 . 976 

33 

575 

9 

759268. 279 

23.4164975 

1000 

41.4213565 

650306 . 432 

33 

620 . 7 1 067 

10 

617175.754 

35.  6938236* 

1000 

41.4213562 

576418.317 

33 

662. 13203 

11 

518439.062 

43.5661001 

1000 

52 . 567 4844 

518439. 062 

~2  "7 

709.  i  26*4.5 

12 

131972. 89 

52.2203924 

2000 

38.7555423 

131972.39 

3*3 

754. 78796 

GIVING  FORCE  COUNTER  BALANCE  OF  751il.il 16#. 


NATION 

1 


INITIAL 

1 

1 . Si 060871 
1 . 94350336 
1 . 9570 1 339 


CALCULATfeO 
1.81 06067 1 
1 . 943“'o336 
1.95701339 
1 . 958305 1. 7 


j AMPLE  SLOPE  STA6ILITV  PROBLEM 

-ACTOR  OF  SAFETV  =  1.95830517  AT  X  =400  V  =500  RADIUS  =450 
EARTHQUAKE  LOADING  FACTOR  =.05 

Results 


Method 
Cross  [3' 
BISHOP  1 


Factor  of  Safety 
1.96 
1.958 
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CHAPTER  V 


FLEXIBLE  PAVEMENT  DESIGN  (AASHTO) 


5 . 1  Problem  Definition 

AASHTO  1  is  a  computer  program  for  the  design  of  flexible 
pavements.  AASHTO  1  is  based  on  the  design  equations  developed 
from  the  AASHO  Road  Test  performed  in  1958  by  the  American 
Association  of  State  Highway  Officials  [1].  (See  Figure  5.1.) 


Layer  # 
1 

2 

3 


4 


a  - 


Asphaltic 

Concrete 


Crushed 

Stone 


Stabilized 

Soil 


Subgrade 


D(i)  -  layer  thickness 


Figure  5.1.  General  Problem  Definition. 


AASHTO  1  is  composed  of  two  basic  programs.  The  first 
program  calculates  the  required  thickness  of  a  layer  given 


the  required  design  parameters  and  the  characteristics  of 
the  other  layers.  The  second  program  calculates  the  number 
of  equivalent  18  kip  wheel  loads  the  pavement  can  endure 
over  a  projected  20-year  service  life.  Both  programs  calculate 
the  structural  number. 

5 . 2  Background  Theory 

5.2.1  General 

The  AASHO  Road  Test  was  performed  in  Illinois  from 
1958  through  1960.  The  test  consisted  of  several  pavement 
types  constructed  in  a  race  track  configuration.  The  pavements 
were  subjected  to  various  wheel  loads  and  durations.  So  much 
data  was  collected  that  it  took  two  years  to  produce  the 
results  and  recommendations  extrapolated  from  the  test. 

The  scope  of  the  test  was  to  produce  a  standardized  design 
based  on  the  useful  service  of  a  pavement  versus  theoretical 
structural  design  criteria.  Rigid  and  flexible  pavements 
were  considered.  This  report  deals  only  with  flexible  pave¬ 
ment  design. 

5.2.2  Flexible  Pavement  Construction 

A  flexible  pavement  consists  of  structural  layers. 

The  capacity  of  a  layer  to  distribute  load  decreases  from 
the  top  layer  down  to  the  subgrade.  In  general,  a  flexible 
pavement  consists  of  a  bituminous  surface  course,  a  base, 
and  a  subbase  which  in  turn  rests  upon  a  soil  subgrade. 


The  surface  course  is  usually  made  of  asphaltic  concrete 
and  is  capable  of  distributing  a  wheel  load  to  the  base 
with  a  minimum  amount  of  distortion  or  consolidation.  By 
distributing  the  constant  pressure  of  a  wheel  load  such 
that  the  pressure  is  reduced  at  some  finite  depth,  the  base 
can  be  made  up  of  a  material  which  is  less  structurally 
capable  than  the  surface  course.  This  rationale  applies 
to  the  subbase  as  well  as  the  subgrade  soil.  From  the 
AASHO  Road  Test,  minimum  thicknesses  for  flexible  pavement 
layers  were  established  as  follows  [1] : 

Surface  Course  2  inches 

Base  Course  4  inches 

Subbase  Course  4  inches 

In  many  cases  a  subbase  course  is  not  used  and  the 

base  course  rests  directly  on  the  subgrade  soil. 

5.2.3  AASHTO  Design  Considerations 

The  principal  factors  of  the  AASHTO  pavement  design 
are  [2]: 

a)  magnitude,  method  of  application,  and  number  of 
wheel  loads 

b)  function  of  pavement  and  base  in  transmitting  the 
load  to  the  subgrade 

c)  measurement  of  the  subgrades  ability  to  support 


the  transmitted  load 


L5 


0 

5. 2. 3.1  Equivalent  18-kip  Loads 

The  AASHTO  design  considers  the  number  of  equivalent 
18-kip  wheel  loads  (E-18's)  applied  to  a  pavement  over  a 
particular  service  life.  Several  methods  of  estimating 
this  quantity  have  been  established  by  various  transportation 
and  highway  organizations.  Methods  range  from  traffic  counters 
with  load  meters  to  extrapolating  historic  data.  In  many 
cases,  organizations  will  expend  more  funds  on  collecting 
this  data  than  can  be  rationalized  for  its  intended  purpose. 
Historic  data  will  usually  suffice  as  a  means  of  projecting 
anticipated  wheel  loads.  This  is  especially  true  when  new 
interstates  are  being  constructed  as  other  factors  will 
influence  usage  such  as  route  direction,  load  limitations, 
city  connections,  and  other  regional  limitations.  Large 
amounts  of  data  exist  for  equating  given  loads  and  load 
configurations  to  18-kip  single  axle  loads.  Service  organiza¬ 
tions  must  adapt  a  method  of  calculating  E-18's  best  suited 
for  their  region. 

5. 2. 3. 2  Soil  Support 

In  the  AASHO  Road  Test,  a  soil  support  value  of  3  was 
established  to  represent  the  subgrade  soil  used  which  was 
an  A-6  soil  [2] .  Crushed  stone  with  a  California  Bearing 
Ratio  (CBR)  of  200  was  assigned  a  soil  support  value  of  10. 

The  soil  support  value  of  various  soils  have  been  estab¬ 


lished  by  different  service  organizations.  Figure  5.2  and 
Tables  5.1  through  5.4  represent  a  selected  group  of  soil 


Table  5.1.  Maximum  Recommended  Soil  Support  Values  [6] 


Classification 

Description 

Upper  Soil 
Support  Value 

A- la 

Largely  gravel  but  can 
include  sand  and  fines 

6.5 

A- lb 

Gravelly  sand  or  graded 
sand;  may  include  fines 

6 

A-2-4 

Sands,  gravels  with  low 
plasticity  silt  fines 

5 

2-2-4 

Micaceous  silty  sands 

2.5  -  3.0 

A- 2- 5 

Sands,  gravels  with 
plastic  silt  fines 

4 

A- 2-6 

Sands,  gravels  with 
clay  fines 

4.0  -  5.0 

A- 2- 7 

Sands,  gravels  with 
highly  plastic  clay 
fines 

4.0 

A- 3 

Fine  sands 

4.5 

A- 4 

Low  compressibility 
silts 

4.0 

A- 5 

High  compressibility 
silts,  micaceous  silts 
and  micaceous  sandy 
silts 

2.5  -  3.5 

A*  6 

Low  to  medium  compress¬ 
ibility  clays 

3.5  -  4.5 

A- 7 

p 

High  compressibility 
ciays,  silty  clays  and 
high  volume  change 
clays 

3-4 
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Table  5.2.  Soil  Support  Values  [2]. 


~r 

7 


~T 

8 


Soil  Support  Value 


support  values  based  on  widely  varying  soil  characteristics. 
Each  organization  must  adopt  a  means  of  evaluating  the  soil 
support  values  which  best  suit  the  soils  encountered  in 
a  particular  region. 

5.2. 3.3  Terminal  Serviceability  Index 

A  significant  concept  developed  by  the  Road  Test  was 
the  Terminal  Serviceability  Index  (P  ) •  The  P  is  a  quali¬ 
tative  measure  of  the  final  condition  of  a  pavement  at  the 
end  of  its  design  service  life.  The  P  scale  is  from  total 
falure  at  1.5  to  outstanding  at  5.0;  3.5  to  5.0  corresponds 
to  poor  to  outstanding  new  construction,  respectively.  For 
design,  a  Pfc  of  2.5  applies  to  the  minimum  serviceability 
of  an  interstate  highway;  similarly,  a  Pfc  of  2.0  applies 
to  secondary  roads. 

5.2. 3.4  Regional  Factors 

In  an  effort  to  apply  the  general  design  equations 
developed  via  the  Road  Test  to  other  regions,  a  ret  onal 
factor  was  introduced  such  that  environmental  factors  not 
encountered  in  Illinois  could  be  incorporated  into  the  design 
The  regional  factor  can  incorporate  the  following  parameters 
[1]  : 

a)  topography 

b)  similarity  to  the  Road  T  st  location 

c)  rainfall 

d)  frost  penetration 

e)  temperature 


f)  groundwater  table 

g)  subgrade  type 

h)  engineering  judgment 

i)  type  of  highway  structure 

j)  subsurface  drainage 

Figure  5.3  is  a  general  guide  for  the  selection  of 
a  regional  factor.  The  scale  varies  from  0.5  to  4.8.  The 
lowest  values  apply  to  permanently  frozen  or  consistently 
dry  roadbed  materials.  The  upper  values  apply  to  severe 
frost  heave  conditions  and  other  mechanisms  which  rapidly 
accelerate  pavement  deterioration.  The  value  used  for  the 
road  test  was  1.0.  AASHTO  recommends  the  following  as 
a  crude  guide. 


Table  5.5.  Regional  Factors  [1], 


Condition 

R  value 

Roadbed  materials  frozen  to  depth  of  5  in.  or  more 
Roadbed  materials  drv,  summer  and  fall 

Roadbed  materials  wet.  spring  thaw 

0.2-1 .0 
0.3-1  .5 

4. 0-5.0 

Historical  data  of  pavement  performance  in  relation  to  the 
number  of  annual  freeze-thaw  cycles,  steep  grades  with  large 
volumes  of  heavy  truck  traffic,  and  areas  of  concentrated 
turning  and  stopping  movements  can  be  useful  in  evaluating 


an  appropriate  regional  factor. 


5.2.3. 5  Structural  Number 

The  E-18s,  regional  factor,  soil  support  value,  and 
terminal  serviceability  index  represent  various  components 
influencing  the  design  of  a  flexible  pavement.  To  combine 
the  above  parameters  into  a  value  which  could  be  translated 
into  a  pavement  design,  a  structural  number  (SN)  was  intro¬ 
duced.  The  structural  number  is  a  value  calculated  from 
the  above  parameters  and  is  the  sum  of  the  layer  thicknesses 
multiplied  by  their  appropriate  structural  coefficients. 


SN  —  a^D^  +  ^2^2 


•  +  aiDi 


where , 

a^  =  structural  coefficient  of  layer  i 
=  thickness  of  layer  i  (inches) 

The  structural  number  is  abstract  but  relates  to  the  strength 
of  the  section. 

5.2.3. 6  Structural  Coefficient 


The  structural  coefficient  is  a  measure  of  a  layers  ability 
to  transmit  load.  As  previously  discussed,  the  surface 
course  will  have  a  higher  structural  coefficient  than  a 
layer  less  capable  of  distributing  a  load  to  a  lower  layer. 
Tables  5.6  through  5.9  relate  structural  coefficients  to 
various  construction  materials.  Many  service  organizations 
develop  their  own  correlations  to  best  suit  their  situation. 


Table  5.6.  Structural  Coefficients  [1] • 


Structural  Layer  Coefficients  Proposed  by  AASHO  Committee  on  Design, 
October  12,  1961 


Pavement  Component 

Coefficient1 

Surface  Count 

Road  mix  (low  stability) 

0.20 

Plantrrax  (high  stability) 

0.44* 

Sand  Asphalt 

0.40 

Base  Count 

Sandy  Gravel 

0.071 

Crushed  Stone 

Cement-Treated  (no  soil-cement) 

Compressive  strength  @  7  days 

650  psi  or  more1  (4.48MPa) 

0.14* 

0.23* 

400  to  650  psi  (2.76  to  4.48MPa) 

0.20 

400  psi  or  less  (2.76MPa) 

BituminouvTreated 

0.15 

Coane-G  faded 

0.34' 

Sand  Asphalt 

0.30 

Lime-Treated 

0.15-0.30 

Subbase  Course 

Sandy  Gravel 

0.11* 

Sand  or  Sandy-Clay 

0.05-0.10 

*  Established  from  AASHO  Road  Test  Data 
1  Compressive  strength  at  7  days. 

1  This  value  has  been  estimated  from  AASHO  Road  Test  data,  but  not  to  the  accuracy  of 
those  factors  marked  with  an  asterisk. 

1  It  is  expected  that  each  state  will  study  these  coefficients  and  make  such  changes  as 
experience  indicates  necessary. 
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Table  5.7.  Selected  Structural  Coefficients  Used  by  Various 
Transportation  Organizations  in  the  AASHO  Interim 
Guide  Design  Method  [6]  . 


Pavement  Component 

Structural 

Coefficient 

Fla. 

Ga. 

Md. 

S.C. 

I.  Surface  and  Binder  Course  (a,) 

■M 

1 

(1) 

Asphalt  Concrete 

fir  BP 

0.44 

0.44 

0.44 

Bituminous  Surfacing 

■am 

' 

0.35 

II.  Base  Course  (a.) 

(4) 

Asphalt  Concrete 

0.21-0.30 

0.30 

0.28 

0.34 

Sand-Asphalt 

0.12 

0.28 

0.20-0.25 

Soil  Cement 

0.22 

0.20 

0.28 

0.20 

Graded  Aggregate 

0.18<1 2) 3 4 5 

0.14 

0.12-0.20 

Cement  Stabilized  Graded 

0.22 

0.28 

0.34 

Aggregate 

Sand  Aggregate 

- 

- 

0.14 

Sand-Clay  CBR  >  49 

0.12 

- 

- 

Limerock  CBR  >  80 

0.15 

- 

- 

- 

Limerock  Stabilized 

0.12 

_ 

Base.  CBR  >  56 

III.  Subbase  (a^) 

Graded  Aggregate 

0.14 

- 

- 

Topsoil  or  Sand-Clay 

0.10 

- 

- 

Gravel  or  Screenings 

- 

- 

0.07 

- 

Soil  Aggregate 

- 

- 

- 

0.08  -  0.12 

Cement  Stabilized  Earth 

0.15 

1.  Georgia  uses  a  coefficient  of  0.44  for  surface  and  binder  to  a  depth  of  4.5  in. 
(110  mm) . 

2.  When  compacted  to  100Z  of  T-180  density. 

3.  Subbase  coefficients  are  used  in  Georgia  below  a  depth  of  12  in.  (300  mm). 

4.  The  Florida  DOT  uses  the  following  structural  coefficients  for  different  base 

mixes:  (1)  Type  I,  500  lb.  Marshall  stability,  a2“0.21;  (2)  Type  II,  750  lb. 

Marshall  stability,  a2*0.25;  (3)  Type  III,  1,000  lb.  Marshall  stability, 

a,  =  0.30. 

5.  The  Florida  DOT  uses  the  following  structural  coefficients  for  different  surface 
mixes:  Type  SI,  1,000  lb.  Marshall  stability,  aj^-0.40;  Type  S2,  1,000  lb. 
Marshall  stability,  a2”0.20;  and  Type  S3,  750  lb.  Marshall  stability,  a^*0.30. 
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Table 

5.8.  Pavement  Coefficients  for  Flexible  Section  Design,  Louisiana 

[1] 

Strength 

Coefficient 

i. 

SURFACE  COURSE 

Asphaltic  Concrete 

- 

Types  1.  2  and  4  BC  and  WC 

1000* 

0.40 

Types  3  WC 

1800+ 

0.44 

BC 

1500+ 

0.43 

n. 

EASE  COURSE 

UNTREATED 

Sand  Clay  Gravel  -  Grade  A 

3.3- 

0.08 

Sand  Clay  Gravel  -  Grade  B 

3.5- 

0.07 

Shell  and  Sand  -  Shell 

2.2- 

0.10 

CEMENT-TREATED 

Soil-Cement 

300  psH- 

0.15 

Sand  Clay  Gravel  -  Grade  B 

500  pst+ 

0.18 

Shell  and  Sand  -  Shell 

500  psr+ 

0.18 

Shell  and  Sand  -  Shell 

650  psr+ 

0.23 

LIME-TREATED 

Sand  Shell 

2.0- 

0.12 

Sand  Clay  Gravel  -  Grade  B 

2.0- 

0.12 

ASPHALT-TREATED 

Hot-Mix  Base  Course  (Type  5A) 

1200+ 

0.34 

Hot-Mix  Base  Course  (Type  SB) 

800+ 

0.30 

HI. 

SUBBASE  COURSE 

Lime-Treated  Sand  Clay  Gravel  -  Grade  B 

2.0- 

0.14 

Shell  and  Sand-Shell 

2.0- 

0.14 

Sand  Clay  Gravel  -  Grade  B 

3.5- 

0.11 

Lime-Treated  Soil 

3.5- 

0.11 

Old  Gravel  or  Shell  Roadbed  (8"  thickness)  (200  mm) 

- 

0.11 

Sand  (R-Value) 

55+ 

0.11 

Suitable  Material-A-6  (PI  *  15-) 

- 

0.04 

Table  5.9.  Recommended  AAS1IT0  Interim  Guide  Structural  Coefficients 
for  Thickness  Design  [6], 


Uns t  abi 1 i zed  S  u  t  id  - 
crushed  scene  blend 


-V  Soil  support  value 


is  20.000- 

s  3  lo.ooo- 

~  2  5000- 


3000  i  9 
2000  m  2 


Design  chart  for  flexible  payments,  p,  =  2.0 


i! 

*  o 

CL  -5 

20.000 

10.000 

5000 

2  | 

1000 

S  y 

■o  o. 

500 

5 

100 

3  ■a 
—  o 

50 

3000  » 

2000  «  ^ 
1000  5  „  =  E 
500  fj  £  I 
'  ^ 

?ll  I 


■2  -0  5 
"a  -  1.0 
9-20 
J?  -50 


Design  chart  for  flexible  pavements,  p,  =  2.5 

Figure  5.4.  AASHO  Flexible  Pavement  Design  Nomographs  [3] 
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5.2.4  AASHTO  Equation 

The  general  equation  developed  by  AASHTO  is  as  follows: 

log  [  (4.2-P.  )/2.7] 
log  Wlg  =  9.36  log  (SN+1)  -  .2  +  - - - - 

*  ^  +  5  Tgi 

(SN+1 )  3  *  ■Ly 


+  log  4  +  -372  (S.  -  3.0) 

K.  1 


where , 

W^g  =  E-18's,  single  axle  loads 
SN  =  structural  number 

=  terminal  serviceability  index 
R  =  regional  factor 
S^  =  soil  support  value 
(log  =  natural  logarithm) 

The  structural  number  is  calculated  by  iterative  trials 
given  W^g ,  Pfc,  R,  and  S^. 

Nomographs  were  developed  for  terminal  serviceability 

indexes  of  2.0  and  2.5  and  are  presented  in  Figure  5.4. 

After  calculating  SN  with  S-  and  W. Q ,  a  weighted  or  design 

1  lo 

SN  is  calculated  using  the  regional  factor. 

Figure  5.5  relates  the  variation  of  wheel  load  type 
to  the  structural  number  with  all  other  parameters 
normalized  to  the  road  test;  i.e.,  S^  =  3.0  and  R  =  1.0. 


(Structural  Number) 


Figure  5.5.  AASHO  road  test  relationship  between  thickness  index  and  axle  loads 
:.t  /.  -  2.5.  [2]  . 

5 . 3  Program  Rationale 

AASHTO  1  is  designed  to  provide  the  user  the  ability 
to  approach  the  design  problem  by  two  avenues.  Given  the 
regional  factor,  soil  support  value,  and  terminal  service¬ 
ability  index,  the  user  can  input  the  pavement  layer  charac¬ 
teristics  and  calculate  the  number  of  equivalent  single 
axle  18-kip  loads  the  pavement  can  expect  to  endure.  Similar¬ 
ly,  the  program  can  calculate  the  required  thickness  of 
a  layer  given  the  regional  factor,  soil  support  value,  terminal 
serviceability  index,  E-18's,  and  the  characteristics  of 
the  other  layers  (see  Figure  5.6). 


Input  —  _ 

GOSUB 

T 

Tables 

Thickness  • _ . 

/  Calculated 
\  Thickness  \— 

1 

\or  E-18  1  s  / 

- E-18  *  s 

Calculate 


Input 

Layer 

Parameters 


Tables 


Input 

Layer 

Parameters 


Tables 


Calculate 


Calculate 

Layer 

Thickness 


Calculate 

E-18's 


Print 


Print 


Figure  5.6.  aafht n  i 

AAbHTO  l  Flow  Chart 
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Originally,  the  author  intended  to  create  a  matrix 
of  possible  layer  thicknesses  but  due  to  the  limited  capacity 
of  the  miniature  computer,  this  idea  was  aborted.  The  pro¬ 


gram  was  then  designed  to  allow  the  user  to  quickly  change 
thicknesses  and  structural  coefficients  in  successive  runs 
without  inputting  repetitive  data. 

The  output  was  designed  to  resemble  a  layered  system 
with  input  variables  and  other  calculated  data  presented 
below  the  output.  This  format  allows  the  user  to  correlate 
data  from  different  runs  with  the  maximum  of  ease. 

5 . 4  Program  Use  and  Limitations 

5.4.1  General 

AASHTO  1  is  user  oriented.  Input  is  prompted  by  state¬ 
ments  and  questions  which  clearly  identify  the  input  require¬ 
ments.  As  with  any  program,  the  user  should  have  the  input 
ready  prior  to  running  the  program.  If  several  runs  are 
to  be  executed,  the  user  should  prepare  a  chart  of  required 
input  variables.  This  method  will  eliminate  most  errors 
associated  with  mistaken  or  improper  variable  input. 

5.4.2  Input 

As  previously  mentioned,  AASHTO  1  is  composed  of  two 
parts.  One  portion  of  the  program  calculates  a  required 
thickness  given  all  other  data,  and  the  second  portion  cal¬ 
culates  the  equivalent  single  axle  18-kip  loadings.  The 
user  will  be  prompted  to  choose  which  option  is  to  be  run. 

The  user  must  input  the  number  of  runs  to  be  made  with  the 
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portion  of  the  program  previously  specified.  After  the 
runs  are  completed,  the  user  may  then  opt  to  run  the  alter¬ 
nate  portion  of  the  program  or  more  runs  of  the  same.  Each 
time  a  new  set  of  runs  are  specified,  all  input  must  be 
re-entered . 

The  following  is  a  list  of  required  input.  The  input 
differences  between  the  two  programs  will  be  discussed  later 
in  this  section. 

a)  soil  support  value 

b)  terminal  serviceability  index 

c)  regional  factor 

d)  E-18's 

e)  number  of  layers 

f)  layer  number  for  thickness  calculation 

g)  layer  number 

h)  structural  coefficient 

i)  layer  thickness 

When  the  user  has  specified  an  E-18  calculation,  input 
d)  and  f)  will  be  emitted;  conversely,  for  a  layer  thickness 
calculation  input  i)  will  be  emitted  for  the  layer  specified 
for  thickness  calculation.  Although  it  is  recommended  to 
number  the  layers  successively  from  #1  for  the  top  layer 
to  #i  for  the  ith  layer,  this  is  not  required.  The  user 
will  find  the  output  format  suited  to  numbering  layers  in 


this  fashion. 


5.4.3  Options 

AASHTO  1  provides  the  user  with  tables  of  typical  values 
of  soil  support,  tej.  Tiinal  serviceability  index,  regional 
factor,  E-18's,  and  structural  coefficients.  The  user  may 
access  these  charts  by  inputting  0  (zero)  instead  of  the 
non-zero  variable  requested.  After  inputting  0,  the  chart 
will  appear  on  the  monitor.  When  the  user  is  finished  re¬ 
viewing  the  chart,  the  user  will  again  be  prompted  to  input 
the  required  variable.  The  charts  are  derived  from  those 
presented  in  Section.  5.2. 

As  previously  described,  the  user  has  the  option  of 
calculating  a  layer  thickness  or  the  number  of  equivalent 
single-axle  18-kip  loads. 

5.4.4  Limitations 

The  user  is  limited  to  specifying  ten  layers.  This 
limitation  is  due  to  memory  capacity.  If  the  user  requires 
more  layers,  there  are  no  constraints  in  the  program  to 
prevent  inputting  more  than  ten.  The  user  should  not  attempt 
to  input  more  than  ten  layers.  If  an  out-of-memory  error 
results,  the  user  must  reduce  the  number  of  layers. 

The  inherent  limitation  of  AASHTO  1  is  the  basis  and 
validity  of  the  AASHO  Road  Test  equations.  The  user  is 


urged  to  study  reference  #4  for  an  evaluation  of  the  guile. 
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iram  List 


HGnE 
SPEED® 
PR  I  NT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


*  AASHTO 


:  PRINT  :  PRINT 
"DANA  K.  EDDY,  573-30-3378 
"60.  INSTITUTE  OF  TECHNOLOGV" 

"SCHOOL  OF  Cl OIL  ENGINEERING" 

"DEPORTMENT  OF  GEOTECHNICAL  ENGINEERING" 

:  PRINT  :  PRINT 

"SYSTEM  HARDWARE:  APPLE  II  PLUS  <64K>" 

"SYSTEM  HARDWARE:  DOS  3.3,  APPLESOFT  BASIC  LANGUAGE" 
"PROGRAM  DATE:  JULY,  1983” 

:  PRINT 
:  PRINT 

PRINT  "AASHTO  IS  BASED  ON  THE  A ASH TO  INTERIM  GUIDE  FOR  THE  DESIGN 
LEGIBLE  PAUEMENT  STRUCTURES,  1972." 


+4 

SPEED®  2! 

46 

HOME 

l  00 

PRINT  "1 

E  SUPPLI 

(O)  IN 

110 

PRINT  : 

1  15 

INPUT  " 

117 

PRINT  : 

i  20 

INPUT  " 

PRINT  :  PRINT  :  PRINT  "< PRESS  THE  SPACE  BAR  TO  CONTINUE >" 
POKE  -  16368,0 
CH  =  PEEK  <  -  16384 > 

IF  CH  <  >  160  GOTO  40 


PRINT  "USER  INPUT  QUESTIONS  AND  COMMANDS  WHICH  ARE  FOLLOWED  BY  <*>  AR 
E  SUPPLEMENTED  BY  LISTS  AND  TABLES.  TO  ACCESS  THE  LIST  OR  TABLE  TYPE 
<  O  >  INSTEAD  OF  THE  APPROPRIATE  UALUE  OF  THE  UARIABLE  REQUESTED." 
PRINT  :  PRINT 

INPUT  "PROBLEM  HEADING  ";H*  ' 

PRINT  :  PRINT 

INPUT  "TO  CALCULATE  THE  THICKNESS  OF  A  PARTICULAR  LAYER,  TYPE  < 1  . > ;  TO 
CALCULATE  THE  MAXIMUM  NUMBER  OF  EQUIVALENT  18-KIP  LOADS,  TYPE  <  2  " 

iBB 

PRINT  :  PRINT 

INPUT  "HOW  MANY  TRIALS  DO  YOU  WANT  TO  RIJN?  "  ;Z 
PRINT  :  PRINT 
IF  BB  =  2  THEN  ISO 

INPUT  "WHAT  IS  THE  CALCULATION  TOLERANCE  <REC'0;  .01).  ";TL 
PRINT  :  PRINT 

INPUT  "SOIL  SUPPORT  UALUE  <*>  =  ";SI 

PRINT  :  PRINT 

IF  SI  >  O  GOTO  240 

GOSUB  200O 

GOTO  190 

INPUT  "TERMINAL  SERICEABILITY  INDEX  <  *  )  =  ">PT 

PRINT  :  PRINT 

IF  PT  >  O  GOTO  290 

GOSUB  300O 

GOTO  240 

INPUT  "REGIONAL  FACTOR  <*)  =  ";RI 
PRINT  :  PRINT 
IF  R I  >  O  GOTO  340 
GOSUB  4000 
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330  GOTO  230 

340  IF  BB  =  2  GOTO  1010 

350  INPUT  "#  OF  EQUIVALENT  SINGLE  PXLE  13-KIP  WHEEL  LORDS  <*>  =  ";WI 

352  PRINT  :  PRINT 

355  IF  WI  >  0  GOTO  370 

360  GOSUB  5000 

365*  GOTO  350 

370  X  =  4 

375  XX  =  1 : VV  =  0:Z2  =  0 
330  PI  =  LOG  <  HI )  /  2.3026 
330  P2  =  .2 

400  P3  =  LOG  <<4.2  -  PT )  ✓  2.7)  /  2.3026 
410  P4  =  LOG  <1  /  RI)  /  2.3026 
420  h5  =  .375  *  < SI  —  3.0) 

430  HP  =  fil  +  P2  -  P4  -  P5 

440  PB  =  <4.065  *  '  LOG  <X>  *  <.4  *  X  ^  5.19  +  1034)  +  P3  *  X  5.13  > 
*  X  5.19  +  1094) 

442  XX  *  XX  +  1 

444  IF  XX  >  500  GOTO  510 

450  IF  PB  >  =  PP  -  TL  /  2  PND  PB  <  =  PP  +  TL  2  GOTO  51 0 

46G  IF  PB  <  PP  GOTO  490 

470  IF  VV  =  0  GOTO  474 

471  IF  VV  =  1  PND  ZZ  =  1  GOTO  474 

472  TL  =  TL  ✓  10 
474  X  =  X  -  TL 
476  ZZ  -  1 

430  GOTO  440 

490  IF  ZZ  =  0  GOTO  494 

491  IF  ZZ  =  1  PND  VV  =  1  GOTO  494 

492  TL  =  TL  /  10 
494  X  =  X  +  TL 
496  VV  =  1 

50O  GOTO  440 
510  SN  =  X  -  1 

512  NN  =  SN  *  100 

513  NN  =  I NT  <  NN  ■ 

514  NN  =  NN  100 


520 

FOR  H 

=  1  TO  Z 

521 

HOME 

524 

PRINT 

1  *  V  V  V  V  V  V  V  H 

^  I 

c  -.=■ 

._j  21 

PRINT 

"TRIPL  #"H 

526 

PRINT 

ii  v V  V  11 

527 

PRINT 

:  PRINT 

er 

PRINT 

"STRUCTURPL  NUMBER  ="NN 

529 

PRINT 

:  PRINT 

530 

INPUT 

" SPEC I FV  #  OF  LPVERS  <MPX-10).  " ;N 

540 

PRINT 

:  PRINT 

550 

INPUT 

"LP’r'ER  #  FOR  THICKNESS  CPLCULPTION. 

560 

PRINT 

:  PRINT 

570 

FOR  I 

=  1  TO  N 

575 

PRINT 

11  ****"*-t-t-1r*-fr*  " 
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577  PRINT  - - - - 

580  INPUT  "LAYER  #.  ";Y 
530  PRINT  :  PRINT 

800  INPUT  "STRUCTURAL  COEFFICIENT  <*>  =  ";ACV> 

610  PRINT  :  PRINT 

620'  IF  AC  V  )  >  0  GOTO  645 

630  GOSUB  6000 

640  GOTO  600 

645  IF  I  =  L  GOTO  700 

650  INPUT  "LAVER  THICKNESS  < INCHES)  =  ";DCV) 

660  PRINT  :  PRINT 

670  IF  OCY)  >  0  GOTO  700 

630  GOSUB  7000 

630  GOTO  650 

700  NEXT  I 

710  SC  =  0 

720  FOR  I  =  1  TO  N 

730  IF  I  =  L  GOTO  750 

740  SC  =  SC  +  DC  I )  *  AC  I > 

750  NEXT  I 

760  IF  SC  <  SN  GOTO  330 

770  PRINT  "**  WARNING  **  INPUT  LAVERS  SATISFY  MINIMUM  STRUCTURAL  REQUIRE 
MENTS. " 

730  PRINT  :  PRINT 

730  INPUT  "TYPE  Cl)  TO  REENTER  UARIABLES,  TYPE  <2)  TO  END  PROGRAM.  " ;OD 

300  PRINT  :  PRINT 

310  IF  DD  =  1  GOTO  520 

320  GOTO  1490  - 

330  DC  L  )  =  SN  /  AC  L  ) 

340  FOR  I  =  1  TO  N 

350  IF  I  =  L  GOTO  370 

360  DC  L  >  =  DC  L  )  -  AC  I  )  *  DC  I  )  /  AC  L  > 

370  NEXT  I 

330  DC  L  )  =  DC  L  )  *  100 

390  DCL)  =  I  NT  CDCL)) 


300 

302 

305 

90S 

DCL)  =  DCL)  '  100 

Li  =  CHRi  C  4 ) :  PRINT  Li; "PR# 

PRINT  :  PRINT  :  PRINT 

PR  T  NT  "  s. s'AV'Z- V V 11 

907 

PRINT 

Hi 

903 

PRINT 

ii  «.•*/  W  W  W V ■-* v  w  v W  V  w •  '  v •••  it 

309 

PRINT 

:  PRINT 

910 

PRINT 

"TRIAL  #"H 

915 

PRINT 

"******•»*" 

920 

PRINT 

:  PRINT 

330 

PRINT 

"  LAVER  "  ,  "  STR  L  COEFF .  "  ,  "  TH I  CM 

950 

FOR  I 

=  1  TO  N 

5  Ay 

PRINT 

SPCC  2  1 1  •  AC  I  )  *  DC  I  > 

97  O 

NEXT 

I 

330 

PRINT 

990 

PRINT 

"STRUCTURAL  HUMBER  ="SN 

932  PR  If 
395  PR  It 

997  PR 1 1 

998  PR  1 1 
1088  ne; 

1005  PR 
1006'  GO 
1010  FO 

1015  HOI 

1016  PR 

1017  PR 

1018  PR 
1013  PR 
1020  IN 
1030  PR 
1040  FG 
1045  PR 
1047  PR 
1050  IN 
I860  PR 
1070  IN 
1080  PR 
1030  IF 
1 1 00  GC 
1110  GC 

1120  ri¬ 
ll  30  PF 
1140  IF 
1150  GC 

1 1 60  GC 

1161  NE 
1163  SN 
1165  FC 
1167  SN 

1169  Nt 

1 1 70  B2 
i 1 80  B3 
I  1 30  B4 
1200  B5 
1210  86 
1220  B7 
1230  B1 
1240  HI 
1250  HI 


PRINT 

PRINT  "EQUIUOLENT  SINGLE  OXLE  1S-KIP  LORO 
PRINT  :  PRINT  :  PRINT 
PRINT  L$i"PR#0" 

NEXT  H 

PRINT  L$;"PR#0" 

GOTO  1370 

FOR  H  =  1  TO  2 

HOME 

print  ny:/x/.y.v:/xn 

PRINT  "TRIOL  #"H 
PRINT  "J00£JS*JS5J5S" 

PRINT  :  PRINT 

INPUT  "SPECIFV  #  OF  LOVERS  CMOX-10)  ";N 
PRINT  :  PRINT 
FOR  I  =  1  TO  N 

PRINT  "-if-*******************" 

PRINT 

INPUT  "LOVER  #  =" ;V 
PRINT  :  PRINT 

INPUT  “STRUCTUROL  COEFFICIENT  <*>  =";flc' 

PRINT  s  PRINT 

IF  fl<  V  >  >  0  GOTO  1120 

GOSUB  6*000 

GOTO  1070 

INPUT  "LOVER  THICKNESS  <  INCHES  )  =  ”;CKV 
PRINT  :  PRINT 
IF  OCV>  >  0  GOTO  1161 
GOSUB  7000 
GOTO  1120 
NEXT  I 
SN  =  0 

FOR  I  =  1  TO  N 
’  SN  =  SN  +  OC  I  >  *  DC  I  > 

NEXT  I 


=  ";V 


PRINT 

INPUT 

PRINT 


DC  I  > 


LOG 

.4  + 


<  C  4 . 1 

1034 


-  PT  > 

-  C  SN  1 


=  LOG  C 1 


=  B2  -  B3 


RI>  /  2.3026 
SI  -  3.0) 

+  B4  /  B5  +  B6 


C  2. 3026 
CHI 


1260  LS  =  CHR$  C  4  > : 

1261  PRINT  :  PRINT 

1262  PRINT 

1263  PRINT  H* 

1264  PRINT  "XXX’.XXr 

1265  PRINT  :  PRINT 
1270  PRINT  "TRIOL  i 


PRINT  L$;"PR#1’’ 
:  PRINT 


,-t •■'»  •''»  ■'« /* 'i <*■  '« <*'»  •’ a * ’•  ••  •  * a •" •  Sm •* •  *  •  *’ a 

H* 

:  PRINT 
"TRIOL  #"H 
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1272  PRINT  "********" 

1280  PRINT 

1285  PRINT  "LOVER" , "STR'L  COEFF. ", "THICKNESS" 

1290  FOR  I  =  1  TO  N 

1295  PRINT  SPC<  3  >1  ,A<  I  >,D<  I  > 

13O0  NEXT  I 

1310  PRINT  :  PRINT 

1320  PRINT  "STRUCTURAL  NUMBER  =”SN 

1330  PRINT 

1340  PRINT  "EGUIUALENT  SINGLE  AXLE  13-KIP  LOADS  ="WI 
1350  PRINT  L$;"PR#0" 

1360  NEXT  H 

1370  Li  =  CHRi  < 4 > :  PRINT  Li;"PR#l" 

1380  PRINT  :  PRINT 

1390  PRINT  "SOIL  SUPPORT  UALUE  ="SI 

140G  PRINT 

1410  PRINT  "TERMINAL  SERUICEABILITV  INDEX  ="PT 
1420  PRINT 

1430  PRINT  "REGIONAL  FACTOR  ="RI 
1435  PRINT  L$;"PR#0" 

1440  INPUT  "DO  VOU  WANT  TO  RUN  ANV  MORE  TRIALS  <0=NO,  1=VES>?  ";Q1 

1450  PRINT  :  PRINT 

1460  IF  Q1  =  O  GOTO  1490 

1470  INPUT  "FOR  A  LAVER  THICKNESS  CALCULATION,  TVPE(  1>,‘  FOR  AN  18-K 
DING  CALCULATION,  TVPE<2>.  ";BB 


i475 

PRINT 

:  PRINT 

1480 

GOTO  140 

1 490 

HOME 

1435 

PRINT 

"THANK  VOU  FOR  USING  AASHTO" 

1500 

PRINT 

1510 

PRINT 

"BVE-BVE" 

1520 

END 

2000 

HOME 

2101 

PRINT 

TAB<  12 >" SOIL  SUPPORT  UALUE" 

2102 

PRINT 

TAB<  17  >"<  AASHO 

2103 

PRINT 

TAB<  12  >"******************" 

2104 

PRINT 

2 1 05 

PRINT 

"  SOIL  DVNAMIC  STATIC  AASHO  MODULUS" 

2106 

PRINT 

"  SUPPORT  CBR  CBR  3  PT 

MR" 

2 1 07 

PRINT 

"  UALUE 

PS  I 

il 

2 1 08 

PRINT 

2109 

PRINT 

TAB<  5  >"8 - 60  - 73 

14.5  — 

NA" 

2110 

PRINT 

TAB<  5  >"7  -  35  -  36  - 

9.75  — 

NA" 

21 1 1 

PRINT 

TAB(  5  >"6 - 17  - 19  - 

6.75  — 

9300' 

2112 

PRINT 

TABC  5  >"5 -  8  - 11  - 

4.5  — 

6400 1 

2113 

PRINT 

TAB<  5  > "  4 -  5  -  6  - 

2.5  — 

4400 

2114 

PRINT 

TAB<  5  ,M3 -  3  - 2.8  - 

1.25  — 

3000 

2115 

PRINT 

TAB*  5  > "  2 - 1.5 1.3  - 

0.  5u  — 

2  1  00 

2116 

PRINT 

TAB-;  5 > "  1 - 0.5 - 0.5  - 

0.25  — 

NA" 

2120 

PRINT 

:  PRINT  :  PRINT 

2122 

PRINT 

:  PRINT 

2130 

IF  Cl 

=  1  GOTO  2300 
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2200 

PRINT  " 

2210 

POKE  - 

2215 

IF  CJ  < 

2300 

HOME 

£301 

PRINT 

2302  * 

PRINT 

2303 

PRINT 

2304 

PRINT 

2305 

PRINT 

2306 

PRINT 

2307 

PRINT 

2308 

PRINT 

2309 

PRINT 

2310 

PRINT 

2311 

PRINT 

2320 

PRINT  : 

:  PRINT 

£322 

IF  Cl  = 

2325 

PRINT  • 

2330 

POKE  ■ 

2335 

IF  CK  < 

2900 

IF  Cl  = 

2920 

HOME 

2925 

INPUT 

2930 

IF  Cl  = 

2935 

L$  =  Cl 

2940 

GOTO  2 

2945 

PRINT  1 

2950 

Cl  =  0 

2955 

RETURN 

3000 

HOME 

3101 

PRINT 

3102 

PRINT 

3103 

PRINT 

3104 

PRINT 

3105 

PRINT 

3 1 06 

PRINT 

3107 

PRINT 

3 1 0i8 

PRINT 

3109 

PRINT 

3110 

PRINT 

3111 

PRINT 

3112 

PRINT 

31 13 

PRINT 

3114 

PRINT 

3115 

PRINT 

3116 

PRINT 

31 17 

PRINT 

3118 

PRINT 

3200 

PRINT 

3205 

PRINT 

SOIL  SUPPORT  UflLUE" 
2. 5-3.0" 
4.0-5.tf" 

2. 5-3.0" 


print  : 


*=;  ii 


2. 0" 
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3207  IF  01  =  1  GOTO  3945 

3218  PRINT  "TO  CONTINUE#  PRESS  SPACE  BAR" 

3220  POKE  -  16368,0: D I  =  PEEK  <  -  16384> 

3225  IF  01  <  >160  GOTO  3220 

3230  HOME 

3235  INPUT  "DO  VOU  WANT  THIS  LIST  IN  HARD  COPY  <0=NO,  1=VES>?"#D1 

3240  IF  D1  =0  GOTO  3950 

3935  L*  =  CHR$  C4>:  PRINT  L$,"PR#1" 

3940  GOTO  3101 


3945 

PRINT 

L$;"PR#0" 

3950 

01  =  0 

3955 

RETURN 

4000 

HOME 

4101 

PRINT 

TABC 

6 >" RECOMMENDED  REGIONAL  FACTORS" 

4102 

PRINT 

TABC 

15 >"FOR  GEORGIA" 

- 

4103 

PRINT 

TABC 

6  > " ***************************** " 

4104 

PRINT 

4105 

PRINT 

T  ABC 

10 >" AREA  1 

FACTOR" 

4106 

PRINT 

4107 

PRINT 

TABC 

6 >" COASTAL  PLAINS 

1.4-1. 7" 

4108 

PRINT 

TABC 

7  )"-SAUANNAH 

1.7” 

4189 

PRINT 

4110 

PRINT 

tabc 

6 >" PIEDMONT 

1.5-1. 3" 

4111 

PRINT 

TABC 

7  )" -ATLANTA 

1.8" 

4112 

PRINT 

TABC 

7  > " -MACON 

1.6" 

4113 

PRINT 

TfiGC 

7  >" -COLUMBUS 

1.8" 

4114 

PRINT 

TAB':: 

7  >" -AUGUSTA 

1.5" 

4115 

PRINT 

4116 

PRINT 

TABC 

6  >"UALLEY  &  RIOGE 

2. 0-2. 2" 

4400 

PRINT 

:  PRINT  :  PRINT 

4410 

PRINT 

:  PRINT 

4415 

IF  El 

=  1  GOTO  4945 

4420 

PRINT 

"TO  CONTINUE,  PRESS  SPACE 

BAR" 

4430 

POKE 

-  163 

68,0: El  =  PEEK  C  - 

1 6384  > 

4440 

IF  El 

<  > 

ISO  GOTO  4430 

4445 

HOME 

4450 

INPUT 

"DO  VOU  NANT  THIS  LIST  IN 

HARD  COPY  <0=1 

4455 

IF  El 

=  0  G 

OTO  4950 

4460 

L$  =  CHR$  c 

4 >:  PRINT  L$j"PR#t" 

4465 

GOTO  4101 

4945 

PRINT 

L*;"PR#0" 

4950 

El  =  0 

4955 

RETURN 

5000 

HOME 

5101 

PRINT 

"  EOIJIUALENT  SINGLE  AXLE  13-KIP  LOADS' 

5102 

PRINT 

TABC 

11 >"20  YEAR  DESIGN 

LIFE" 

5 1 0.' 

F'R  I  NT 

■*~******************************** 

5104 

PRINT 

5105 

PRINT 

"  CLASS  TYPE  PUNT 

EuUIUALENT" 

5106 

PRINT 

TABC 

23  >" 18-KIP  LOAD*" 

5107 

PRINT 

"  LIGHT  PARKING,  CITY" 

5 1 08 

PRINT 

TABC 

1 1 >" STREET ,  RURAL 

22000-428000 " 
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51Q9  PRINT  TABC  1 1 )" ROADS" 

5110  PRINT 

5111  PRINT  "  MEDIUM  SECONDARV  420000-3000000" 

5112  PRINT  TABC  11  )"HIGHHAY" 

5113'  PRINT 

5114  PRINT  "  HEflUV  INTERSTATE  3000000-10000000" 

5115  PRINT  TABC  1 1  )“ HIGH WAY" 

5116  PRINT  :  PRINT 

5117  PRINT 

5113  PRINT  TABC  6:)"*  BASED  ON  580.7  13-KIP  LOADS  PER" 

5119  PRINT  TABC  7 >"1000  TRUCKS" 

5120  PRINT  :  PRINT 

5125  IF  FI  =  1  GOTO  5945 

5130  PRINT  "CTO  CONTINUE#  PRESS  SPACE  BAR >" 

5140  POKE  -  16368#0:FI  =  PEEK  <  -  16384) 

5150  IF  FI  <  >  160  GOTO  5140 

5160  HOME 

5170  INPUT  "DO  VOU  HANT  THIS  LIST  IN  HARD  COPV  C0=NO,  1=VES)?  "#F1 

5180  IF  FI  =  0  GOTO  5950 

5190  Li  =  CHRi  <  4  >:  PRINT  Li;"PR#l" 


5200 

GOTO  5 

101 

5945 

PRINT 

Li# "PR#0" 

5950 

FI  =  0 

5955 

RETURN 

6000 

HOME 

6101 

PRINT 

TABC  5) "SELECTED  STRUCTURAL  COEFFICIENTS 

6102 

PRINT 

TABC  5  ) " -iHHHHUHE-  ************************* 

6103 

PRINT 

6104 

PRINT 

II 

PAUEMENT  COMPONENT 

_) 

cc 

i — 

to 

COEFF 

II 

• 

6105 

PRINT 

TABC  22) "FLA  GA  MD 

SC" 

6 1  OG¬ 

PRINT 

ll 

SURFACE  COURSE" 

GI  07 

PRINT 

II 

ASPHALT  CONCRETE  .2-. 

4  •  44* 

.44 

.44" 

6103 

PRINT 

6 1 09 

PRINT 

II 

BASE  COURSE" 

6110 

PRINT 

II 

ASPHALT  CONCRETE  .2-. 

3  . 30 

•  lLO 

.34" 

6111 

PRINT 

II 

SAND-ASPHALT 

.  12 

.23 

.25" 

6112 

PRINT 

11 

SOIL-CEMENT  .22 

.20 

o 

• 

.20" 

61 13 

PRINT 

ll 

GRADED  AGGREGATE  - 

.  13 

.  14 

.  15" 

6114 

PRINT 

II 

CEMENT  STABILIZED  - 

.22 

.  28 

.34" 

6115 

PRINT 

II 

GRADED  AGGREGATE" 

6116 

PRINT 

6117 

PRINT 

It 

SUBBASE " 

6113 

PRINT 

II 

GRADED  AGGREGATE 

.  14 

- 

_  ll 

6119 

PRINT 

II 

TOPSOIL  OR  SAND- 

.  10 

- 

_  II 

6120 

'PRINT 

ll 

CLAV" 

6122 

PRINT 

TABC  8)"*  MAXIMUM  DEPTH 

OF  4.5 

IN.  " 

6200 

PRINT 

6202 

IF  G1 

= 

1  GOTO  6945 

6205 

PRINT 

II 

CTO  CONTINUE#  PRESS  SPACE  BAR) 

ll 

6210 

POKE 

- 

16368 #0:61  =  PEEK  C  - 

16384) 

6215 

IF  GI 

<’* 

>  160  GOTO  6210 

6220  HOME 

SS  Til  ™  SOTot950LIST  IN  HflRD  C0PV  <0=N0'  1=VES>?“iS‘ 
6935  L$  =  CHRf  <4>:  PRINT  L*i"PR#l« 

6940  GOTO  6101 
6945  PRINT  L$;“PR#0 11 
6950  61  =  0 
6955  RETURN 
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5.6  Variable  List  (AASHTO  1) 


Input 


H$  =  Heading 

Z  =  #  of  trials 


TL  =  Tolerance 


SI  =  Soil  support  value 

PT  =  Terminal  serviceability  index 

RI  =  Regional  factor 

WI  =  Equivalent  18-kip  single  axle  loads 
L  =  Layer  #  for  calculation 
Y  =  Layer  # 

A  ( Y)  =  Structural  coefficient 
D  ( Y)  =  Layer  thickness 
N  =  #  of  layers 
D (L)  =  Layer  thickness 
Flow  Control 

BB  =  Calculate  thickness  of  E-18's 
XX  =  Calculate  thickness  of  E-18's 
YY  =  Calculate  thickness  of  E-18's 
ZZ  =  Calculate  thickness  of  E-18's 
Q1  =  Table  listing 
Cl  =  Table  listing 
Dl  =  Table  listing 
El  =  Table  listing 
Fl  =  Table  listing 
G1  =  Table  listing 


rft 
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Counters 
H  =  Run  # 

I  =  Layer  # 

Miscellaneous 

X  =  Structural  number  +  1 

AA  =  Intermediate  equation  values 

AB  =  Intermediate  equation  values 

NN  =  Rounded  off  structural  number 

SN  =  Structural  number 

SC  =  Structural  number  (with  N-l  layers) 
A ( 1 ) -A ( 5 )  =  Intermediate  equation  values 
B  (1 ) -B ( 7 )  =  Intermediate  equation  values 


«  _  V  .SjbVaVj 
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5 . 7  Program  Verification 


3  RUN 

*  AASHTO  * 


DANA  K.  EDDV ,  578-30-8378 
GA.  INSTITUTE  OF  TECHHOLOGV 
SCHOOL  OF  CIUIL  ENGINEERING 
DEPORTMENT  OF  GEOTECHNICAL  ENGINEERING 


SYSTEM  HARDWARE:  APPLE  II  PLUS  <6410 

SYSTEM  HARDWARE:  DOS  3.3,  APPLESOFT  BASIC  LANGUAGE 

PROGRAM  DATE:  JULY,  1983 


HASH TO  IS  BASED  ON  THE  PASHTO  INTERIM  GUIDE  FOR  THE  DESIGN  OF  FLEXIBLE  PAGE ME NT 
STRUCTURES,  1973. 


PRESS  THE  SPACE  BA  TO  CONTINUE) 

uSER  INPUT  QUESTIONS  AND  COMMANDS  WHICH  ARE  FOLLOWED  BY  <*>  ARE  SUPPLEMENTED  BY 
LISTS  AND  TABLES.  TO  ACCESS  THE  LIST  OR  TABLE  TYPE  <0>  INSTEAD  OF  THL  APPROPRIh 
TE  UALUE  OF  THE  VARIABLE  REQUESTED. 


PROBLEM  HEADING  AN  EXAMPLE  PROBLEM 


TO  CALCULATE  THE  THICKNESS  OF  A  PARTICULAR  LAYER,  TYPE  < 1 >; 
I MUM  NUMBER  OF  EQUIVALENT  18-KIP  LOADS,  TYPE  <2  ■.  1 


TO  CALCULATE  THt 


MAX 


HOW  MANY  TRIALS  DO  YOU  WANT  TO  RUN?  1 

WHAT  IS  THE  CALCULATION  TOLERANCE  (REC'D;  .01).  .01 
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SOIL  SUPPORT  UfiLUE  <*>  =  6 


SOIL  SUPPORT  UfiLUE 
<  fifiSHO  > 

-*•■*■■*-*■ ■*  ■$•**•*••*•■*■*•*■■*■£ 


SOIL  DVNANIC  S 

TflTIC  flflSHO 

MODULI! 

IJPPORT  CBR 

UfiLUE 

CBR  3  PT. 

MR 

PS  I 

S - ™  GO  - 

78  -  14.5 

—  Nfl 

7 - 35  - 

38  -  9.75 

—  Nfi 

6 - 17  - 

19  -  8.75 

—  930O 

5 -  8  - 

11  -  4.5 

—  6400 

4 -  5  - 

6  -  2.5 

—  4400 

3 -  3  - 

2.3  -  1.25 

-  3000 

2 - 1.5 - 

1.3  -  O.50 

—  2100 

1 - 0.5 - 

0.5  -  0.25 

—  Nfi 

(TO  CONTINUE  LIST ,  PRESS  SPfiCE  BOR)  _ 

SOIL  SUPPORT  UfiLUE  * 

< GO  DOT  > 

«r  ■*  *  ■*  *«*■• *  #•  *  ■*■  *  -Jr  #■  *  ** •fr 

/  '  * 

REGION  SOIL  SUPPORT  UfiLUE 

PIEDMONT  2. 5-3.0  * 

COftSTfiL  PLAIN  4. 0-5.O 

UflLLEV  &  RIDGE  2. 5-3.0 

I 


TO  CONTINUE#  PRESS  SPfiCE  BAR 

00  VOU  WANT  THIS  LIST  IN  HARD  COPV  (0=N0,  1=VES- 
SOIL  SUPPORT  UfiLUE  ( * >  =  3 


,'RD-A129  271 


UNCLASSIFIED 


COMPUTER  APPLICATIONS  TO  GEOTECHNICAL  ENGINEERINGS) 
AIR  FORCE  INST  OF  TECH  WRIGHT-PATTERSON  AFB  OH 
D  PC  EDDV  AUG  82  AFIT7CI/NR-82-8ST 


s  / 


TERMINAL  SERICEABILITY  INDEX  <*>  =  0 


TERMINAL  SERVICEABILITY  INDEX#  PT 
< AASHTO > 

*************************************** 


CLASSIFICATION 

PT 

PRIME  ROUTES,  MAJOR 

ARTER I ALS , EXPRESSWAYS 

2.5 

PRIME  SECONDARY  ROUTES, 

IND.  &  COMM.  STREETS 

2.25 

MINOR  SECONDARY  ROUTES, 
RESIDENTIAL  STREETS, 
PARKING  LOTS 

2.0 

FAILURE,  DEFINED  BY 

AASHTO 

1.5 

TO  CONTINUE,  PRESS  SPACE  BAR 

DO  YOU  WANT  THIS  LIST  IN  HARD  COPY  <©=NQ, 

TERMINAL  3ERICEABILITY  INDEX  <*>  =  2.5 


REGIONAL  FACTOR  <*>  =  0 


RECOMMENDED  REGIONAL  FACTORS 
FOR  GEORGIA 


************ 


AREA 

FACTOR 

COASTAL  PLAINS 

1.4-1. 7 

-SAUANNAH 

1.7 

PIEOMONT 

1.5-1. 8 

-ATLANTA 

1.8 

-MACON 

1.6 

-COLUMBUS 

1.8 

-AUGUSTA 

1.5 

UALLEY  &  RIDGE 

2. 0-2. 2 

i 


Is  YE 
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TO  CONTINUE#  PRESS  SPACE  BAR 

DO  YOU  WANT  THIS  LIST  IN  HARD  COPY  <0=NO#  1=YES>?  6 

REGIONAL  FACTOR  <*)  =  1.5 


»  OF  EQUIUALENT  SINGLE  AXLE  18-KIP  WHEEL  LOADS  C*)  =  O 


EQUIUALENT  SINGLE  AXLE  1S-KIP  LOADS 
20  YEAR  DESIGN  LIFE 
*********************************** 


CLASS  TYPE  PUMT  EQUIUALENT 

18-KIP  LOAD* 

LIGHT  PARKING#  CITY 

STREET,  RURAL  22000-420000 
ROADS 


MEDIUM  SECONDARY 
HIGHWAY 


420000—3000000 


HEAUY  INTERSTATE  3000000-10000000 
HIGHWAY 


*  BASED  ON  580.7  13-KIP  LOADS  PER 
1000  TRUCKS 


CTO  CONTINUE#  PRESS  SPACE  BAR) 

DO  YOU  WANT  THIS  LIST  IN  HARD  COPY  <0*NO#  1*YES>?  9 

»  OF  EQUIUALENT  SINGLE  AXLE  18-KIP  WHEEL  LOADS  <*>  =  8590000 


TRIAL  #1 


STRUCTURAL  NUMBER  =5.71 


SPECIFY  #  OF  LAVERS  C MAX- 18).  3 


LAVER  #  FOR  THICKNESS  CALCULATION.  1 


******************** 


LOVER  #.  1 

STRUCTURAL  COEFFICIENT  C*>  =  8 

SELECTED  STRUCTURAL  COEFFICIENTS 

PAUEMENT  COMPONENT  STR'L  COEFF. 


FLA 

GA 

MD 

SC 

SURFACE  COURSE 

ASPHALT  CONCRETE  .2-.  4 

.44* 

.44 

.44 

BASE  COURSE 

ASPHALT  CONCRETE  .2-.  3 

.30 

.28 

.34 

SAND-ASPHALT 

.12 

.28 

.25 

SOIL-CEMENT  .22 

.20 

.28 

.20 

GRADED  AGGREGATE 

.18 

.14 

.15 

CEMENT  STABILIZED  - 

.22 

.28 

.  34 

GRADED  AGGREGATE 
SUBBASE 

GRADED  AGGREGATE  -  .14 

TOPSOIL  OR  SAND-  -  .10 

CLAV 

*  MAXIMUM  DEPTH  OF  4.5  IN. 

<  TO  CONTINUE,  PRESS  SPACE  BAR ) 

00  VOU  WANT  LIST  IN  HARD  COPV  <0=NO,  1=YES>?  0 
STRUCTURAL  COEFFICIENT  <*>  =  .44 


it ir * 

LAVER  #.  2 


STRUCTURAL  COEFFICIENT  <*>  =  .14 


LAVER  THICKNESS  < INCHES)  =  12 


•Hr  * <8*  &  ■£ *r 4? & •& *8“ & ■£  *£* &  & ■£ 

LAVER  #.  3 


STRUCTURAL  COEFFICIENT  <*>  =  .11 


LAVER  THICKNESS  CINCHES)  =  13.5 


AH  EXAMPLE  PROBLEM 

V  V V V  W V V vvv V V V V V  V V V V V V V V  V 


TRIAL  #1 


LAYER 

STR'L  COEFF. 

THICKNESS 

1 

.44 

5.8 

2 

.  14 

12 

.11 

13.5 

STRUCTURAL  HUMBER  =5.71399934 


EQUIUALENT  SINGLE  AXLE  18-KIP  LOAD  =8569000 


SOIL  SUPPORT  UALUE  =3 

TERMINAL  SERUICEABILITV  INDEX  =2.5 

REGIONAL  FACTOR  =1.5 

DO  VOU  WANT  TO  RUN  ANV  MORE  TRIALS  'C0=NO,  1=YES>?  1 

FOR  A  LAYER  THICKNESS  CALCULATION,  TYPE< 1 >#  FOR  AN  18-KIP  LOADING  CALCULATION, 
YPE'-'  2  >.  2 

HON  MANY  TRIALS  DO  VOU  WANT  TO  RUN?  1 

SOIL  SUPPORT  UALUE  <*>  =  3 

TERMINAL  SERICEABILITV  INDEX  <*)  =  2.5 


REGIONAL  FACTOR  <*>  =  1.5 


TRIAL  #i 


SPECIFY  #  OF  LAVERS  (MAX-10)  3 

it -£■¥: if  it 

LAVER  #  =1 

STRUCTURAL  COEFFICIENT  (*>  =.4 

LAVER  THICKNESS  < INCHES )  =  5.8 

LAVER  #  =2 

STRUCTURAL  COEFFICIENT  (*>  =. 1 

LAVER  THICKNESS  (INCHES)  =  12 

■fr  &  4? -ar  && ie-ie & & ■£ & & & & ■£  & 

LAVER  #  =3 

STRUCTURAL  COEFFICIENT  <*>  =. 1 


LAVER  THICKNESS  (INCHES)  =  13 


■>  ■?  ■>  ■ 


-r .  f  r  .  ■  r  l.*1.' 


IV  Tr^  '  (T" 


wvvvwvw  vwwwwwww 

AN  EXAMPLE  PROBLEM 


TRIAL  #1 

******** 


LAVER 

STR'L  COEFF. 

THICKNESS 

1 

.44 

5.8 

2 

.  14 

12 

3 

.11 

13.5 

STRUCTURAL  NUMBER  =5.717 

EGUIUALENT  SINGLE  AXLE  18-KIP  LOADS  =8421531 


SOIL  SUPPORT  UALUE  =3 

TERMINAL  SERUICEABILITV  INDEX  =2.5 

REGIONAL  FACTOR  =1.5 

DO  VOU  WANT  TO  RUN  ANV  MORE  TRIALS  <0=N0,  1=VES>?  8 


THANK  VOU  FOR  USING  AASHTO 


*  ~  *  \  *  V* 


BVE-BVE 


Lift 

UK  «Z 

Lift  *  3 

Soil.  SuPPof^T  VALUE  =  3  (  5_/ ) 

7 ~£KrYiiHflt_  Sb(1\J\CE  flriu.  tri  lAJb£?c  =  2  .  S' 

Becj/o'JtL  Fflcroz.  -  /.S'  (ft) 

E  -/S'j.  =  8,  sco,  ooo  (E)t) 


Tz y  Sa/  =  s.q 


J2o^  (kj-t)  =  f-36  ^  ( S>/)  -.21- 


4.2  -  IS 


Z.~i 


/094 


J&Q 


'3  ' 


■4  f 

^  r  .372  (3-  3-o) 


~  .  4  , 

=  /(,. 77  -.2,  ^ 

<  SOOO^f 


-  /S\Z<4:T 

-  ^ir7,  2476  (E-/8i) 


7>y  Sa/  -1.0 

~^°o  ~  0^*0 


z  + 


-y>6 

<  4  1- 


4.Z-Z.S-\ 
2.^  y1 
/  044 
(1H)S-'C 


J_ 

/.S' 


,  S72  (3  -  ?.  0) 


--  74,46, 


Z  + 


-,46, 

■  4zzs 


- ,  4cr 


n  ,i(,q 

fhiii.i'n 


[£  -/S  ,) 


r~» 


7/z.V  SN--  6 . 0 


j2o£  (U-t)  -  Zb  (lc  ri)  -,  Z  -f - 


/Q04 
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=  /S’,  732,  311  (E-/%'s) 


Tky  Sa/  =  ■5'.  5" 
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Instructions  for  Using  AASHTO  Interim  Guides  for  Pavement  Structure  Design 
in  Georgia  -  81ank  form  on  page  number  6 


1.  Project  Number 

\ 

2.  County 

* 

3.  Description:  Describe  the  project  giving  length,  termini,  number  of 

lanes,  new  or  existing  location,  widening  and  resurfacing,  overlay, 
or  any  other  significant  information. 

4.  Type  of  Adjoining  Pavement:  Metal  surface,  bituminous  surface  treat¬ 

ment,  asphaltic  concrete,  or  P.C.  concrete. 

5.  Traffic  Data:  Derive  the  mean  traffic  in  VPD  for  one  uirection  during 

the  design  period.  The  design  period  for  Interstate  and  Primary 
projects  is  20  years  and  for  Secondary  and  other  projects  the 
design  period  is  15  years. 

6.  Design  Loading:  Distribute  the  mean  AAC7  for  one  direction  into  the 

highest  design  lane  traffic  with  a  division  between  trucks  ard 
other  vehicles.  Use  a  truck  classification  of  multiple  units  and 
single  units  if  available.  Use  estimated  lane  distribution  factors 
if  the  project  is  a  multi-lane  facility  (page  *7  ). 

Multiply  the  number  of  vehicles  in  each  vehicle  classification  by 
an  appropriate  18  kip  single  axle  equivalent  load  (18^  S.A.E.L.) 
factor  (page  8  ).  The  load  factors  may  be  derived  as  required  for 
specific  cases;  i.e.,  a  road  leading  to  a  pulpwood  yard  or  a  quarry. 

Sum  the  daily  18  S.A.E.L.  and  multiply  by  the  number  of  days  in 
the  design  period  to  derive  the  total  design  period  loading. 

7.  Design  Data:  Set  the  design  terminal  serviceability  (P^)  at  2.5  for 

Class  III  or  higher  projects  and  2.0  for  Class  IV  or  less  projects. 
The  soil  support  value  is  furnished  in  the  soil  survey  but  if  a 
soil  survey  is  not  conducted  the  soil  support  value  may  be  estimated 
in  conjunction  with  the  district  soils  engineer  and  the  soil  labora¬ 
tory  at  Forest  Park.  The  regional  ""actor  is  also  given  in  the  soil 
survey  but  may  be  estimated  from  the  attached  chart  if  a  soil  survey 
is  not  conducted. 

8.  Recommended  Flexible  Pavement  Structure:  Enter  the  respective  monograph 

for  terminal  serviceability  of  2.0  or  2.5  with  the  soil  support 
value,  18^  single  axle  equivalent  loads  and  the  regional  factor  to 
find  the  required  weighted  structural  number  of  the  pavement  struc¬ 
ture  to  be  designed.  This  weighted  structural  number  is  then  matched 
by  a  design  structural  humber  which  satisfies  thefequation  SN  -  aj  di^ 
a2  ^2  +  a3  ^3  +  an  dn  where  a  =  structure  layer  coefficient  and  d  = 
depth  of  each  respective  layer  in  inches.  Engineering  judgment  ard 
knowledge  of  local  materials  must  be  used  to  set  depths  ard  types  of 
materials  in  the  pavement  structure;  the  district  construction  engi¬ 
neer,  district  materials  engineer,  or  the  state  bituminous  construc¬ 
tion  engineer  should  be  consulted  if  required. 
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COMMITTEE  ON  ROADWAY  PAVEMENT  STRUCTURES 

PROCEDURE  FOR  ESTABLISHING  PAVEMENT  STRUCTURE  DESIGNS 
FOR  SECONDARY  PROJECTS,  AUTHORITY  PROJECTS  AND  COUNTY  CONTRACT  PROJECTS 

Secondary  Projects  ’ 

> 

The  procedure  for  making,  processing,  and  reviewing  these  designs  shall 
be  as  follows: 

A.  Projects  Designed  in  District  Offices 

1.  The  Designs  will  generally  be  made  for  the  District  Engineer 
by  his  Design  personnel  in  cooperation  with  the  Assistant  Dis¬ 
trict  Engineer,  Pre-Construction  anc  the  District  Materials 
Engi neer. 

2.  The  District  Engineer  shall  transmit  this  design  to  the  State 
Road  and  Airport  Design  Engineer  listing  the  thickness  and 
materials  to  be  used  in  each  layer  of  the  pavement  structure 
down  through  the  subgrade  stabilizer  aggregate  or  select  ma¬ 
terial,  if  these  are  required.  The  District  Engineer  shall 
also  include  information  as  to  the  availability  of  local  ma¬ 
terials  used  in  the  design,  and  the  availability  of  alternate 
materials.  This  submission  shall  be  made  as  early  as  possible 
in  the  plan  development  process. 

3.  The  State  Road  and  Airport  Design  Engineer  shall  check  the  design, 
reconcile  any  differences  of  opinion  with  fhe  District  Engineer 
and  the  State  Materials  and  Research  Engineer,  and  if  the  road¬ 
way  is  Class  I  or  II,  he  shall  submit  the  design  to  the  Committee. 

If  lower  than  Class  II,  the  design  need  not  be  submitted  to  the 
Committee  if  concurred  in  by  the  State  Pond  and  Airport  Design 
Engineer  and  the  State  Materials  and  Research  Engines:*-. 


4.  After  consideration  by  the  Committee,  the  State  Road  and 
Airport  Design  Engineer  shall  inform  the  District  Engineer 
that  either  the  design  is  approved  or  that  the  Committee 
desires  specific  changes  in  the  design,  and  that  the  approved 

t 

or  modified  design  may  be  u^ed  for  the  project  plans. 

5.  If  the  District  Engineer  wishes  to  question  modifications 
made  by  the  Committee,  he  may  do  so  through  the  State  Road 
and  Airport  Design  Engineer. 

B.  Projects  Designed  in  the  General  Office 

1.  These  designs  shall  be  prepared  in  accord  with  policies  of 
the  Committee. 

2.  If  these  roadways  are  Class  I  or  Class  II,  they  shall  be 
submitted  to  the  Committee  in  regular  manner. 

3.  If  the  roadway  is  lower  than  Class  II,  the  design  need  not 
be  submitted  to  the  Committee  if  concurred  in  by  the  State 
Urban  and  Multi-Modal  Design  Engineer  or  the  State  Road  and 
Airport  Design  Engineer  (whichever  has  responsibility  for 
the  plans)  and  the  State  Materials  and  Research  Engineer. 

II.  Authority  Projects 

A.  When  requested  by  the  Director  of  Operation.,  the  Committee  will  ,-ev 

and  make  recommendations  concerning  Authority  Projects. 

III.  County  Contract  Projects 

A.  When  requested  by  the  Director  of  Operatic  ns,  ..he  committee  will  re; 

and  make  recommendations  concerning  county  contract  proje  r-lS. 


PAVEMENT  STRUCTURE  DESIGN  COMMITTEE 

Consists  of: 

State  Highway  Engineer  -  Chairman 

State  Materials  and  Research  Engineer  -  Secretary 

Director  of  Operations 

Director  of  Pre-Construction 

Director  of  Construction  and 
State  Construction  Engineer 

State  Road  and  Airport  Design  Engineer 

State  Urban  and  Multi-Modal  Design  Engineer 

State  Maintenance  Engineer 

Project  Review  Engineer 

FHWA  Representative 

Reviews:  All  projects  over  1000  VPD  (future  traffic)  ana  any 
others  requested. 

Design  Methods  Approved: 

(A)  Flexible  Pavement 

(a)  over  1000  VPD  -  AASHO  Interim  Guide 

(b)  under  1000  VPD  -  generally  use  ultimate  strengtn 

(B)  Rigid  Pavement  AASHO  Interim  Guide 

(C)  Overlay 

(1)  Flexible  over  flexible 

(a)  AASHO  structural  coefficients  assigned 

(2)  Flexible  over  rigid 

(a)  AASHO  structural  coefficients  assigned 

(b)  Corps  of  Engineers 

(c)  layered  elastic  theory 

(3)  Rigid  over  flexible  -  AASHO  Interim  Guide 

(4)  Pigid  over  rigid 

(a)  Corps  of  Engineers 

(b)  AASHO  Interim  Guide 


The  percent  over-under  design  should  be  computed  by  dividing 
the  difference  between  the  weighted  structural  number  and  the 
total  SN  and  dividing  by  the  weighted  structural  number.  The 
percent  over-under  design  should  be  limited  to  about  fifteen 
percent. 


Any  additional  pertinent  information  should  be  entered  under 
"remarks".  The  submissiOR  of  the  form  should  follow  T.  D. 
Moreland's  instructions  erf  January  6.  1970. 


Beginning  of: 
End  of  Projer 


> 

Type  of  Adjoining  Pavement: 


Traffic  Data:  2h  Hr.  Truck  Percentage 

AADT  Beginning  of  Design  Period 

v?n 

Year 

AAIZP  End  of  Design  Period 

V?  D 

"""  Year 

Mean  AADT  (One  Way) 

VPD 

Design  Loading: 

Design  Lone  Traffic 

lS*  Axle 

Zq.  Load 

Y  - 

x  ° 

x  z 

'  Total  Diily  Leading; 

Total  Design  Period  Loading  =  _ _ 


Design  Data;  Serviceability  (P^)  _ 

(Prom  Soil  Survey)  Soil  Support  Value  (S) 
Recommended  Flexible  Pavement  Structure: 


Regional  Fhctor  ( R) 


Type  of  Material 

Thickness 

Coefficient 

SN 

» 

I 

i 

_ 

_ ! _ 

Weighted  Structural  Value  (3N)  ( Prom  Nomograph)  - 
Actual  Design  Life  (Years) 

[Remarks : 


Tot  o.l  3N 


ereent  Cvrr-’Jr.der  Design 


Prepared  By: 


Bjbir.i  tied  By: 

Recommended 

.Lite 

District  tVy-Constructior  Lr.y it. 

L-jtO 

Approved 

State  Road  Decign  Eng  s.vr.-r 

Chairman,  Committee  on  load  way  Pavement  Si '■*ieT.u~.-. 

.  v.e 

Table  4.  Lane  Distribution  Factors 


FACILITY 


Four  lane  rural  freeway 

Four  lane  urban  freeway 

Six  lane  rural  freeway 

Six  lane  urban  freeway 

Four  lane  rural  highway-free  acces< 

Four  lane  urban  free  access 

Two  lane  highways  and  ramps 


Design  Lane 
Distribution  Factors 

-  percent  of  One  Way  Truck 
in  the  Heaviest  Lane 


100 


ICO 


18*  Single  Axle  Equivalent  Loads  for  Flexible  Pavement 


All  Multiple  Unit  Trucks 

Average  18k  S.A.E.L. 
.1.4 

5  Axle  M.U.  Trucks  Only  > 

1.7 

4  Axle  M.U.  Trucks  Only 

1.1 

3  Axle  M.U.  Trucks  Only 

0.9 

All  Single  Unit  Trucks 

0.4 

3  Axle  S.U.  Trucks  Only 

0.9 

2  Axle  S.U.  Trucks  Only 

0.2 

Facility 


%  M.U.  Trucks 


%  S.U.  Trucks 


Average  18*  S.A.E.L 
Flexible  Pavement 


Interstate 

Routes 

100 

90 

80 

0 

10 

20 

1.4 

1.3 

1.2 

Primary  System 

70 

30 

1.1 

Heavy  State  Routes 

60 

40 

1.0 

50 

50 

0.9 

Medium  State 

40 

60 

0.8 

Routes 

30 

70 

0.7 

20 

80 

0.6 

Light  State  Routes, 
Secondary  System, 
City  Streets 


10 

0 


90 

100 


0.5 

0.4 


STATE 


HIGHWAY  DEPAR 
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Figure  71-1  Design  Chart  for  Flexible  Pavements,  ] 
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STRUCTURAL  COEFFICIENTS  FOR 
PAVEMENT  DESIGN 


APPLICABLE  DEPTH  BELOW 
SURFACE,  INCHES 


Surface  Course1 

> 

Asphaltic  Concrete  Surfacing  and  Binder 
Base  Courses 
Asphaltic  Concrete 

Graded  Aggregate  and  Crushed  Limestone 
(Compacted  to  modified  density) 
Graded  Aggregate  and  Crushed  Limestone 
(Compacted  to  standard  density) 
Topsoil  or  Sand  Clay  Bases 
Topsoil  or  Sand  Clay  (stabilized  with 
(150  Ibs./sq.  yd.  x  6"  stabilizer 
aggregate) 

Cement  Stabilized  Graded  Aggregate 
Soil -Cement 
Sand  Bituminous 

Subbase  Courses 

Graded  Aggregate  or  Crushed  Limestone 
Topsoil  or  Sand  Clay 
Sandy  Gravel  (Float  materia!  which 
fails  graded  aggregate) 

Crushed  Aggregate  Subbase 

Subqrade  Courses 

Class  I  Soil 
Class  II  Soil 


Coef . 
Per  Inch 


Below  12 


FLEXIBLE  PAVEMEJT  DESIGN  ANALYSIS 

- - -  20  7 

(Based  on  AASHO  Interim  Guide  for  the  Design  of  Flexible  Pavement  Structures) 

Project:  County:  _  _ 


Description:  Two  lane 


Type  of  Adjoining  Pavement:  Beginning  of'Project:  asnha [tic  cone 


End  of  Project:  asphal tic  corcrett 

Traffic  Data:  24  Hr.  Truck  Percentage  _ iO^ _ 

AADT  Beginning  of  Design  Period  _ 400 _  VPD  _ 

AADT  End  of  Design  Period  _  900 _  VPD  _ 

Mean  AADT  (One  Way)  (400+900)  x  1  x  r  =  325  _ 


Design  Loading: 

Design  Lane  Traffic 

325  VPD  x  10%  Trucks  x  100%  Lane  D i st  =33  True l s 


l3k  Axle 
En.  Load 


Total  Design  Period  Loading  =  14  x  365  x  15  years 

Design  Data:  Serviceability  (P  )  2 . 0 

(From  Soil  Survey)  Soil  Support  Value  (s)  _ 3 . 0 

Recommended  Flexible  Pavement  Structure: 


Type  of  Material 

_ Asohaltic  Concrete  E _ 

_ Asphaltic  Concrete  A  or  9  _ _ 

_ Graded  Aggregate _ _ _  ... _ 


Total  Drily  Loading; 

m 52 _ 

Regional.  Factor  (?.) 


;  Thi  v.nr-.c  s 

I  '  ' , . 


Weighted  Structural  Value  (SN)  (From  Nomograph)  -  2 . 5 

Actual  Design  Life  (Years)  1 1 _  Percent  CtWT 

Remarks : 


Frepared  By 

Recommended 

Approved 

State  Road  Design  Frcinec 

Approved 

Chairman ,  Committee  or  Roadway 

.  >*;  t*1  ,  ■  /  .  i 

H38WITN  1V8fUdpHIS  031HDI3M 


Figure  1 1 -2  Design  Chnrl  for  Flexible  Piivomcnts 
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Instructions  For  Using  Ultimate  Strength  Pavement  Structure  Design  In  Georgia 


Personnel  other  than  those  assigned  to  the  Materials  arid  Testing  Laboratory 
are  not  authorized  to  design  pavement  structures  using  ultimate  strength  methods 
Ultimate  strength  is  used  when  a  thinner  pavement  section  is  obtained  than 
using  the  AASHTO  Interim  Guides.  This  normally  occurs  at  below  1000  VPD  with 
10%  trucks.  % 

The  method  is  based  on  using  a  design  criteria  of  preventing  shear  failure 
by  limiting  the  imposed  vertical  stress  on  the  subgrade  to  a  safe  limit  below 
the  bearing  capacity  of  the  subgrade. 

The  steps  in  design  area  as  follows: 

1)  Determine  total  loading. 

2)  Obtain  bearing  capacity  from  the  soils  laboratory. 

3)  Determine  the  safety  factor  from  page  ,  1')  . 

4)  Compute  the  safe  bearing  capacity  =  bearing  capacity  . 

safety  factor 

5)  Find  total  pavement  thickness  from  page  2 O  ■ 

6)  Use  engineering  judgment  and  experience  to  set  surface 
courses . 


Ultimate  Strength  Sample  Problem 
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1.  Total  Loading  -  Assume  same  as  page  13 

Total  Loading  =  7^.650  18*4  S.A.E.L. 

* 

2.  From  soil  survey  -  bearing  capaci-ty  of  subgrade  soil  =•  Bl  psi 

3.  Set  design  terminal  serviceability  =  2.0 

4.  From  page  19  -  safety  factor  =  1.6 

5.  Allowable  subgrade  bearing  capacity 

=  bearing  capacity  of  subgrade  soil  -r  safety  factor  -  3 1  psi  4-  !.C 

6.  From  page  20  -  a)  depth  required  for  stabilized  base  =  7?  inches 

use  minimum  design 

I3"  asphaltic  concrete  E 
2"  asphaltic  concrete  A  or  B 
6"  cement  stabilized  base 

b)  depth  required  for  non-s tab  I  I i zed  base  =  I3j  inches 
use  1^"  asphaltic  concrete  E 

2"  asphaltic  concrete  A  or  B 
10"  graded  aggregate 

c)  cement  treated  base  with  thin  surfacing  acceptable 
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CHAPTER  VI 


SUMMARY  AND  CONCLUSIONS 


6 . 1  Scope 

6.1.1  General 

This  report  presents  four  programs  designed  to  solve 
specific  engineering  problems.  The  report  organization 
is  based  on  the  theme  of  a  service  that  would  be  provided 
an  engineering  client  acquiring  engineering  software.  Each 
program  is  supported  with  background  theory,  programming 
rationale,  and  a  user's  guide;  additionally,  a  program  list, 
example  problems,  and  hand  verified  solutions  have  been 
included. 

The  programs  include  the  solution  of  the  embedded  post 
subject  to  lateral  loads  above  grade  (SIGNPOST  1) ,  the  solution 
of  the  cantilevered  sheet  pile  wall  (CANTWALL  1) ,  the  limit 
equilibrium  analysis  of  slope  stability  by  the  Bishop  method 
(BISHOP  1) ,  and  the  design  of  flexible  pavement  based  on 
the  AASHTO  Interim  Guide,  1972  ( AASHTO  1). 

6.1.2  Hardware 

Programming  was  performed  on  an  Apple  II-Plus  personal 
computer  with  64K  storage.  The  programming  language  was 
Applesoft  Basic  operating  on  DOS  3.3.  The  programs  were 


stored  on  5V  disks.  Peripheral  equipment  included  two 
disk  drives,  a  green  screen  monitor,  and  a  thermal  printer. 

6 . 2  Personal  Computers  in  Engineering  Practice 

Personal  computers  (PC)  have  become  quite  economical 
and  extremely  popular  in  the  last  five  years.  Businesses 
of  all  types  rely  upon  PC's  for  processing  routine  data, 
storage  and  retrieval,  word  processing,  cost  accounting, 
and  repetitive  problem  solving. 

When  General  Electric  introduced  the  first  main  frame 

computer  (MACH  I)  in  the  early  1950's,  computation  speed 

was  approximately  three  calculations  per  second.  MACH  I 

2 

required  1500  ft  of  floor  space,  a  considerable  air  condi¬ 
tioning  system,  and  constant  maintenance  replacing  vacuum 
tubes.  Today's  PC  requires  six  square  feet  of  desk  space, 
standard  environmental  controls,  and  a  minimal  maintenance 
program.  Computation  speeds  range  up  to  hundreds  of  calcu¬ 
lations  per  second.  Although  technological  advances  have 
drastically  increased  the  speed  of  main  frame  computing 
facilities,  the  size  remains  a  physical  constraint.  Many 
businesses  rely  upon  hard  wire  connections  through  telephone 
lines  to  utilize  main  frame  computers.  Costs  include  installa 
tion  and  maintenance  fees  as  well  as  charges  assessed  on 
compilation  and  computing  time. 

A  construction  cost  estimating  service  provided  by 
McGraw-Hill  Information  Systems  charges  $300  to  compile 


a  "quantity  take  off"  for  an  average  20,000  ft  ,  one  floor 
office  building.  The  cost  of  "time  sharing"  main  frame 
computing  facilities  can  be  large  and  subsequently  prohibi¬ 
tive  for  small  businesses. 

The  principal  advantage  of  the  PC  over  main  frame  time 
sharing  is  its  easy  access  and  low  relative  cost.  For  account 
ing  purposes,  the  purchase  of  a  PC  and  software  represent 
a  one-time  fixed  cost.  With  exception  of  maintenance  costs, 
the  value  of  a  PC  is  amortized  through  depreciation. 

The  absence  of  variable  costs  directly  influence  the 
cost  of  a  service;  consequently,  bid  prices  or  negotiated 
costs  can  be  lower.  Lower  service  costs  directly  influence 
the  volume  of  services  provided.  Firms  which  utilize  personal 
computers  can  maintain  a  sharp  competitive  edge. 

As  a  result  of  increased  personal  computer  use,  the 
demand  for  software  has  become  enormous.  In  response  to 
this  demand,  hundreds  of  small  software  service  companies 
have  been  formed.  This  is  particularly  evidenced  by  the 
advertisements  in  trade  magazines  and  professional  publica¬ 
tions  for  software  services.  At  this  time,  the  demand  appears 
to  be  almost  limitless. 

Many  firms  have  been  disappointed  with  the  services 
provided  when  purchasing  software  packages.  Although  the 
actual  programs  are  delivered  as  promised,  the  supporting 
documentation  has  left  the  user  deserted  and  basically  help¬ 
less.  It  is  with  this  in  mind  that  the  author  chose  the 


organization  and  theme  of  this  report.  The  bulk  of  this 
text  is  devoted  to  educating  the  user  about  the  programs 
use  and  limitations. 

6 . 3  Program  Applications 

This  section  presents  some  classic  applications  for 
which  the  programs  are  best  suited.  Iteration  of  the  back¬ 
ground  theory  is  avoided  as  each  program  is  supported  within 
its  respective  chapter.  Typical  applications  of  each  Program 
are  as  follows: 

6.3.1  SIGNPOST  1 

a)  Highway  signs  and  markers  subjected  to  wind  loads 

b)  Pole-type  buildings  which  resist  wind  loads  through 
embedded  post  columns 

c)  Commercial  signs  and  billboards  subjected  to  wind 
loads 

d)  Utility  poles  (power  and  telephone)  subjected  to 
cable  loads,  guy  wire  loads,  and  wind  loads 

6.3.2  CANTWALL  1 

a)  Shallow  excavation  (<15'-20')  when  surface  deflec¬ 
tions  are  tolerable 

b)  Marine  applications 

6.3.3  BISHOP  1 

a)  Earth  dams 

b)  Highway  cuts  and  fills 

c)  Slopes  near  or  under  structures 

d)  Any  slope  whose  failure  can  be  approximated  by  a 


circular  failure 


6.3.4  AASHTO  1 


a)  Flexible  pavement  design 

b)  Analysis  of  existing  pavement 

c)  Economic  feasibility  study 
6 . 4  Recommendations  for  Future  Work 

The  ultimate  goal  is  to  produce  a  program  capable  of 
solving  a  problem  given  an  endless  spectrum  of  varying  condi¬ 
tions  and  parameters.  Additionally,  a  program  should  be 
foolproof.  Experience  proves  these  goals  serve  as  sound 
guidance  but  are  impossible  to  totally  achieve.  Man-hour 
and  computer  memory  constraints  force  programming  efforts 
to  be  concentrated  around  specific  tasks  with  definable 
end  goals. 

As  with  any  program,  the  programs  presented  in  this 
report  can  be  improved.  The  limitations  listed  in  each 
chapter  can  best  serve  as  a  basis  for  improvement. 

6.4.1  SIGNPOST  1 

The  amelioration  of  this  program  might  include  a  routine 
to  compare  the  post  diameter  to  the  foundation  volume  (as 
a  function  of  the  required  depth) .  The  volume  of  a  wood 
pole  or  of  concrete  is  directly  related  to  cost.  The  user 
could  make  an  economic  judgment  as  to  the  best  diameter 
and  depth  of  a  foundation. 


6.4.2  CANTWALL  1 


This  program  can  be  enhanced  by  providing  the  ability 
to  specify  a  surcharge  above  the  upper  soil.  Soil  layering 
at  the  option  of  the  user  would  be  desirable.  A  major  improve 
ment  would  be  a  routine  to  calculate  the  bending  moment 
in  the  wall  considering  moment  distribution.  This  could 
be  complemented  with  a  section  modulus  selection. 

6.4.3  BISHOP  1 

Enrichment  of  BISHOP  1  could  include  a  search  routine 
in  which  the  minimum  factor  of  safety  is  calculated  without 
user  intervention.  Another  improvement  would  increase  the 
maximum  number  of  points,  lines,  and  soil  types  the  user 
can  specify. 

6.4.4  AASHTO  1 

AASHTO  1  could  be  ameliorated  by  providing  two  types 
of  economic  analysis.  The  first  anlysis  would  contrast 
the  use  of  alternate  construction  materials  on  a  unit  cost 
basis  such  as  dollars  per  square  yard  of  pavement.  The 
second  option  would  provide  a  life  cycle  cost  analysis. 


